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Die 84-Minuten-Abstimmung beim Kreisel- und beim Raumkompaf 
Von P. Christoph 


1. Einleitung. Fiir die Navigation ist es seit langem das letzte Ziel, nicht nur eine saubere 
Nordrichtungsanzeige an Bord zur Verfiigung zu haben, sondern méglichst auch eine zuverlassige 
Fahrt-iiber-Grund-Anzeige. Aus der wahren Fahrtgeschwindigkeit und dem Kurs des Fahr- 
zeuges kénnte dann in einem Rechengerat der jeweilige Standort des Fahrzeuges nach geographi- 
scher Breite und Lange gewonnen werden, wenn die Daten der Position am Ausgangspunkt der 
Reise gegeben sind. Es ist zwar méglich, die Fahrtgeschwindigkeit als Zeitintegral aus der ge- 
messenen Beschleunigung zu erhalten. Um jedoch Schwerebeschleunigung und Fahrtbeschleuni- 
ging voneinander sauber trennen zu kénnen, muB vor allen Dingen eine einwandfrei horizontierte 
Plattform als Trager fiir die BeschleunigungsmeBorgane vorhanden sein; die Plattform mu 
méglichst genau parallel zur wahren Horizontalebene liegen. Die drei Probleme, das der Fahrt 
tiber Grund, das des kiinstlichen Horizontes und das der wahren Nordrichtung, hangen eng mit- 
einander zusammen und bilden insgesamt das Problem der Inertial-Navigation. Es soll gezeigt 
werden, dafs man prinzipiell mit einem einzigen Kreisel die kreiseltechnische Aufgabe lésen kann. 


2. Die Bewegungsvorgiange auf der Erde infolge Erddrehung und Fahrtgeschwindigkeit. Zunichst 
sollen die kinematischen Vorginge betrachtet werden, wie sie bei beliebiger Bewegung eines Punktes 
auf der kugelférmig angesetzten Erdoberflache auftreten. Bezeichnet Si den vom Erdmittelpunkt 
ausgehenden Radiusvektor eines Erdoberflachenpunktes P, so ist seine durch den Erddreh- 
geschwindigkeitsvektor ll bedingte Geschwindigkeit im Raum gegeben durch 


dh 
a UX. 


Bewegt sich der Punkt P noch zusatzlich gegeniiber der Erde auf der Erdoberflache mit der Relativ- 
geschwindigkeit %, so gilt 
Hm _Uxh +8. (1) 
t 
Fiihrt man ein erdbezogenes rechtshandiges, orthogonales Dreibein i, j, f ein (Einheitsvektoren 
in der Nord-, Ost- und Nadirrichtung!, so kann man gemab 
KR =— RE (2) 
auch schreiben 
df ¥ B 
yo Uixt—_ (3) 
mit ; - 
eee | (4) 
U=i(-i)+(U- f =U (icosp — fsing) 
(U Betrag der Erddrehgeschwindigkeit, p geographische Breite, R Erdhalbmesser, x Kurswinkel.) 
Wegen f- f= 1 ist 


f-— =0, (5) 


df/dt liegt also wie B und WU x f in der Horizontalebene. Um di/dt und dj/dt durch 1, j, f, Wound & 
darzustellen, geht man von den Gleichungen aus 


t=O, = 05. 4°{=1, text, l= const. (6) 


1 Der Einheitsvektor £ hat die Richtung der reinen Gravitationskraft ausschlieBlich der Fliechkraft in- 
folge der Erdrotation. 
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Durch Differentiation ergibt sich 


dba hd Saat dj _ 
eee ar = poe dae | ae ee : 
. dj dimes! Melee, (7) 
Lp aoa da adh as ds 


Der Vektor dj/dt laBt sich durch die drei aufeinander senkrecht stehenden Vektoren U, j, U x j 


darstellen: 
t=4U + ej +e, xj. 


Skalare Multiplikation mit U, j und f ergibt nach den Gleichungen (7) 


C= 05 €, = 0 5 At aa REE 
Zur Berechnung von c, benutzt man die erste Gleichung aus (7), die mit (3) ergibt 


dt decid! eae tec 
pod j-%= j-(uxt R) 


Fir den Koeffizienten c, ergibt sich schlieBlich 


dj ; ¥ pe ae Me 
Ro j (tx te Weex poh wits 
Ae Re Wei Ui Re ies 
Somit ist eee: 
ext ana 
7 =(l teat j)Wxi= (+R )xi-i Rui 3 
Rectan Ols ae pie (8) 
=(u4 Ra SSR )xi- 
Aus der goeeuien Gleichung (7) lat sich mit (3) und (8) di/dt berechnen: 
Soh ae) RF ele AG. Boe 
g(t tg )x i+i SS =(U+-R pe )xi. (9) 
Nunmehr kann der Geschwindigkeitsvektor eines jeden erdfesten Vektors © = e, 1 + ej + eg f, 
mit 
dey de,” . de, _ 
eo: apa = 


angegeben werden. Es ist 


dé edi =u ee te OE] x6. 


Pre miele nil ap te °F 


I R IT 

Andererseits dreht sich ein inertialraumfester Vektor 3 = i,i + i,j + i, f im Erdsystem des 
Beobachters mit der Geschwindigkeit? 

5 fo di, _ ae (UL + 8) (B -}) d 

en ain ae 3x [4 rE RU-i |. ic ak 
Allgemein gilt fiir den Differentialquotienten eines beliebigen Vektors Q{ 

dX d : ; ; a a .d d di dj df 

@ ait hit) Sete ee tena a 

(Uf) (B +f) 
a + lu 4 oc oh. (10) 


Der Vektor 


ur =U+=Re 4 1 OPE (11) 


1 Mit einem Punkt sei der zeitliche Differentialquotient fiir Veranderungen relativ zur Erde bezeichnet! 
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ist der momentane Gesamtdrehgeschwindigkeitsvektor — von den Drehbewegungen des Schlingerns 
und Stampfens abgesehen —, den das auf der Erdoberflache befindliche Fahrzeug infolge Erd- 
drehung und EKigengeschwindigkeit im Inertialraum besitzt. 1(* laBt sich allein durch die Einheits- 
vektoren und ihre Ableitungen ausdriicken: 


di la. dt fk dh\. 
u = Exe + Ei ar) ee (i aE (12) 


dabei ist £ x df/dt die in der Horizontalebene liegende Komponente von l\* und f (j-di/dt) die 
Vertikalkomponente von \|*. 


3. Die Bewegungsgleichungen des Kreisels und Kompafkreisels auf der Erde. Ist der betrachtete 
Vektor der Drehimpulsvektor % eines erdfest gehaltenen Kreisels!, so andert sich seine Lage im 
Inertialraum gemaB 


G5 tre ap 
Gt x8 


Soll der Kreisel seine erdfeste Orientierung behalten, so mu man auf ihn das Drehmoment 
Mt = U* x B ausiiben. Dabei erzeugt er auf seine Umgebung das entgegengesetzt gerichtete 
Drehmoment der Kreiselwirkung, das verallgemeinerte Erdachsenrichtmoment § x 11*. Zeigt 
z. B. ein Kreisel die Vertikale an: 8 = B f, so muB man auf ihn das Drehmoment 


SO Sui (GL oR NES oh : B 
M = B(U ae. ¢ Rit) X $= —B(iV eos + 


einwirken lassen, damit er weiterhin die Vertikale anzeigt. AuBer der Fahrtgeschwindigkeit und 
der geographischen Breite muf also fiir die richtige Drehmomentenbildung auch noch die Ost- 
richtung j bekannt sein. 
Allgemein kann man die raumliche Veranderung des Kreiseldrehimpulses 8 gemaB (10) aus- 
driicken 
= S+UxB 7 (13) 
dt d 
Die Bewegungsgleichung des ungediémpft schwingenden KompaBkreisels ist durch den Dreh- 
impulssatz gegeben: 


dB * df aie seh S : df 
Babel karledpoa-Ealtereg) 


Hierbei ist G der vom Unterstiitzungspunkt des Kreiselsystems ausgehende Radiusvektor des 
Kreiselsystem-Schwerpunktes, multipliziert mit dem Gewicht des Systems, und g die Schwere- 


2 
beschleunigung. Es sei.S = |G]. Ferner ist — ne die Beschleunigung des Kreiselsystems im 
Inertialraum: 
ast PE ple dx atte hd ¥ 
Th ee ee Cae Re (ux t— | 
dt 1 dB\_ oh On een 
=— R(X — Rg) = —RUx (xt R) +B +UxB 
; 2 —(L- B)(B ei g 
= rit xs F2UXB+B+|BXt co =| ao 
Radar 
Endlich ist © x £ das durch die Erdschwere entstehende Drehmoment, © x Sia das durch 


die Beschleunigung entstehende Moment. 6-@ = 0 sei durch Konstruktion und Justierung 
erfillt. , 


4. Die beschleunigungsunabhingige azimutale und horizontale Gleichgewichtslage des Kompag- 
kreisels. Eine notwendige Bedingung dafiir, daf} die stationare Gleichgewichtslage des Kreisel- 
drehimpulses erreicht ist, ist das Verschwinden des verallgemeinerten Meridianrichtmomentes, der 
Vertikalkomponente des verallgemeinerten Erdachsenrichtmomentes 


t: (Be xt) =0 (15) 


1 Es werden nur schnelle Kreisel betrachtet (vgl. z. B. R. Grammel, Der Kreisel, 2. Aufl., Bd. 1, Berlin/ 
Gottingen/Heidelberg 1950). 
16* 
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oder auch 
Me Pe eet ine dea 
t [Box (EX =F ti z))| =f (Be x (§xz)|=2- 
Mit dem Entwicklungssatz fiir das zweifache Vektorprodukt ergibt sich weiter und wegen (5) 
a be Chi GN ope dt 
t-(E (Be: G)— gBo" D) = Beg = 0 (16) 


Der Vektor 8 mu8 also in der Gleichgewichtslage Bg auf dem in der Horizontalebene liegenden 
Vektor df/dt senkrecht stehen. Ergainzt man noch die notwendige Bedingung (16) durch den An- 
satz Bo: f= 0, d.h. Bg liegt in der Horizontalebene, so kann man insgesamt fiir die Gleich- 
gewichtslage Bg setzen ca 
Be ae A at c (17) 
Fir © folgt daraus und wegen §-© = 0 

Se == S c . (17a) 


In der Tat ist der Ansatz (17) und (17a) erfiillbar, wenn (14) erfiillt werden kann. Es miiBte somit 
d2 2 2 
Se Age es ee 


ay  g dt? g di? * 


gemaB (14) sein 


Diese Gleichung kann erfiillt werden, wenn 4 = S R/g gemacht wird. Sorgt man also dafiir, daB 
der Drehimpuls des Kompafkreisels stets durch die Gleichung 


SR df S R |dt SR BR 
De Wigs eee atl 5S ae 
B=|8|=[txG == z [ilxt R (18) 
gegeben ist, d. h. proportional zur Gesamtgeschwindigkeit des Fahrzeugs im Raum 4 =R | ue 
t t 
gehalten wird, und ist die momentane Richtung seines Vektors durch 
SR dt SR, ; 2h SIR fe Aa A V 
yas =, so Penney ee era: . ‘ 
Cee 2 1 z ix(Uxt ae - [i Veosp +i p sin x igo): (19) 
die des Schwerpunktvektors durch 
Ge =15 f (19a) 


gegeben, so ist B, und Gg eine stationdre Gleichgewichtslage. Wie der Aufbau der Gleichung (19) 
2 
zeigt, ist sie véllig unabhaingig von Beschleunigungen Ree - Die aus (19) folgende zeitliche 


Anderung von Bg 


dBe_ d (SR yy dt) _ SRy at 
de dt\ g Gey ace ae 


wird exakt durch das Beschleunigungsdrehmoment 
Rd R df SRY 
Sex (Et; ga) = SEX(E ; a) ar ix 


erzielt. Die zeitliche Anderung von Wg relativ zur Erde ergibt sich aus (13) zu 


. aS ; Ss : 
Bet (UA eee 
c tu toe) : te 


Sie liegt in der Horizontalebene. Befindet sich der Drehimpulsvektor % nicht in der Gleichgewichts- 
lage Be, so filhrt er Schwingungen um diese aus, die bei kleinen Auslenkungen bekanntlich har- 
monisch sind. Nach Kinfiihrung Eulerscher Winkel kann man (14) im (i, j, )-System koordinaten- 
mafig darstellen. Die so entstehenden Koordinatengleichungen sind ziemlich uniibersichtlich und 
sollen deshalb nicht angeschrieben werden. Fiir verschwindende F ahrt V = 0 und bei Beschran- 
kung auf kleine Auslenkungen erhalt man fiir die harmonischen Schwingungen die bekannte 
Formel fiir die Meridianschwingungsdauer T'1 des Kreiselkompasses 


B 
T = = 
27 SUcosp ’ 


1 Vgl. z. B. R. Grammel, Der Kreisel, Bd. 2, 2. Aufl., S. 109. Berlin/Gottingen/Heidelberg 1950. 
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die wegen der aus (19) folgenden Beziehung 


ih SRUcosp 
§ 
in den fiir alle Breiten konstanten Wert 
Dee oat i 84 Minuten (20) 
§ 


ibergeht.+ 


Aus (19) ergibt sich, daB der Drehimpulsvektor nicht genau in die wahre Nordrichtung i zeigt. 
Diese Auslenkung von %, aus der Nordrichtung ist der von geographischer Breite ~, Fahrtgeschwin- 
digkeit V und dem Kurs x abhangige Fahrtfehler 6: 


= COS % 
R 
tgd = (21) 
U cos y + R sin % 
oder, wenn man von dem um den Fahrtfehler gefalschten KompaBkurs Z = x — 6 ausgeht, 
: V cosZ 
sin 0 Tecaaois (21a) 


Der Fahrtfehler wird hervorgerufen durch das in (19) enthaltene Drehimpulsglied ~ SS ote. Da 


seine Richtung B x f£ stets nach Backbord gerichtet ist, wird es von Glitscher® als ,,natiirlicher 
Drall nach Backbord“ bezeichnet. 


Um die Beschleunigungseinwirkung besser zu erkennen, mége der beschleunigungsproportionale 
Vektor d*f/dt? entsprechend seiner Einwirkung auf den Kreisel in orthogonale Komponenten auf- 
geteilt werden: 

d*f 


df d ; é d 
= eG tbixgtet Fete ig Ae (22) 


dt 
Die zweite Beschleunigungskomponente von (22) b f x df/dt, die die Richtung des Kreiseldreh- 
impulses hat und die besagt, daB sich df/dt um die Vertikale dreht, ergibt mit dem Schwerpunkt- 


Vektor 6 = S f ein Drehmoment Ja £x (0 Ex) senkrecht zum Drehimpuls und lat nach 


dem PrazessionsgeSetz den Kreiseldrall sich um die Vertikale £ drehen, entsprechend dem Fahrt- 
fehler, der sich wegen der Drehung von df/dt um die Vertikale ebenfalls andert. Die erste Be- 
schleunigungskomponente von (22) a df/dt, die die Richtung von df/dt hat und die besagt, das der 
Vektor df/dt seinen Betrag andert, erzeugt ein Drehmoment parallel zum Kreiselimpuls; er andert 
seine Richtung nicht, jedoch seinen absoluten Betrag gemaB (18). Das durch Veranderung des 
Drehimpulsbetrages am Kreiselgehause wirksam werdende Reaktionsmoment ist entgegengesetzt 
gleich dem Beschleunigungsdrehmoment S f x peer Vektor © dreht sich nicht um die freie 
Kreiseldrehimpulsachse, sondern bleibt weiterhin vertikal gerichtet. Auf der Beachtung dieses 
Abgleichs der um die Kreiselachse wirksamen Drehmomente beruht die Raumkompah-Eigenschaft. 


Die vertikale Beschleunigungskomponente c f tragt wegen S f x —c f = 0 zur Momentenbildung 
§ 


iiberhaupt nichts bei. 


Man kann zusammenfassen: Ein Kompafkreiselsystem, das die Forderung (18) erfiillt und sich 
in der durch (19) und (19a) gegebenen Gleichgewichtslage befindet, bleibt, total unabhangig von 
Beschleunigungen, in dieser. Seine azimutale Gleichgewichtslage ist durch den Fahrtfehler gegeben, 
seine horizontale Gleichgewichtslage ist der wahre Horizont, mit anderen Worten: ein solches 
Kreiselsystem ist ein ,,Raumkompab“. 


5. Die verschiedenen prinzipiellen Méglichkeiten zur Realisierung des Raumkompasses. Die 
Erfiillung von Gleichung (18) gelingt am einfachsten, wenn man einen einzigen Kreisel benutzt 


1 M. Schuler, Physik. Z. 24 (1923) S. 344. 
2 K. Glitscher, Wissensch. Veroffentlichungen aus den Siemens-Werken 19 (1940) S. 191. 
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und seine Drehzahl @ in Abhingigkeit von y, V und x variiert gemaB der Gleichung 


Ae 2V Ucosy 
R R 


(sin x — 1) (23) 


SR : 
B=00=", |/(Ueosy | 


wobei 9 das Trigheitsmoment des Kreisels ist. 

Eine weitere Méglichkeit besteht darin, daB man den Drehimpuls des Kreiselsystems aufteilt 
in zwei Kreisel mit gleichbleibender Kreiseldrehzahl, deren Gehaduse auf motorischem Wege in 
einem Rahmen um parallele Achsen, z. B. die Hochachse, um entgegengesetzt gleiche Winkel 
gedreht werden kénnen. Durch die GréBe des Spreizwinkels 2 ist der Gesamtdrehimpuls gegeben. 
Die gemaf (18) bei horizontalen Beschleunigungen senkrecht zur Gesamtdrehimpulsachse erforder- 
liche Drehimpulsinderung wird durch von auBen gesteuerte motorische Veranderung des Spreiz- 
winkels 2 9} erzwungen. Das bei dieser zwangsweisen Gesamtdrehimpulsanderung sich entwickelnde 
Kreiselwirkungsmoment — parallel zum Gesamtdrehimpuls — wird wieder aufgenommen durch 
das Moment der am Kreiselsystem-Schwerpunkt angreifenden Beschleunigungskraft. 

Beim Anschiitz-Raumkompa®! wird Gleichung (18) automatisch durch einen Zweikreisel- 
verband? mit passender Abstimmung der Fesselungsmomente erreicht. Im Folgenden soll unter- 
sucht werden, welcher Winkelfunktion sie geniigen miissen. 


Fiir die beiden Kreisel mit den Drehimpulsen §, und %, gilt 


d&. yr 
Mn, = ze od es (24) 


Mt, und Nt, sind die auf Kreisel 1 und 2 insgesamt einwirkenden Drehmomente, also einschlieBlich 
der zwischen beiden Kreiseln wirksamen Reaktionsdrehmomente. Nach dem Prazessionsgesetz 
prazedieren die Kreisel im Raum mit den Winkelgeschwindigkeiten X$, und YW, : 


Seas eons We = 2S. x Be. (25) 


YW, und YW, kann man aufteilen in den Drehgeschwindigkeitsanteil Xp des Kreiselverbandrahmens 
und in die Relativdrehgeschwindigkeiten der Kreisel gegenttber dem Rahmen: 


XS, = Wp + BWira> . Wo = Wr + Wea: (26) 
Wegen der Kopplung der Kreisel gilt 
Beret = — Wire - (27) 
Unter Beachtung von (25), (26), (27) erhalt man durch Summation und Differenzbildung 
Wy + My = Wr X (By, + Be) + Wiser X (Bi — Be) » 
MN, — My = Wr x (Bi — Be) + Wire X (Bi + By) - 


SchlieBen die beiden Kreiselachsen den Winkel 2? miteinander ein, und bezeichnen r, q. 5) die 
orthogonalen Einheitsvektoren der Kreiselresultantenachse, der Querachse und der Hochachse des 
Kreiselrahmens, so kann man schreiben 


Bra = 9. (29) 
Durch skalare Multiplikation von (28,2) mit  vereinfacht sich diese Gleichung unter Beachtung 


von (29) zu 
D- Oa — We) =F: (Wr X (B, — B,)). 


In der Hochachse ist zwischen den Kreiseln in der Differenz der Momente nur das Federmoment 
wirksam: 


(28) 


F(8) = h- (My — M,) = isc (Wr xo — 8) . (30) 
Fir die Drehgeschwindigkeit des Kreiselrahmens im Raum kann man schreiben 


wenn fv die Drehgeschwindigkeit des Kreiselsystemrahmens relativ zur Erde bezeichnet. Somit 
wird aus (30) mit (31) 


F(3) = })- ((w + W*) x (B, — B,)) - (32) 


ae 1 J. W. ee ate 6 (1935) S. 229. oe Geckelerschen Gleichungen fiihren immer noch zu kleinen 

Wwingungen bel bahrzeugbewegungen. Der Grund liegt darin, daB er bei der Ableit i 

das Glied 2 By 6 d (sin e)/dt unberiicksichtigt gelassen he aan ee eae 
2 Vel. z. B. P. Christoph, Dtsch. Hydrograph. Z. 6 (1953) S. 207, 
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Da in der Gleichgewichtslage sich der Rahmen nur um die Vertikale drehen und fh) die Vertikale 
anzeigen soll, fallt ty aus (32) heraus. Somit vereinfacht sich (32) zu 


£-(U* x (6, — %,)) = F(9). (33) 
Andererseits mu, wenn 
8=8,+8,—"2tx4 umd S=SH=St 
et ae! po Un =Sh= (34) 
sein soll, 6, — B, die Richtung von df/dt haben: 
a a dt 
3 — B,=M—. (35) 
Gleichung (33) kann somit weiter umgeformt werden zu 
: df 
F(9) = ¥- (cus x Mz) (36) 
oder wegen (12) 
df\2 
x M(t x 7) = Fo). (37) 
Es muB also die vektorielle Differenz der Drehimpulse die Form haben 
Seas repeal pes ACES (38) 


dt df \2 
( a a) 


Durch Betrachtung der Betrage der Gleichungen (34) und (38) 
%, + B, = 2 By, cos ?, 
%,— 8, =2 By sin? , 


dt Chi 
A Aer 


9 


it x 


wobei B, der Drehimpuls der einzelnen Kreisel ist, ergibt sich die bekannte Forderung fiir das 
Federgesetz der Kreiselfesselung beim Anschiitz-RaumkompabB: 


F(8) = 4 Basin 9 cos 3 = (39) 


Jeder Kreisel muB mit dem Moment = F(#) in die statische Gleichgewichtslage ? = 90° vom 


Kreiselrahmen zuriickgezogen werden. Ist diese Forderung erfiillt, dann hat der Kompaf die 


Eigenschaft des Raumkompasses. Die Drehzahl der beiden Kreisel bleibt konstant, der Gesamt- 
drehimpuls des Kreiselsystems stellt sich automatisch auf den geforderten Wert on £ x ein. 
Eine horizontale Beschleunigungskomponente senkrecht zum Gesamtdrehimpuls, welche ein Dreh- 
moment um die Gesamtdrehimpulsachse erzeugt, laBt die beiden Kreisel gegeneinander im Innern 
des Kreiselrahmens prazedieren, so daB der Betrag des Gesamtdrehimpulses sich gemaB (18) andert. 

Bei dem KompaB von Glitscher! befindet sich im Kreiselsystem auBer dem Hauptkreiselsystem 
noch ein Hilfskreiselsystem, das den Drall nach Backbord kiinstlich erzeugen soll. Der Drehimpuls- 
betrag des Hilfskreiselsystems wird proportional der Fahrzeuggeschwindigkeit und seine Richtung 
auf motorischem Wege stets nach Backbord orientiert gehalten. Ist ~ der wahre Kurs des Fahr- 
zeuges, so soll also der Gesamtdrehimpuls des Kreiselsystems durch 


) SV ; 
ya see | “> (tsinz —q c0s 2) (40) 


gegeben sein. Der Vektor t gibt die Achsenrichtung des Hauptkreiselsystems, q die Querachse an. 
Es 1a8t sich leicht zeigen, daB die allgemeine Kompabschwingungsgleichung 


et R at 
8 _§4utxB=Sx(t+2 4) | 
tenet ; 2 B 2 Bxt U-hs-4 
=6x|f+ 2urxuxyt+2exu—S—1 +A | ay 


1 K.Glitscher, DRP. Nr. 599596, 1932. Kreiselkompa® fiir Fahrzeuge. Er beachtet jedoch nicht die 
g-Abhangigkeit des Drehimpulses des Hauptkreiselsystems. 
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mit dem Ansatz (40) die stationare Gleichgewichtslage 
te=t, Ge=i, He=t (42) 
enthilt. Infolge des hier kiinstlich eingegebenen Dralls nach Backbord entfallt die Auslenkung 


des Haupt-Kreiselsystems aus der wahren Nordrichtung um den sonst vorhandenen Fahrtfehler. 
Auch hier ist der KompaB ein Raumkompab. 


6. Die Erfiillung der Bedingung bei den handelsiiblichen Kreiselkompassen. Die strenge Vor- 
schrift fiir den Drehimpuls des Raumkompab-Kreiselsystems 


SIR SIC dee ERIS ’ ay ld 
8] == |x G1= 5 uxt— =" (Ucosy + zsin x) + 


B 
R 


ist bei den handelsiiblichen Kreiselkompassen nur annahernd erfillt. Es wird &/R gegen U x f 
vernachlassigt, vielfach auch die durch 11 x f =i U cosy gegebene Abhangigkeit von der geo- 
graphischen Breite g. Nur bei gewissen Kriegsmarine-Kompassen ist eine stufenweise Beriick- 
sichtigung der Breitenabhangigkeit fiir ziemlich groBe Breitenbereiche méglich. Infolgedessen 
miissen bei den handelsiiblichen Kompassen sowohl in der Azimut- wie in der Horizontanzeige 
kleinere Fehler auftreten. 


7. Die technischen Aussichten zur Realisierung des Raumkompasses. In den dreifiger Jahren 
versuchte die Firma Anschiitz, ihren Raumkompa® in einem Modellgerat! zu verwirklichen. Die 
gewonnenen experimentellen Ergebnisse entsprachen nicht den Erwartungen, da die theoretischen 
Forderungen wohl nicht voll erfiillt werden konnten. Es ist zu vermuten, dai bei der Lagerung 
der beiden Kreisel um die Hochachse im Innern der Kreiselkugel nicht die Reibungs- und Stér- 
momentenfreiheit erreicht wurde, wie sie angesichts der winzigen Federfesselungsmomente (39) 
erforderlich ist. Es ist méglich, daB auch magnetisch bedingte Stérmomente mitgespielt haben. 

Es erscheint aussichtsreich, die RaumkompaBeigenschaften tiber die Kreiseldrehzahlsteuerung 
zu gewinnen, allerdings nicht mit einem, sondern zwei drehzahlgesteuerten Kreiseln, die durch eine 
starkere Fesselung an eine praktisch konstante Gleichgewichtslage im Kreiselverband-Rahmen 
gebunden sind?. Bei Verwendung eines einzelnen Kreisels wiirde sich die Drehzahlsteuerung wegen 
der Tragheit des Umformers, falls man keinen Wechselrichter benutzt, nicht verzégerungsfrei 
durchfiihren lassen. Das wiirde insbesondere zu einem Schlingerfehler fiihren. Bei einem Kreisel- 
verband wiirde sich aber eine geringe Verzégerung infolge der Integrationseigenschaft der Kreisel 
kaum auswirken. Prinzipiell kénnte jedes handelsiibliche KreiselkompaBsystem, das punktweise 
gelagert ist (Anschiitz, Arma, Plath), hierfiir verwendet werden. Die Schwierigkeit des Anschiitz- 
Raumkompasses, geniigend reibungs- und stormomentenfreie Lager im Kreiselsystem-Innern her- 
zustellen, wiirde entfallen. Dafiir ware die Kreiselantriebsfrequenz und -spannung in Abhangigkeit 
von gy, V und x zu steuern. Es diirfte fiir viele Zwecke geniigen, die Fahrt mit einem iiblichen Fahrt- 
messer zu messen. Den Kurs liefert der KompaB selbst. Die nord-siidliche Fahrtkomponente 
B-i = Vcosx gibt die Breitenanderung 


dp _Vcosx 
di ai ae 


und das Integral die ganze Breitendnderung Aq, die sich gegeniiber der Breite des Ausgangspunktes 
der Reise eingestellt hat: 


Ap => | V cos x dt. 


Ein Rechengerat, dem V und x laufend eingegeben werden miBten, hatte den Sollwert fiir die 
Antriebsfrequenz zu liefern. 


8. Der RaumkompaB als Vorstufe fiir ein Inertial-Navigationssystem. Mit dem Raumkompa8 
kann aber auch ein ganzes Inertial-Navigations-MeBsystem aufgebaut werden, das auBer der Nord- 
richtung und dem wahren Horizont auch die Fahrt iiber Grund? und den Standort anzeigt. Hierfiir 
wire erforderlich, im Kreiselsystem in den Richtungen der beiden horizontalen Hauptachsen, der 
Kreiseldrehimpulsachse und der Querachse, je ein BeschleunigungsmeBorgan einzubauen, ihre MeB- 


: Vel. R. Grammel, Der Kreisel, 2. Aufl., Bd. 2, S. 140. Berlin/Gottingen/Heidelberg 1950. 
* Beim Arma-KompaB werden die beiden Kreisel bereits mit verschiedener Drehzahl stufenweise in Ab- 


hangigkeit von der Breite betrieben, vgl. A. L. Rawlings, The Theory of ; i i 
Deviations, S. 73, New York 1944, ai ebaan ccna aan ik Cs Ae 


° Naheres bei P. Christoph, Dtsch. Hydrograph. Z. 11 (1958) im Druck. 
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werte nach aufen zu iibertragen und da dem Rechenageregat zur Verarbeitung zuzufiihren, das 
dann bei gegebener Breite nach Integration und nach Abzug des durch die Erddrehung gegebenen 
Anteils die Fahrt iiber Grund, aber auch die Breitenanderung und die Langenanderung, vom Aus- 
gangspunkt gemessen, liefert. Die Steuerung der Kreiseldrehzahl miiBte allein von dem MeBwert 
des BeschleunigungsmeBorganes besorgt werden, das die Beschleunigungskomponente in Richtung 
der Querachse miBt, da die Drehbeschleunigung der Kreisel proportional zu dieser Beschleunigungs- 
komponente zu halten ist. 


9. Zusammenfassung. Bekanntlich muf bei allen Kreiselkompabtypen zur Vermeidung von 
Beschleunigungseinfliissen auf die KompaBanzeige bei Fahrt- und Kursanderungen die Meridian- 
schwingungsdauer des Kreiselsystems auf 84 Minuten abgestimmt werden. Die diesbeziigliche 
Bedingungsgleichung wird in Strenge fiir den schnellen Kreisel abgeleitet. Der Kreiseldrehimpuls 
ist danach durch das Schwererichtmoment vorgeschrieben, aber auch noch durch die geographische 
Breite, die Fahrtgeschwindigkeit und den Kurs. Wenn es gelingt, die Bedingungsgleichung exakt 
zu erfillen, gelangt man zum sogenannten RaumkompaB, welcher nicht nur unter Beachtung des 
Fahrtfehlers die Nordrichtung, sondern auch den wahren Horizont bei beliebigen Fahrtbewegungen 
des Fahrzeuges auf der Erdoberflache exakt anzeigt. Die Diskussion der Bedingungsgleichung 
fihrt iiber den bekannten RaumkompaBb von Anschiitz hinaus zu weiteren Verfahren, den Raum- 
kompaB zu verwirklichen. Es wird darauf hingewiesen, da bei den handelsiiblichen Kreiselkom- 
passen die Fahrtgeschwindigkeitsabhangigkeit in der Bedingungsgleichung stets vernachlassigt 
wird und die Abhangigkeit von der geographischen Breite meistens nur fiir eine Breite erfiillt ist, 
und daB infolgedessen mit kleineren Fehlern fiir Azimut und Horizontlage bei den Fahrzeug- 
manévern zu rechnen ist. SchlieBlich wird darauf hingewiesen, da mit der Schaffung des Raum- 
kompasses auch die Aufgabe der Inertial-Navigation, die Messung der Fahrt iiber Grund, im Prinzip 
zu lésen ist. 


(Eingegangen am 14. Oktober 1957.) 
Anschrift des Verfassers: Dr. Peter Christoph, GroBhansdorf bei Hamburg, Barkholt 3. 
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Ein Verfahren zur Druckverteilungsrechnung 
an geraden und radialen Schaufelgittern * 


Von H. Kriger 


1. Einleitung. Zur theoretischen Behandlung wird die Strémung durch Leit- und Laufrader 
von Turbomaschinen in bekannter Weise durch Abwicklung und Verebnung der FluBflachen in die 
Umstrémung eines geraden (bei Axialmaschinen) oder eines kreisférmigen (bei radial oder diagonal 
durchstrémten Maschinen) ebenen Schaufelgitters iiberfiihrt. Die potentialtheoretische Druck- 
verteilung dieser ebenen Strémung gibt, solange Grenzschichtablésung, Kompressibilitats- und 
Kavitationseffekte unberiicksichtigt bleiben kénnen, einen brauchbaren Anhalt fiir die tatsach- 
lichen Druckverhaltnisse und Auftriebskrafte am Fliigelschnitt; sie ist Ausgangspunkt fiir die 
Grenzschichtrechnung und damit fiir die Abschatzung des Widerstands- und Ablésungsverhaltens. 

Zur Berechnung der potentialtheoretischen Geschwindigkeit am Einzelprofil — der Druck- 
verlauf ist mit ihr durch die Bernoulligleichung verkniipft — existieren einfache, hinreichend 
genaue Verfahren. Liegt das Profil im Verband eines Schaufelgitters, so werden die Zusammen- 
hange kompliziert und genaue Rechnungen langwierig. Es ist hier von Vorteil, wenn man sich auf 
einmal durchgefiihrte Berechnungen stiitzen kann. Diese Arbeit enthalt ein Verfahren zur Ge- 
schwindigkeitsbestimmung an diinnen, gewélbten Profilen im geraden und radialen Gitter mit 
Hilfe vorausberechneter, tabellierter Geschwindigkeitsverteilungen. AuBer dem Parallelanstrom 
des geraden und dem Wirbelquellanstrom des Kreisgitters kann auch dessen Rotation erfaBt werden. 

Das ebene Schaufelgitter ist durch Belegung seiner Skelettlinien mit Wirbeln, Quellen und 
Senken darstellbar, deren Stérungsfeld den Anstrom zum UmflieBen der Profilkontur zwingt. Die 
Stérgeschwindigkeiten normal zur Kontur sind durch Gittergeometrie und Anstrémung bestimmt; 
die Tangentialkomponenten ergeben zusammen mit denen der Anstrémung die Geschwindigkeits- 
verteilung am Profil. 

Fir schwachgewélbte Schaufeln im geraden Gitter ist der Zusammenhang zwischen Singulari- 
tatenbelegung und Stérgeschwindigkeiten praktisch der gleiche wie beim Gitter mit geraden Skelett- 
linien. Beschrankt man sich weiterhin auf Wirbelbelegungen, so werden die Schaufeln unendlich 
diinn, und Normal- und Tangentialgeschwindigkeiten sind miteinander verkniipft wie bei Gittern 
aus ebenen Platten. 

Das vorliegende Verfahren behandelt das ,,direkte Gitterproblem“ (Gitterform und Anstré- 
mung vorgegeben, Geschwindigkeiten am Profil gesucht) und geht von den Normalgeschwindig- 
keiten aus. Fiir einfache, den ersten drei Gliedern einer Potenzreihe entsprechende Normal- 
geschwindigkeitsverlaufe sind die zugehérigen Tangentialgeschwindigkeiten an verschieden ge- 
staffelten und geteilten Plattengittern tabelliert. Sie entsprechen den bekannten Birnbaumschen 
Normalverteilungen der Einzelplatte! und ergeben durch Uberlagerung Lésungen des direkten 
Problems fiir diinne, maBig gewolbte Profile mit S-Schlag im geraden Gitter sowie — nach einer 
Transformation — fiir diinne, wenig von der logarithmischen Spirale abweichende Schaufeln im 
Kreisgitter. 

Damit ist das Verfahren Gegenstiick und Erginzung zu den Gitterberechnungsmethoden von 
Schlichting? und Scholz’. Schlichting lést das direkte Problem fiir makig gewdlbte und maBig dicke 
Profile im geraden Gitter, indem er Wirbel- und Quellbelegung durch je drei Glieder einer Birn- 
baum-Glauertschen Reihe darstellt, deren Stérgeschwindigkeiten fiir Plattengitter tabelliert sind. 
Die entstehende Kontur soll an drei Profilabszissen mit der vorgegebenen iibereinstimmen. Man 
erhalt ein lineares Gleichungssystem, aus dem sich die Anteile der einzelnen Belegungen und damit 
die Faktoren zum Aufbau der Geschwindigkeitsverteilung berechnen lassen. 


__™ Gekiirzte Fassung der Dissertation an der T. H. Braunschweig; Berichter Prof. Dr. H. Schlichting ; 
Mitberichter: Prof. Dr. A. Betz und Dozent Dr. N. Scholz. Die Arbeit wurde im Max-Planck-Institut fiir 
Strémungsforschung, Gottingen, durchgefiihrt. 

+ W. Birnbaum, Die tragende Wirbelflache als Hilfsmittel zur Behandlung des ebenen Problems der 
Pe ae Z. angew. Math. Mech. 3 (1923) S. 290. 

* HT. Schlichting, Berechnung der reibungsfreien inkompressiblen Stré e fur ei g 
Schaufelgitter. VDI-Forschungsheft 447 (1985). Fe ea ae 


foes, Scholz, Kin Berechnungsverfahren zum Entwurf yon Schaufelgitterprofilen. VDI-Forschungsheft 442 
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Ausgangspunkt zur Lésung des ,,Umkehrproblems“, Profilgitter von maBiger Wélbung und 
Dicke mit bestimmten aerodynamischen Eigenschaften zu erzeugen, sind auch bei Scholz nach den 
ersten Gliedern der Birnbaum-Glauertschen Reihe gestaltete Wirbel- und Quellbelegungen gerader 
Skelettlinien. Im Zusammenwirken mit wenig von der Sehnenrichtung abweichenden Parallel- 
str6mungen ergeben sie profilierte Schaufeln, deren Form aus Stérgeschwindigkeitstabellen leicht 
zu berechnen ist. 

Die Brauchbarkeit des vorliegenden Verfahrens fiir potentialtheoretische oder ihnen nahe- 
kommende wirkliche Strémungsverhiltnisse wird bestatigt durch Nachrechnung eines ebenen 
geraden Gitters mit bekannten theoretischen 
Geschwindigkeitsverteilungen und durch Ex- 
perimente an einem die ebene Strémung ver- 
wirklichenden Radialgitter. 


2. Gerade Gitter. a) Zusammenhang 
zwischen Stérgeschwindigkeiten, Pro- 
filform und Anstrémung. Bringt man 
ein ebenes, gerades Gitter aus diinnen, zu- 
nachst durchlassig gedachten Schaufeln in 
eine Parallelstrémung der Geschwindigkeit 
W.. mit den Komponenten W,,,, in und W,, 
senkrecht zur Gitterrichtung, so wird die 
Schaufelkontur mit normalen und tangen- 
tialen Geschwindigkeitskomponenten wy und 
wy durchsetzt (Abb. la). Um die Schaufeln 
undurchlassig, also zu Stromlinien zu machen, 
hat man die wy durch entgegengesetzt gleiche 
Normalgeschwindigkeiten w, zu kompensie- 
ren, die man sich von einer Wirbelbelegung! 
ys», wo s die Lauflinge der Schaufel ist, 
auf der Kontur induziert denken kann. Die Abb. lau. b. Anstrom-, Stérgeschwindigkeiten und Geschwindigkeits- 

dreiecke am Gitter aus gewélbten Profilen. 
beim hier behandelten ,,direkten Gitter- 
problem“ gesuchte Geschwindigkeitsverteilung Wy setzt sich aus den Tangentialkomponenten des 
ungestérten Anstroms und der von den Wirbeln herriihrenden Stérstrémung zusammen: 


Wr =u, + u,. (1) 
Wir machen die Schaufelkoordinaten x in und y, senkrecht zur Sehnenrichtung mit der Sehnen- 
lange / dimensionslos: 


2% 
ae 
Ursprung ist die Sehnenmitte; Vorder- und Hinterkante liegen dann bei § = 7 1, 7, =0. Fir 
maBig gewélbte Schaufeln 7,¢&>, deren Sehne mit der Gitterfront den Winkel fg einschlieBt, ist 


Uny=W, (cos Bs — a sins] Wu (sin bs + a cos Bs}; 


wr = W, (sin Bs + a cos Bs) —Wou (cos bs — al sin). 
Wenn wir voraussetzen, daB das Profil durch den Potenzansatz dritten Grades 
Noa ae Bee At Bs (3) 
gut darstellbar ist (A ist fiir die Wélbung, B fiir den S-Schlag maBgebend)*, folgt ein Normal- 
geschwindigkeitsverlauf vom zweiten Grade 


WwW, == —_ Wn = We (Goa A Gok = Cone) + We u (Gy, ae Ci4.8 =— Ge, é*) 2 (4) 


1 ¢.+-) bedeutet: Funktion von .. 
2 4 und B lassen sich entweder aus den Endneigungen gegen die Sehne oder aus zwei Ordinaten, z. B. 
bei 25 und 75% der Schaufeltiefe bestimmen (Abb, 2): 


2 
SS i= 1, 


(2) 


_ (dns ae dns ea dis dns 3 
- A= ("s)o,5 ar (Ns)_ 0,5? a B= ("s)o,5 ik (7s)_ O55)* (3 b) 


Die Zeiger V, H und spater M gelten fiir Vorder-, Hinterkante und Schaufelmitte. 
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dessen Faktoren G nur von Profilform und Staffelung des Gitters abhangen: 


Go. = B sin Bs — cos Bs, Gy, = — B cos f, — sin Bg, | 
G,, =—2Asin fg, G,, =2AcosB,, (5) 
G,, =—3 Bsin fg, G,, = 3 B cos pe: 


Als Zusammenhang zwischen den so gegebenen Normalgeschwindigkeiten w, und den gesuchten 
Tangentialgeschwindigkeiten w, der Stérstromung wird der Zusammenhang zwischen den Auf- 
wartsgeschwindigkeiten v und den Langsgeschwindigkeiten u des von den Profilsehnen gebildeten 
Plattengitters benutzt. 

b) Geschwindigkeits- und Druckverteilung, Schaufelkraft und Umlenkwirkung. 
Fiir konstante, linear und parabolisch verlaufende, mit dem Hinterkantenwert normierte Auf- 
wartsgeschwindigkeiten 

Vv, v Vv 

po ee (6) 

“0H “1 V2 ‘ 
sind die Langsgeschwindigkeitsverlaufe fiir Plattengitter verschiedener Teilungsverhaltnisse t/1 
und Staffelwinkel £, in Ziff. 7 (Tab. 1 bis 5) mitgeteilt; sie seien nach den Exponenten als nullte, 
erste und zweite Grundverteilung 


7 Ug = Cb The)» 


Abb. 2. Zur Bestimmung der Profilkonstanten. 


bezeichnet. Mit ihnen lautet die Geschwindigkeitsverteilung an der Kontur der diinnen Gitter- 
schaufel: 


W, = »(Go. ux + G,, uf + G,, us + sin Bs + ee cos) 
S 
dis. 
+ Wou(Gou us + G,,, ut + Cay.uf — 008 Bs + Ge sin Ba) (7) 
Die Druckverteilung folgt aus ihr nach der Bernoulligleichung 
p ++ Wk =konst., (8) 


wo @ die Dichte des strémenden Mediums ist. Die Schaufelzirkulation setzt sich aus den ebenfalls 
tabellierten Zirkulationen 


iy 1p. Le 
eines Eres oF et 
der einzelnen Grundverteilungen zusammen: 
ISUW, Goal’ + Cal On, 1S) UU Wew (Gon lp Cre ly Goes (9) 
Damit betragt nach dem Kutta-Joukowskyschen Satz die senkrecht zu W., gerichtete Schaufelkraft 
Bp Wale (10) 


Die Geschwindigkeiten W, weit vor und W, weit hinter dem Gitter sind von der Geschwindig- 
keit W,, durch den Zusatz =F ['/2t in Gitterrichtung unterschieden (Abb. 1b). 


c) StoB®freier Anstrom. In wirklicher Strémung zeichnet sich der Fall stoBfreien Eintritts 
durch geringe Verluste aus. Er liegt beim unendlich diinnen Profil dann vor, wenn die Geschwin- 
digkeiten an der Vorderkante endlich bleiben, d. h. wenn sich die unendlich groBen Geschwindig- 
keiten der nullten Grundverteilung und die der zweiten Grundverteilung gegenseitig aufheben. . 
Vom Verfassert ist gezeigt worden, daB sich die zweite Grundverteilung aus einem Anteil ohne 
Gesamtzirkulation mit endlicher Vorderkantengeschwindigkeit? und einem Anteil (/°*/I"*) u* der 
nullten Grundverteilung zusammensetzt, der fiir die Zirkulation J’, aufkommt?. Damit folgt aie (9) 


' H. Kriger, Ein Verfahren zur Druckverteilun h r r 
5 gsrechnung an geraden und radialen Schaufelgittern. 
Bericht 56/B/10 des Max-Planck-Instituts fiir Strémungsforschung, Gottingen. ‘ 
Es Bei Birnbaum ist dieser Anteil als Normalverteilung angegeben. 
Der wichtige Quotient [¥/I’, ist bei Klingemann (FuBnote 1 von Seite 254 mit S bezeichnet). 
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nach Weglassen der an der Vorderkante endlichen Glieder die Bedingung fiir stoBfreien Anstrom: 


ore (Go. EF af Coq fT) = We u stfr. (Cou i 8 Coy EP) = 0 . (11) 
d) Bezugnahme auf die Zustrémung. In praktischen Aufgaben ist gewohnlich die Zu- 
strémung W, mit den Komponenten W, = W, sin £,, W,,, = — W, cos fb, vorgegeben. Aus den 


Geschwindigkeitsdreiecken ergibt sich daraus fiir W,, die Komponente in Gitterrichtung 


io Wie 2 fy ; (12) 


(Die Komponente senkrecht zum Gitter bleibt aus Kontinuitatsgriinden ungeandert.) 


e) StaffelungseinfluB bei gewélbten Schaufeln. Beim Plattengitter mit 6,—=90° und 
fs = 180° betragen ebenso wie bei der Einzelplatte die von einer Wirbelbelegung y<x> induzierten 
Langsgeschwindigkeiten an Ober- und Unterseite u,, = + y/2, und die mittlere Langsgeschwindig- 


un 
keit auf der Kontur ist Null, was anschaulich aus der Symmetrie der Profilanordnungen und damit 
auch der Stérstrémungen folgt. Beim gestaffelten Plattengitter ist die mittlere Langsgeschwin- 
digkeit u = (u,,-+-u,,) von Null verschieden, wie man z.B. aus Abb. 6 ersieht. Sie wird, da u 
an der Einzelplatte verschwindet, von den Belegungen der anderen Schaufeln, des sogenannten 
Restgitters, induziert und Andert sich in nachster Umgebung der Platte nicht wesentlich. Ein 
an die Stelle der Platte gebrachtes diinnes, gewélbtes und durchlassiges Profil 7,¢&) wird von der 
Stérstrémung des aus Platten bestehenden Restgitters mit den Geschwindigkeiten normal und 
tangential zur Kontur durchsetzt: 
dns 


, 
wn = —u-= w 


dé ? = ie (13) 


WY 


Wahrend die Tangentialkomponenten gegeniiber dem Fall verschwindender Wélbung praktisch 
ungeandert bleiben, werden die Normalkomponenten wy des ungestérten Anstroms durch die 
Induktion des gestaffelten Restgitters in Verbindung mit der Profilwélbung um den (gewohnlich 
kleinen) Zusatz wy vermehrt, der durch eine aus unseren Grundverteilungen aufgebaute zusatz- 
liche Wirbelbelegung abgesattigt werden mu8. Dazu wird wy mit Hilfe der Konstanten A, B und 
der tabellierten u-Werte an Vorder-, Hinterkante und in Profilmitte durch einen Ausdruck vom 
zweiten Grade angendahert: 


wy = — wy = W, (Coo t+ Cia’ + Geek) + Wau (Cou + Gink + G2. 8). (14) 
Die Koeffizienten G’ sind durch die Gittergeometrie vollstindig bestimmt; sie lauten 
Goa =G) tim B+ C,,tfuy B+ 6,,ujm B, 
Be Cp, [ahr (4 > Bf aty (4B) — Cy, [eta (4+ B) Et (AB) 
— G24 [win (A+B) + uzy(A—B)], (15) 
Goa = Coa [wiv (A—B) —uj} ou B—upy(A + B)] + Gia [Ui y (A—B) —tlo B— un (4 + B)] | 
+ G24 [wy (A— B) —uim B—uyu(A + B)). 


Zur Bestimmung von Gy, Gi,, G2, ist der Zeiger a durch wu zu ersetzen. In den Gleichungen (7) 
fiir die Konturgeschwindigkeit, (9) fiir die Schaufelzirkulation und (11) fiir den stoBfreien Anstrom 
ist zur Beriicksichtigung des Staffelungseinflusses lediglich Gy + 6), G, + G1, Gy + Gi, anstelle 
von G,, G,, G, zu schreiben. 


3. Radiale Gitter. a) Konforme Abbildung eines radialen auf ein gerades Gitter, 


geometrische und Geschwindigkeitszusammenhange. Ein in einer t(=r e'°)-Ebene 
konzentrisch zum Ursprung liegendes n-schaufliges Kreisgitter mit den Innen- und AuBenradien 


rj, und dem mittleren Halbmesser ry = Vry rz wird durch die konforme Abbildung 


Pe ee (16) 
oD rM 


Ty 


in ein gerades Gitter der Teilung t mit dem Schaufelwinkel /, in einer z(=x-+iy)-Ebene iber- 
fiihrt (Abb. 3). Insbesondere gehen tiber: 

1. die Kreise vom Halbmesser ry, ry, rz in die Geraden in Gitterrichtung durch Vorderkante, Sehnen- 
mitte und Hinterkante, 
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2. logarithmisch-spiralige Schaufeln mit dem Neigungswinkel f, gegen Umfangsrichtung in Platten- 
profile der Neigung J, gegen die Gitterachse, 
3. der Radmittelpunkt 7 = 0 und das unendlich Ferne t =o in die unendlich fernen Punkte 
z = foo beiderseits des geraden Gitters. 

Eine von der logarithmisch-spiraligen abweichende diinne Schaufel sei gegeben durch das 
Halbmesserverhaltnis r,/ry, die Winkeldifferenz o—oy der Fahrstrahlen zur Aus- und Eintritts- 
kante und den Verlauf der Neigungswinkel /, der Kontur gegen die Umfangsrichtung mit den 


Endwerten fy, By. Die Vorder- und Hinterkante verbindende logarithmische Spirale wird zu 
einer Platte im geraden Gitter vom Tei- 


fs lungsverhAltnis 


nr Schaufeln ae Es el PATE, = Fa 
<= nym airy) + (or — 27) 
elogari 
[7] < ee : mit dem Schaufelwinkel 


=> pelt ee 
Bs = — are ctg ae (18) 
Das Schaufelbild wélbt sich tiber dieser 
Sehne mit den Anfangs- und Endtangenten 


dns dns 
se =bv—bs» Fe=Pba—Bs> (19) 


wird also durch Gleichung (3) angenahert 


mit 


<a 


DANES pa Nae 

4 B=2fs—fy —fy- (20) 
Ein Punkt der logarithmischen Spirale 
mit dem Abstand r vom Ursprung hangt 
mit dem entsprechenden Plattenpunkt zu- 
sammen durch 


rd 


ae BE (21) 


'M Ue 


Bei einer Kreisgitterschaufel, die nur 
die Enden mit der logarithmisch-spiraligen 
gemeinsam hat, ist der Zusammenhang 
komplizierter, weil nicht Punkte mit glei- 
chem Radius, sondern mit kiirzestem Ab- 


Abb.3a u. b, Transformation des Kreisgitters in ein gerades Gitter 
Lig ent ee, stand einander zugeordnet sein sollen. In 


durch ze 


ne der z-Ebene liegen zugeordnete Punkte 
von Sehne und Schaufelbild senkrecht zur Sehne iibereinander, wahrend die Linien konstanter 
Radien die Richtung der Gitterachse haben. Eine solche Linie durch den Punkt (x, y,) der Bild- 
schaufel schneidet die Sehne im Punkt x-+ y, ctg Bg (Abb. 3b). Fir stark von 7/2 abweichende 
Schaufelwinkel ist die Abszissenverschiebung y,ctg fg, auch bei geringen Abweichungen der 
Schaufel von der logarithmischen Spirale betrachtlich. An Stelle von (21) tritt der allgemeine 
Zusammenhang 


pe E+ et hg (44+ BE— AB BED] 
= =) (22) 
M V 

Fur die Ubertragung der Schaufelgeschwindigkeiten ist der Abbildungsmafstab an der Kontur 
maBgebend : 
dt 


pees == 9p 
dz\x 


nh Tp 


Ww, = W . 
z % nt 


(23) 


Tt 


b) Geschwindigkeitsverteilung beim Anstrom aus einer Wirbelquelle. Durch die 
Transformation (16) geht die ungestérte Strémung einer im Ursprung der t-Ebene gelegenen 
Wirbelquelle der Ergiebigkeit FE und der Wirbelstarke J", in eine Parallelstrémung in der z-Ebene 
mit den Komponenten W, = E/nt, W.,=/'./nt tber. Der Geschwindigkeitsverlauf (7) an 
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der Bildschaufel wird mit (23) auf die Kreisgitterschaufel iibertragen: 


Ua Deane L nat (Goa Ug ale Ci, ut +. Cy, Ug —- sin b,) 
te : ; 
ere rie (Go, us +, ux + 6,,, u* —cos f,). . (24) 


Die Schaufelzirkulation bleibt von der Transformation unberiihrt, betrigt also 


El 
ee Ga ire (Gels Gy. Tt 4G, 1s). 25) 

c) Rotierende Radialgitter. Beim mit konstanter Winkelgeschwindigkeit @ in unendlich 
ausgedehnter Fliissigkeit umlaufenden Radialgitter kreist die ungestérte relative Anstrémung 
(Schaufeln durchlassig gedacht) als starre Scheibe mit —@ um den Radmittelpunkt. Die Normal- 
komponenten kénnen wir, um die Schaufeln zu Stromlinien der Relativstrémung zu machen, 
mit Hilfe von Wirbelbelegungen kompensieren. Das von ihnen induzierte Stérfeld, dessen Ge- 
schwindigkeiten an der in Drehrichtung vorangehenden Schaufelseite von dieser weg, an der anderen 
auf sie zu gerichtet sind, ist identisch mit dem Absolutstrombild, einer Momentaufnahme des 
instationaéren Vorgangs. 

Ubertraigt man die Normalkomponenten wy = —rw sinf, der ungestérten Anstrémung auf 
das gerade Gitter, so erhalt man mit (22) und (23) die quadratische, an Schaufelmitte und -enden 
genau giiltige Naherung fiir die zu erzeugenden Aufwiartsgeschwindigkeiten (als charakteristische 
GréBe wird entsprechend E und /',, das mit der Umfangsgeschwindigkeit im Radius rj, gebildete 
Linienintegral benutzt) : 
ao 2a ri, o 


nm 


(Go w | CG, mes + Go w &) (26) 


nt 


mit den durch Gittergeometrie und Profilform bestimmten Konstanten 
1 
7 bv —Pa) ete Bb, 


r 
Ho 
Go op & 


. 1 > A) 
sin 7- (6 bs —PBy — Px) » | 
V 


rae Tot), 

2615 = ney — sin py » 
V H 
ripe Pier 

2G. = —4 sin By ae sin py —2 Go. - 
V H 


Die Geschwindigkeitsverteilung an der Schaufel, die im allgemeinen Fall zu (24) hinzuzufiigen 
ist, lautet 


is 2 
Wr y= yy 10) “u Go 7) ux = Gy w us Se Gy 7) us ai a cos b, > (28) 
M 
die Schaufelzirkulation betragt fiir den Rotationsanteil der Strémung 
1 2 l aa : 
Py =~ (Cool § + Galt + Gol). (29) 


d) Druckverteilung, Kraft- und Umlenkwirkung eines rotierenden, aus einer 
Wirbelquelle angestrémten Kreisgitters, An die Stelle der Bernoulligleichung tritt im 
rotierenden System die Energiegleichung der Relativstr6mung 


pt+ + [Wk — (r w)?] = konst. (30) 


Drehmoment M und theoretische Férderhéhe H,, folgen aus der Schaufelzirkulation und der 
Quellstarke bzw. der Drehgeschwindigkeit zu 


M=o 3-0, (31) 
Hy =5— Tl. (32) 


Der Drall wird durch das Gitter vom Eintrittswert /’, auf den Austrittswert /', geandert gemaB 


Sah es ow & 
y Bie yaaa +>: (33) 
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e) StoBfreier Anstrom. Analog (11) miissen bei stoBfreiem Hintritt Quellstarke, Wirbel- 
stirke und Winkelgeschwindigkeit der Bedingung gentigen 
Eosifr. (Goa 1% + G2, 1%) cee Ree (Goud + G2, IF) hs 2 7 Ti Wsifr. (Gog ls + Gs FF) 0G (34) 
f) Bezugnahme auf die Zustrémung. Aus den gewéhnlich vorgegebenen GréBen E, [4, w 
berechnet sich mit (33) die Drallkomponente der ungestérten Anstrémung zu 


2t 
EN(Gy ede oi Gioia Coto ae DL + 2272, 0 (6, LF Gy Lee Gye TF) 
Ss — —— ; (35) 
——— 


SE CA ee aati 


E und w bleiben ungeandert. 

g) Staffelungseinflu®. Bei gréBeren Abweichungen von der logarithmisch-spiraligen 
Schaufelform macht sich der unter 2e) beschriebene Einflu® der Staffelung geltend. Man hat 
dann die Konstanten G,, G,, G, in den Gleichungen fiir Konturgeschwindigkeit, Schaufelzirku- 
lation und stoBfreien Anstrom durch die korrigierten Werte G, + Gi, G,,+ Gi, G, + Gi, zu ersetzen. 


Gi, und Gi, folgen aus (15), wenn statt des FuBzeichens a u bzw. w geschrieben wird. 


4. Berechnungsbeispiele. a) Gerades Gitter, Vergleich mit vorliegenden theoretischen 
Ergebnissen. In Abb. 4 ist ein gerades Gitter t/] = 1,0, 6; = 135° aus diinnen Schaufeln der 


t/l=40 
1959 ——— unten 
a ABs ——tKeration | 2. 
ida o Bericht 


Abb, 4, Geschwindigkeiten am Parabelgitter, Vergleich mit Truckenbrodt*. 


a ee 
“a ee Oy & * By B,* Bs [Sx Ba Ban* 
ee ee 
1 0 0 —13,6:) 13,6 || 1214+) J21,4 121,4 | 121,4 | 121,4 | 121,4 
2 0,591 | 0,586 0,1 0 135,1 | 135,0 121,8 | 121,9 | 129,0 | 129,0 
3 1,0 1,0 8,2 8,2 | 14319) |" 143.9 122,429.49 1347 | siaa7 
4 1,151 Se 11,0 = 146,0 sia 122.8 = 136,8 = ws 
4a = 1,174 a ite ee 146,3 = 122,6 nt 1369 } OUES 
5 1,510 | 1,5 16,9 Moe |] TEL || isn PEM |] TIPE || aT ee || nan 


(Alle Winkelangaben in Grad) 


* Bei dem zum Vergleich herangezogenen Bericht handelt es sich um die in FuBnote 1 auf S, 249 angefithrte Arbeit. 
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Form 7 = 0,2(1—&?) mit seinen Geschwindigkeitsverteilungen skizziert. Kreise bezeichnen Re- 
sultate nach dem Schlichtingschen Verfahren, die von Truckenbrodt! mitgeteilt sind; die Kurven 
sind fiir jeweils gleiche W., und 8. nach unserem Verfahren berechnet. Der StaffelungseinfluB 
mu beriicksichtigt werden; die durch ihn bedingten Zusatzgeschwindigkeiten betragen Frets 
8% der Zustrémgeschwindigkeit. Die Ubereinstimmung ist gut. Kine Vergleichstabelle zeigt weit- 
gehendes Ubereinstimmen auch der tibrigen aerodynamischen Parameter. Dabei ist zu bedenken 
daBs heide Naherungsverfahren das Problem von verschiedenen Seiten angreifen: bei Schlichting? 
wird die Wirbelbelegung yy, hier der Aufwindverlauf v aus einfachen Funktionen aufgebaut. 


50 100 150mm, 


-O4 a 


Abb, 5a u. b. Radialgitter aus geraden diinnen Schaufeln; Abb. 6. Grundverteilungen, Mittelwerte und Zirkulationsbeiwerte 
Annaherung der Schaufelform nach Gleichung (3). fiir ¢/1 = 1,18, Bg = 120°. 


Die den StaffelungseinfluB léschende zusatzliche Wirbelbelegung gilt genau nur fiir gerade 
Skelettlinien. An einem gewélbten Skelett erzeugt sie auBer den gewiinschten neue, durch die 
Konturneigung gegen die Sehne bedingte Normalgeschwindigkeiten. Genau genommen sind daher 
iterativ weitere Zusatzbelegungen aufzubringen, doch erweist sich schon die erste als praktisch 
bedeutungslos. Fiir c, = 0 ist die hierdurch bewirkte Anderung der Geschwindigkeitsverteilung 
mit eingetragen; sie macht maximal 1% der Zustrémgeschwindigkeit aus. 


b) Radialgitter, durchgerechnetes Beispiel. Vorgegeben ist ein Radialgitter aus diinnen, 
geradlinigen Schaufeln nach Abb. 5: Schaufelzahl n = 12, AuBenhalbmesser ry = 409 mm, Innen- 
halbmesser r, = 279 mm, Winkeldifferenz zwischen Hinter- und Vorderkante og — oy = 0,222. 
Gesucht ist: 


1 FE. Truckenbrodt, Bericht S 52/18 des Lehrstuhls fir Stromungsmechanik der TH Braunschweig. 


2 Siehe FuRBnote 2 von Seite 242. 
17 
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Tabelle la. Rechnungsgang fiir ein Radialgitter 


Spalte 1 2 3 | 4 5 6 
ays Teilung Schaufelwinkel Profilkonstanten 
Halbmesser- ittel- d hériven 
ie em: werhaltnie® | halbniesse: parcicn Cie Wolbung S-Schlag 
berechnet aus rylty Vry "iy (17) (18) oe (20) 
Zeichen Palty "wr afl | Bs A B 
Zahlenwert 1,465 337 mm. TN 120,0° 0,0555 0,0038 
Spalte 9f | 8 | 9 | 10 11 | 12 13 | 14 | 15 
Benennung Koeffizienten der Grundverteilungen 
berechnet aus (5) (27) 
Zeichen G a | G, a CG, a G, u G, u G, u Gy o G, o G, wo 
Zahlenwert 0,497 | —0,096 | 0,010 | —0,868 | — 0,056 0,006 0,857 0,395 0,084 


: Mittelwerte der nullten Mittelwerte der ersten 
LORS aE: ETON OAR Grundverteilung Grundverteilung 
berecunet aus ib? fe 1,182; Ps = 120,0° aus den Tabellen interpoliert 
Zeichen TS We ea ip Wit was, uty it* yr | ee 
TZablonwertns. Auftragung in Abb. 6 £015 0,19.) 04g 0a noe 0,17 


=a 


Mittelwerte der zweiten 


Zirkulationsbeiwerte 
Benennung Grundverteilung der Grundverteilungen 
berechnet aus (ie aus den Tabellen interpoliert 
Zeichen usp | Us Us ry lie | rs 
Zahlenwert — 0,12 | 0,06 | 0,19 2,36 1,36 23 


1. Geschwindigkeitsverteilung an der Schaufel 


a) W, bei Anstrémung aus einer Quelle E, 
b) Wp, bei Anstrémung aus einem Wirbel J}, 
c) W,, bei Rotation mit der Winkelgeschwindigkeit , 


2. Summe der Schaufelzirkulationen 


a) n I’, bei Anstrémung aus einer Quelle E, 
b) nI'p, bei Anstrémung aus einem Wirbel J}, 
c) nI’, bei Rotation mit der Winkelgeschwindigkeit «. 


Tab. 1 gibt zusammen mit Abb. 6 (interpolierte Grundverteilungen) und Abb. 7 (Geschwindig- 
keiten an der Schaufel) den Rechnungsgang wieder. Um zu zeigen, inwieweit der geradlinige 
Schaufelverlauf durch Gleichung (3) angenahert wird, ist in Abb. 5b die danach berechnete Schaufel- 
form gezeichnet. 


5. Vergleich mit Versuchsergebnissen. In einem Versuchsstand, der die Verwirklichung der 
ebenen Strémung zulieB, wurden u. a. die Druckverteilungen an einem dem Berechnungsbeispiel 
entsprechenden Radialgitter (Abb. 8) beim Anstrom aus einer Wirbelquelle gemessen., Abweichend 
von 4b) waren die Schaufeln nicht Platten, sondern symmetrische verallgemeinerte Joukowsky- 
Profile? mit 10,4°% Dicke und 31% Dickenriicklage. Bei Berechnung der theoretischen Druck- 


1 Siehe FufBnote 1 von Seite 244. 
2 A. Betz, Kine Verallgemeinerung der Joukowskyschen Fligelabbildung. ZFW 15 (1924) S. 10. 
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Tabelle Lb. Rechnungsgang fiir ein Radialgitter (Fortsetzung) 


Spalte | 31 32 | 33 34 | 35 | 36 37 38 39 


Benennung Zusatzkoeffizienten der Grundverteilungen zur Beriicksichtigung des Staffelungseinflusses 
berechnet aus (15) ; Fas aye 
Zeichen Coa Cia Cie Cou Gin Ce AE: Cie Chon 
Zahlenwert — 0,000 | —0,006 | —0,013 0,001 0,011 0,027 |—0,001 —0,011 —0,035 


1,114 | 0,760 | 0,896 | 0,444 
1,160 | 0,876 | 0,904 | 0,426 
1,210 | 1,000 | 0,912 | 0,408 


’ 


’ 


Benennung korrigierte Koeffizienten der Grundverteilungen 
berechnet aus Spalte 7 plus Spalte 31 usf. . 7 
Zeichen Goa ete Coa Cia .c Gy a | G, a a Coa | Cou =P Cou Gin sr Giu | Gey al Gis Coot Goo Grt+G on GoatGoe 
Zahlenwert 0,497 — 0,102 | — 0,003 | — 0,867 | —0,045 0,033 0,856 0,384 | 0,049 
Spalte | 49 | 50 51 52 | 53 | 54 
Drallkomponente Ergebnis: 
der ungestorten Anstromung Summe der Schaufelzirkulationen 
Benennung : ; 
Quellanstrom |Drallanstrom Rotation Quellanstrom |Drallanstrom Rotation 
=o —=90) (E=o=—0)| (B=j{f,=0) ((£,-=o= 0)|\F=e0=0)) =F, = 0) 
berechnet aus (35) (33) 
Zeichen TjE Doi, leh t Noe 2h nT/E nTjP, nI/20r,o 
Zahlenwert 0,232 0,534 0,588 0,464 0,932 1,176 
Spalte | 55 | 56 57 | 58 | 59 60 | 61 | 62 
sae oes Ergebnis: 
Linea | Rotiger | Lauf- | _ ortlicher | | Geschwindigkeiten am Radialgitter 
Benennung g alb- | ange chaulelwinke ; 
geraden messer | Quellanstrom|Drallanstrom| Rotation 
Gitter der Radialgitterschaufel CO SO=0)/ GS @=0) | Gah, = 
berechnet aus | gewahlt (22) Gittergeometrie (24) | (28) 
Zeichen E r[ryq X/L sin B, | —cos B, | 2ar,y W,/E 2 yg W ,jT| Wolryy © 
| trragung 1 om 
Taio boet Bia 0,826 0,000 | 0,802 0,598 Auftragung in Abb 
— 0,8 0,855 0,080 0,816 0,978 
— 0,6 0,885 0,164 | 0,830 0,558 
— 0,4 0,919 0,255 | 0,844 0,538 | 
— 0,2 0,953 | 0,346 0,856 0,518 
— 0,0 0,990 | 0,444 | 0,867 0,499 
0,2 1,030 | 0,545 | 0,877 0,479 
0,4 1,070 | 0,647 0,887 0,461 
0,6 
0,8 
1,0 


» 


verteilung wurde der Dickeneinfluf in vereinfachender Weise nach der Kontinuitatsgleichung 
beriicksichtigt und der Geschwindigkeitsabfall an den Staupunkten durch den Riegelsfaktor+ 


dargestellt 


W profit Z tg (36) 


Ty ns >) 


a on t 
W piace y/1 ae (dD/dL)? oO 


wo D die Schaufeldicke, L die Lauflange ab Vorderkante, t, die Umfangsteilung und D, die Schaufel- 
dicke in Umfangsrichtung bedeutet. 


1 F. W. Riegels, Das Umstrémungsproblem bei inkompressiblen Potentialstromungen, I. und II. Mittei- 


lung. Ing.-Arch. 16 (1948) S. 373 und 17 (1949) S. 94 sowie 18 (1950) S. 359. er 
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Abb. 7. Geschwindigkeiten am Radialgitter von Abb. 5. 


Einlauf - 
“gleichrichter. 


Haltestabe 


Abb. 8. Ausgefiihrtes Radialgitter mit geraden profilierten Schaufeln 
(MaBe in mm). 


Aufgetragen wurde die Druck- 
differenz gegen den bei ungestérter 
Quellstrémung im Halbmesser ry 
herrschenden Wert, dimensionslos 
gemacht mit dem Staudruck 


x Q(E/2 0 ry)’ 


der ungestérten Quellstrémung im 
Halbmesser r;,. Abb. 9 zeigt den Fall 
stoBfreien Anstroms, einen Fall mit 
kleinerem und einen mit gréBerem 
Kintrittsdrall. Die gerechneten Druck- 
verteilungen sind fiir Platten ge- 
strichelt, fiir die profilierten Schau- 
feln ausgezogen. Punkte stellen MeB- 
werte dar. Der tatsachliche Druck- 
verlauf — auch seine Uberschneidung 
bei kleinem Eintrittsdrall — wird 
durch die Rechnung richtig wieder- 
gegeben. Fir die Abweichungen 
zwischen Messung und Rechnung 
ist die vereinfachte Abschatzung des 
relativ groBen Dickeneinflusses, eine 
Formungenauigkeit an der Schaufel- 
nase und die auch mit der Faden- 
sonde festgestellte Grenzschichtablé- 
sung im hinteren Profildrittel ver- 
antwortlich. 


6. Die Grundverteilungen. a) Zur 
Berechnung der Grundvertei- 
lungen. Die durchlassig gedachten 
Schaufeln eines ebenen, geraden Plat- 
tengitters kann man mit einer an der 
Hinterkante verschwindenden Wirbel- 
verteilung so belegen, daB die indu- 
zierten Aufwartsgeschwindigkeiten v 
mit der n-ten Potenz der von Plat- 
tenmitte aus gerechneten Lauflange 
gehen. Derzugehérige Langsgeschwin- 
digkeitsverlauf u,, wird als n-te Grund- 
verteilung bezeichnet; nullte, erste 
und zweite Grundverteilung sind fiir 
verschiedene Teilungen und Staffe- 
lungen in den Tabellen 1 bis 5 von 
Ziff. 7 mitgeteilt. 

Zu ihrer Berechnung wurde das 
Plattengitter mit der Weinigschen 


Abbildung? 
w= =o * PS Oy Stg G - 

gs Nr Ctg & e*) (37) 
auf den Kreis transformiert (Bezeich- 


nungen folgen in b)). Betrachtet man 
die komplexe Geschwindigkeit u—iv 


* F. Weinig, Die Strémung um die Schaufeln von Turbomaschinen, Leipzig 1935. 
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e 

R— ; n=% 4 
° Oberseite By=120° 8 5 nae Os fs= 120 
¢ Unierseite 7 ¢Unterseite 

F = 046 F = 66 

(stoBtreier Anstram) 
P-Py 
ao 


Te) 
| 


n= 
° Oberseite  ig=120° 
+ Unterseite 7 

F = 986 


Abb. 9a—e. Druckverteilungen an der Radialgitter- 
schaufel. 


gestrichelt; diinne, ausgezogen: profilierte Schaufel, 
Punkte: Messung. 
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als Potential, so hat es an entsprechenden Stellen von Kreis und Platte den gleichen Wert. Fiir 
nullte und erste Grundverteilung sind die Potentiale durch einfache Anordnungen symbolischer 
Singularitaten darstellbar! und fiihren auf ibersichtliche Zusammenhiange. 

Zur Ermittlung der zweiten Grundverteilung wurden die v mit Hilfe des Abbildungsmafstabs 


als Radialgeschwindigkeiten auf den Kreis transformiert; die zugehérigen Umfangsgeschwindig- 
keiten werden aus ihnen mit dem bekannten 


Poissonintegral? numerisch berechnet und dann 
zuricktbertragen. 

b) Geometrische Beziehungen zwischen 
Plattengitter und Bildkreis. Die konforme 
Abbildung (37) fihrt das Gitter der Teilung t, 
der Plattenlange | mit dem Schaufelwinkel /, in 
der z(=x-+iy)-Ebene (Platte in x-Richtung) 
iiber in den Kreis vom Halbmesser a in der 
Z(=aev” + ‘?)-Ebene, wobei die unendlich fernen 
Punkte beidseits des Gitters 2=+oco in 
Z=+ae% und die Plattenenden z = + 1/2 in 

—+ae?’H = Tae'%’vy tbergehen (Abb. 10). 
Ein Schaufelpunkt x ist durch den Bildkreis- 
winkel m bestimmt: 

cos @ sin p 


TUX ° 
oe —— To yh Bs Mr Tq Tj wo SKS) Bs arc tg Sit YW . 
(38) 


Mit « = 1/2 und y = Py gibt diese Gleichung das 
Teilungsverhaltnis an, wahrend fiir den Schaufel- 
winkel gilt 


ctg Bs = — tg py Sig y- (39) 
Differentiation von (37) fiihrt auf den Abbildungs- 
mafstab 
2aa\dz| al dé 
t dZ K Lie dp 
Abb. 10a u, b. SR arte etn eel Saga ee td = cos Bs Sin Wy COS D Be sin Bs Cof Wo sin y (40) 
D 


Coj2 py cos 2 p 
mit dem ein Plattenelement |dz| = dx in ein Konturelement |dZ| = a dp des Kreises iitberfiihrt wird. 
Beim Ubergang zum LEinzelprofil (t/l>0oo) wird aus (37) die Joukowskysche Abbildung 
2=Z+07/Z oe der Plattenlange /=4a, der Punktzuordnung 2x/l=cosq und dem Abbildungs- 
a|dz| d& ax 
I \dZ\, dp 
c) Die Grundverteilungen, ihre Grenzfalle und Zirkulationsbeiwerte. Zu kon- 
stanter Aufwartsgeschwindigkeit vy = vp q gehért die nullte Grundverteilung 


2 
mastab sing. 


sin 2 Oy — tg a 
ie = = 
v Sin 2 wy, 5 ca 
Zu linear verlaufender Aufwartsgeschwindigkeit v, = vn gehort die erste Grundverteilung 
Det cos i 
us = ——(cos fy Ur F ? i eset 
5 ak Bs Ue Tq Core + sin fg arc tg San (42) 


Zu parabolisch verlaufender Aufwartsgeschwindigkeit v, = v2 q& gehért die zweite Grundver- 
teilung uy, die mit dem Poissonintegral numerisch berechnet wurde. 
Im Grenzfall der Einzelplatte entstehen die Birnbaumschen Normalverteilungen 


_1/l1— 
uy = —tg t= F = 
uy = sin @ == Vea (43) 
, lip Tiers 
u3 =— 7 (sin29 + ts F) = +E rarea| "1 2s 


1 W. Klingemann, Berechnung der theoretischen K inl i 
Ree UE Bear ona enaLiane SOS ec von Turbomaschinen. Ing.-Arch. 11 (1940) 
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bei sehr dicht stehenden Gittern erhalt man 


uf =—etgy, uf =—EFetgf,, uf ——Setg fy. (14) 
Aus der Betrachtung des unendlich dichten Gitters als Wirbelband mit in Gitterrichtung kon- 
stanter Intensitat — hei ihm fallen die induzierten Geschwindigkeiten in Gitterrichtung — folgt 
auch anschaulich u=—vetg fs. Wie die in Abb. 6 aufgetragenen Grundverteilungen zeigen, wird 
der Charakter der bekannten Birnbaumverteilungen beibehalten. Durch den unter Ziff. 5 beschrie- 
benen RestgittereinfluB ist der Mittelwert w = 3(u,, + u,,,) im allgemeinen von Null verschie- 
den, insbesondere ist an der Hinterkante u > 0 fiir 8, > 2/2. Mit abnehmender Teilung riickt, 
wie ein Vergleich zeigt, der Wirbelschwerpunkt der nullten Grundverteilung zur Vorderkante. 
Die bekannten Zirkulationsbeiwerte fiir die Einzelplatte mit konstanter, linearer und para- 
bolischer Aufwartsgeschwindigkeit 
14 7 
Pe ors Lt ers (45) 
werden im Gitterfall abgedndert auf 
aya 
eects. ‘ Lana OS : 
i ee T Sindy, VGoj 2 Wy — cos 2 My, (46) 


24 e 
re = (7) —In Cig yo - | (47) 


Fir [’¥ ist von Klingemann! eine Reihenentwicklung angegeben worden; in vorliegender Arbeit 
ergaben sich die Zahlenwerte /’* aus der Hinterkantenbedingung im Verlauf der numerischen 
Berechnung der zweiten Grundverteilung. Die Gleichungen (46) und (47) lauten fiir sehr dicht 
stehende Schaufeln 
2t 
Lee, * 
Ly =ly I sin Bs ° 


Nach Klingemann gilt dieser Wert auch fiir /*. 


d) Tabellen. Die in Ziff. 7 mitgeteilten Tabellen 1 bis 5 enthalten fiir Plattengitter der Tei- 
lungen t/l = 0,5; 1,0; 2,0; co und Schaufelwinkel 6, = 90°; 120°; 135°; 150°; 180° — die wp und 
Py Sind unten links angegeben — in Spalte 1 den vom Hinterkantenpunkt ab gerechneten Bild- 
kreiswinkel ~—qqy, in Spalte 2 die zugehérige Lauflinge & aus (38), in Spalte 3 den Abbildungs- 
maSstab (40), der zwar fiir das vorliegende Berechnungsverfahren nicht benétigt wird, aber fiir 
die Behandlung anderer Gitterprobleme niitzlich ist; in den Spalten 4, 5 und 6 die nullte, erste 
und zweite Grundverteilung mit ihren Mittelwerten an-den Plattenenden und in Plattenmitte und 
den Zirkulationsbeiwerten nach 6c). Die Grundverteilungen sind iiberall dort anwendbar, wo zu 
analytischen Funktionen mit konstant, linear und quadratisch am Streckenprofilgitter verlaufen- 
dem Real- bzw. Imaginarteil der andere Teil gesucht ist, z. B. bei der Erzeugung dicker Profil- 
konturen aus dem Plattengitter durch Verschiebungsvektoren, deren Komponenten in und senk- 
recht zur Plattenrichtung Real- und Imaginarteil analytischer Funktionen der Gitterebene sind?. 
Hierzu geeignete Tabellen, die auBer den behandelten auch die nullte und zweite Grundverteilung 
bei verschwindender Gesamtzirkulation umfassen, sicehe FuBnote 1 von Seite 244. Es gelten fol- 
gende Symmetrien: 


(48) 


wi Ti Bs ) Sr vod — Bs > (49) 
Tir Bs) pees oa 73 wt —s ? (50) 


ECT Asi P— 9m) = EC 73 7—Bss 22 — (Pn) ). (51) 
dé | t 


73% —Pss 22 —(—$u) ) (52) 


PL 75 Bes 90H) = Ppa Bs 20—(P — 9m) S- (54) 


1 Siehe FuBnote 1 von Seite 254. dee ae 
2 H. Kriiger, Ein Typ von dicken Gitterprofilen mit leicht berechenbarem Geschwindigkeitsverlauf. 
Bericht 55/B/28 des Max-Planck-Instituts fiir Stromungsforschung, Gottingen. 
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7. Tabellen 


Tabelle la. Konforme Abbildung auf den Einheitskreis und Grundverteilungen fiir 


Plattengitter 
pecan ee etree Wert Se a ee 
P?—?H 2 ml d& ak u* u* 
[°] S cre dp 0 1 2 
ee 2 ee ae es 

0 1.0000 0,0000 0,0000 0,0000 0,0000 
5 0,9691 — 11,6150 —0,0038 — 0,2513 — 0,2793 
10 0,7399 — 9,2648 —0,0076 — 0,3530 — 0,4510 
15 0,6291 — 6,9784 —0,0113 —0,3974 — 0,4809 
20 0,5431 — 5,5169 — 0,0152 — 0,4212 — 0,4032 
30 0,4149 — 3,8980 — 0,0232 — 0,4456 — 0,3222 
40 0,3193 — 3,0672 —0,0315 — 0,4574 —0,2696 
50 0,2415 — 2,5875 — 0,0403 — 0,4642 — 0,2176 
60 0,1740 — 2,2950 —0,0499 —0,4685 —0,1761 
10 0,1129 — 2,1184 —0,0605 — 0,4708 —0,1360 
80 0,0556 — 2,0229 —0,0725 —0,4721 —0,0978 
90 0,0000 —1,0925 —0,0865 —0,4725 —0,0613 
100 —0,0556 20229 —0,1030 — 0,4721 —0,0266 
110 —0,1129 —2,1184 —0,1235 —0,4708 0,0055 
120 — 0,1740 — 2,2950 —0,1497 —0,4685 0,0346 
130 —0,2415 — 2.5875 —0,1854 — 0,4642 0,0575 
140 — 0,3193 — 3,0672 =i) o%5 —0,4574 0,0789 
150 — 0,4149 — 3,8980 — 0.3227 — 0,4456 0,0769 
160 — 0,5431 —5,5169 — 0,4903 Se 0,044.7 
165 — 00,6291 — 6,9784 — 0,6567 —0,3974 0,0072 
170 —0,7399 — 9,2648 — 0,9882 — 0,3530 — 0,2552 
175 —0,9691 —11,6150 —1,9801 —=0201'5 aa bl rir 

180 —1,0000 0,0000 fore) 0,0000 fore) 
185 — 0,9691 11,6150 1,9801 0,2513 I PATEH 
190 — 0,7399 9,2648 0,9882 0,3530 0,2552 
195 — 0,6291 6,9784 0,6567 0,3974 —0,0072 
200 — 0,5431 5,5169 0,4903 0,4212 — 0,0447 
210 —0,4149 3,8980 073229 0,4456 —0,0769 
220 —0,3193 3,0672 0,2375 0,4574 —0,0789 
230 — 00,2415 2,5875 0,1854 0,1642 —0,0575 
240 —0,1740 2,2950 0,1497 0,4685 — 0,0346 
250 —0,1129 2,1184 0,1235 0,4708 —0,0055 
260 —0,0556 2,0229 0,1030 0,4721 0,0266 
270 0,0000 1,0925 0,0865 0,4725 0,0613 
280 0,0556 2,0229 0,0725 0,4721 0,0978 
290 0,1129 2,1184 0,0605 0,4708 0,1360 
300 0,1740 2,2950 0,0499 90,4685 0,1761 
310 0,2415 2,5875 0,0403 0,4642 0,2176 
320 0,3193 3,0672 0,0315 0,4574 0,2696 
330 0,4149 3,8980 0,0232 0,4456 023222 
340 0,5431 59,9169 0,0152 0,4212 0,4032 
345 0,6291 6,9784 0,0113 0,3974 0,4809 
350 0,7399 9,2648 0,0076 0,3530 0,4510 
355 0,9691 11,6150 0,0038 0,2513 0,2793 
360 1,0000 0,0000 0,0000 0,0000 0,0000 

Tr ut ws 

Tr = 0,5 
V 0,0000 0,0000 0,0000 
Bs = 90° M 0,0000 0,0000 0,0000 
H 0,0000 0,0000 0,0000 
Wo = 0,0864 

Wor eS a as Le 

0,9983 0,7804 0,7080 
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Tabelle 1b. Konforme Abbildung auf den Einheitskreis und Grundverteilungen fiir 


Plattengitter 
ee 
PRO. ? al dé wk uk * 
[°] Ss gras. 0 1 Mg 
psc Ne ee LS ol A a a a ee a 
0 1,0000 0,0000 0,5744 0,3855 0,3856 
5 0,9120 — 17,9672 0,5710 0,1218 — 0,0454 
10 0,8118 — 6,3311 0,5670 —0,0019 — 0,1462 
5; 0,7338 —4,9935 0,5630 —0,0767 — 0,2082 
20 0,6710 — 4,1002 0,5598 —0,1294 — 0,2188 
30 0,5735 — 3,0405 0,5518 — 0,2032 —0,1996 
40 0,4978 — 2,4662 0,5438 — 0,2558 —0,1995 
50 0,4344 — 2,1255 0,5358 — 0,2975 —0,1805 
60 0,3785 —1,9176 0,5261 — 0,3329 —0,1593 
70 0,3271 —1,7954 0,5165 — 0,3644 —0,1379 
80 0,2782 — 1,7364 0,5045 — 0,3936 —0,1145 
90 0,2302 —1,7321 0,4917 — 0,4216 —0,0919 
100 0,1816 — 1,7820 0,4757 — 0,4493 — 0,0654 
110 0,1345 — 1,8881 0,4565 —0,4750 — 0,0367 
120 0,0757 — 2,0824 0,4308 —0,5067 — 0,0094 
130 0,0139 — 2,3964 0,3972 — 0,5386 0,0149 
140 — 0,0593 — 2,9215 0,3468 —0,5745 0,0560 
150 — 0,1524 — 3,8843 0,2659 — 0,6168 0,1009 
160 —0,2855 — 6,0255 0,1057 — 0,6689 0,1180 
165 — 0,3833 — 8,3194 — 0,0553 — 0,6990 0,1156 
170 —0,5277 — 13,1120 — 0,3788 —0,7265 0,0191 
175 — 00,7929 — 20,7755 —1,3238 — 0,6739 —0,5274 
180 —1,0000 0,0000 fore) — 0,3855 co 
185 —0,9120 71,9672 2,4722 | —0,1218 1,2310 
190 —0,8118 6,3311 155192 0,0019 0,4930 
195 — 0,7338 4,9935 1,2037 0,0767 0,2220 
200 —0,6710 4.1002 1,0427 0,1294 0,1054 
210 —0,5735 3,0405 0,8825 0,2032 0,0225 
220 — 0,4978 2,4662 0,8016 0,2558 — 0,0264 
230 — 0,4344 2AQ55 0,7520 0,2975 — 0,0353 
240 —0,3785 1,9176 0,7175 0,3329 — 0,0308 
250 —- 0,3271 1,7954 0,6919 0,3644 — 0,0202 
260 — 0,2782 1,7364 0,6727 0,3936 — 0,0016 
270 — 0,2302 BPA! 0,6567 0,4216 0,0189 
280 —0,1816 1,7820 0,6439 0,4493 0,0475 
290 —0,1345 1,8881 0,6319 0,4750 0,0811 
300 — 0,0757 2,0824 0,6222 0,5067 0,1192 
310 — 0,0139 2,3964. 0,6126 0,5386 0,1694 
320 0,0593 2.9215 0,6046 0,5745 0,2291 
330 0,1524 3,8843 0,5966 0,6168 0,3228 
340 0,2855 6,0255 0,5886 0,6689 0,4422 
345 0,3833 8,3194 0,5854 0,6990 0,5457 
350 0,5277 13,1120 0,5814 0,7265 0,6637 
355 0,7929 20,7755 0,5774 0,6739 0,7489 
360 | 1,0000 0,0000 0,5744 0,3855 0,3856 
Ths whe ux 
+ 05 
V — 0,2100 — 0,3855 —0,2010 
ole 02 M 0,4980 0,0000 0,1143 
A 0,5744 0,3855 0,3856 
Wo = 0,0622 
Pq = 9,0359 Vs Uy Ty 
1 V525 0,8841 0,7736 
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Tabelle le. Konforme Abbildung auf den Einheitskreis und Grundverteilungen fiir 


Plattengitter 
BE —E—E——————— 
Y—?H - al d& * * u* 
Cy é ees ati} it 2 
ow Sh ee ee 
0 1,0000 | 0,0000 0,9975 0,6467 0,5259 
5 0,904.4 — 71,4429 0,994.4 0,3266 0,1095 
10 0,8166 — 5,4409 0,9911 0,1871 0,0215 
15 0,7505 —4,1817 0,9880 0,0990 — 0,0863 
20 0,6984 — 3,3893 0,9847 0,0348 — 0,1097 
30 0,6181 — 2,4971 0,9780 —0,0581 —0,1309 
40 0.5560 — 2,0206 0,9710 — 0,1266 — 0,1333 
50 0,504.2 —1,7399 0,9634 — 0,1822 — 0,1314 
60 0,4584 —1,5678 0,9554 — 0,2301 — 0,1242 
70 04164 — 1,4676 0,9465 — 0,2735 —0,1110 
80 0,3764 — 1,4198 0,9363 —0,3141 — 0,0970 
90 0,3371 —1,4160 0,9246 —0,3535 — 0,0860 
100 0,2974 —1,4564 0,9106 — 0,3930 — 0,0679 
110 0,2557 —1,5427 0,8934 — 0,4336 — 0,0508 
120 0,2107 —1,7057 0,8712 —0,4771 — 0,0241 
130 0,1601 —1,9651 0,8411 —0,5249 — 0,0004 
140 0,0999 — 2,4026 0,7971 —0,5805 —0,0195 
150 0,0232 — 3,2126 0,7252 — 0,6489 0,0571 
160 — 0,0896 — 5,0060 0,5837 — 0,7358 0,0939 
165 Osis —7,1247 0,4432 — 0,8055 0,0980 
170 —0,2984 | —12,1712 0,1635 — 0,8887 0,0810 
175 —0,5641 | —31,6412 — 0,6874 —0,9907 — 0,1932 
180 — 1,0000 0,0000 co — 0,6467 co 
185 -— 0,9044 7,4429 2,6088 — 0,3266 1,1602 
190 — 0,8166 5,4409 1,8316 —0,1871 0,5059 
195 ==), 70 4,1817 1,5519 — 0,0990 0,2675 
200 — 0,6984 3,3893 1,4114 — 0,0348 0,1681 
210 —0,6181 2,4971 1,2699 0,0581 0,0523 
220 —0,5560 2,0206 1,1980 0,1266 0,0092 
230 — 0,5042 1,7399 1,1540 0,1822 —0,0118 
240 — 0,4584 1,5678 1,1239 0,2301 —0,0185 
250 — 0,4164 1,4676 1,1017 0,2735 ~@ | —0,0136 
260 — 0,3764 1,4198 1,0845 0,3141 — 0,0040 
270 —(eRi7il 1,4160 1,0705 0,3535 0,0056 
280 — 0,2974 1,4564. 1,0587 0,3930 0,0250 
290 — 0,2557 1,5427 1,0486 0,4336 0,0467 
300 — 0,2107 1,7057 1,0397 0,4771 0,0817 
310 —0,1601 1,9651 1,0316 0,5249 0,1192 
320 — 0,0999 2,4026 1,0241 0,5805 0,1620 
330 — 0,0232 3,2126 1,0171 0,6489 0,2403 
340 0,0896 5,0060 1,0104 0,7358 0,3617 
345 0,1703 71,1247 1,0071 0,8045 0,4518 
350 0,2984 12,1712 1,0040 0,8887 0,6084 
355 0,5641 31,6412 1,0007 0,9907 0,8660 
360 1,0000 0,0000 0,9975 0,6467 0,5259 
us ait ut 
Tie 0,5 
V — 0,6000 — 0,6467 — 0,4765 
ee BSS M 0,8590 0,0000 0,1769 
H 0,9975 0,6467 0,5259 
Yo = 0,0364 
Pq = 0,0364 ry rt Ts 
1,4096 1,0540 0,8845 
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Tabelle 1d. Konforme Abbildung auf den Einheitskreis und Grundverteilungen fiir 


Plattengitter 
a a 
| 
%— 9H L dé 
ry i ae = u* ue ux 
0 1,0000 0,0000 Vo732t 1,0659 0,7507 
2 0,9450 5,2482 1,7313 0,8954 0,3015 
5 0,8678 — 17,5594 1,7303 0,5715 0,1900 
10 0,7864 — 4,6491 1,7282 0,4003 0,0554 
15 0,7319 — 3,3456 1,7260 0,2945 — 0,0081 
20 0,6910 — 2,6306 1,7242 0,2177 — 0,0336 
30 0,6299 — 1,9026 1,7198 0,1057 — 0,0485 
40 0,5838 —1,4859 1,7154 0,0230 — 0,0866 
50 0,5459 —1,2636 1,7110 — 0,044.4 — 0,0899 
60 0,5129 —1,1295 1,7058 — 0,1027 — 0,0895 
10 0,4827 — 1,0496 1,7000 —0,1555 — 0,0868 
80 0,4543 — 1,0089 1,6939 —0,2051 — 0,0850 
90 0,4264. —1,0005 1,6864 — 0,2534 — 0,0776 
100 0,3990 —1,0231 1,6776 — 0,3020 —0,0705 
110 0,3693 —1,0804 1,6667 —0,3518 — 0,0562 
120 0,3380 —1,1825 1,6531 — 0,4054 — 0,0436 
130 0,3031 —1,3509 1,6342 —0,4647 | —0,0279 
140 0,2620 — 1,6329 1,6065 — 0,5338 —0,0085 
150 0,2104 — 2,1456 1,5617 — 0,6198 0,0263 
160 0,1377 — 3,2684 1,4730 — 0,7382 0,0724 
165 0,0849 — 4,5004 1,3852 —0,8219 0,0752 
170 0,0067 — 17,2482 1,2099 — 0,9408 0,0983 
175 —0,1481 | —18,6268 0,6860 —1,1494 0,1670 
178 — 0,4861 — 71,4878 — 0,8845 —1,1611 — 0,2012 
180 — 1,0000 0,0000 fore) — 1,0659 lore) 
182 — 0,9450 59,2482 4,3487 — 0,8954 1,6505 
185 —0,8678 17,5594. 2,7782 —0,5715 0,8589 
190 | —0,7864 4,6491 2,2543 — 0,4003 0,3266 
195 | —0,7319 3,3456 2,0791 —0,2945 0,1739 
200 —0,6910 26306 1,9912 | —0,2177 0,1040 
210 — 0,6299 1,9026 1,9025 —0,1057 0,0456 
220 — 0,5838 1,4859 1,8577 — 0,0230 — 0,0132 
230 — 0,5459 1,2636 1,8300 0,044.4 — 0,0285 
240 —0,5129 1,1295 1,8112 0,1027 — 0,0352 
250 — 0,4827 1,0496 Nes) 0,1555 — 0,0366 
260 —0,4543 | 1,0089 1,7866 0,2051 — 0,0372 
270 — 0,4264. 1,0005 1,7760 0,2534 — 0,0314 
280 — 0,3990 1,0231 1,7703 0,3020 — 0,0227 
290 — 0,3693 1,0804 1,7642 0,3518 — 0,0059 
300 — 0,3380 1,1825 1,7585 0,4054 0,0107 
310 — 0,3031 1,3509 1,7532 0,4647 0,0334 
320 — 0,2620 1,6329 1,7488 0,5338 0,0648 
330 — 0,2104 2,1456 1,7444 0,6198 0,1215 
340 —0,1377 3,2684 1,7400 0,7382 0,2100 
345 — 0,0849 45004 1,7383 0,8219 0,2972 
350 — 0,0067 7,2482 1,7361 0,9408 0,3695 
355 0,1481 18,6268 Ihnssole 1,1494 0,7071 
358 0,4861 71,4878 1,7329 1,1611 1,1478 
360 1,0000 0,0000 | Iisyall | 1,0659 0,7507 
ue u* ux 
7 = 0,5 
V — 0,9100 — 1,0659 —0,6110 
(oie = AUX M 1,4520 0,0000 0,2580 
A 1,7321 1,0659 0,7507 
Wo = 0,0119 
Pq = 90,0197 rs vy Ty 
1,9999 1,4224 1,0307 
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Tabelle 5. Konforme Abbildung auf den Einheitskreis und Grundverteilungen fiir 
die Einzelplatte 


a 


ee E Bas u* u* u* 
[°] ‘ dp 0 1 2 

0 1,0000 0,0000 0,0000 0,0000 0,0000 
10 0,9848- | —0,1737 | —0,0875 | —0,1737 | —0,2148 
20 0,9397 — 0,3420 — 0,1763 — 0,3420 — 0,4096 
30 0,8660 — 0,5000 — 0,2680 — 0,5000 — 0,5670 
40 _— 0,7660 — 0,6428 — 0,3640 — 0,6428 — 0,6744 
50 0,6428 — 0,7660 — 0,4663 — 0,7660 — 0,7256 
60 0,5000 — 0,8660 —0,5774 — 0,8660 —0,7217 
70 0,3420 — 0,9397 — 0,7002 — 0,9397 —0,6715 
80 0,1737 — 0,9848 — 0,8391 — 0,9848 —0,5905 
90 0,0000 — 1,0000 — 1,0000 —1,0000 — 0,5000 
100 —0,1737 — 0,9848 —1,1918 — 0,9848 — 0,4249 
110 — 0,3420 — 0,9397 — 1,4282 — 0,9397 — 0,3927 
120 — 0,5000 — 0,8660 —1,7321 — 0,8660 — 0,4331 
130 . — 0,6428 — 0,7660 — 2,1445 — 0,7660 —0,5799 
140 —0,7660 | —0,6428 — 2,7475 — 0,6428 — 0,8814 
150 — 0,8660 — 0,5000 — 3,7321 — 0,5000 — 1,4331 
160 — 0,9397 — 0,3420 — 95,6713 — 0,3420 -—2,9142 
170 — 0,9848 —0,1737 — 11,4301 —0,1737 — 5,5441 

180 —1,0000 0,0000 foe) 0,0000 lo-e) 
190 — 0,9848 OMT Si 11,4301 OATS 7 0,044] 
200 — 00,9397 0,3420 59,6713 0,3420 2,9142 
210 — 0,8660 0,5000 Seoul 0,5000 1,4331 
220 — 0,7660 0,6428 2,7475 0,6428 0,8814 
230 — 0,6428 0,7660 2,1445 0,7660 0,5799 
240 — 0,5000 0,8660 Irae 0,8660 0,4331 
250 — 0,3420 0,9397 1,4282 0,9397 0,3927 
260 —0,1737 0,9848 1,1918 0,9848 0,4249 
270 0,0000 1,0000 1,0000 1,0000 0,5000 
280 0,1737 0,9848 0,8391 0,9848 0,5905 
290 0,3420 0,9397 0,7002 0,9397 0,6715 
300 0,5000 0,8660 0,5774 0,8660 0,7217 
310 0,6428 0,7660 0,4663 0,7660 0,7256 
320 0,7660 0,6428 0,3640 0,6428 0,6744 
330 0,8660 0,5000 0,2680 0,5000 0,5670 
340 0,9397 0,3420 0,1763 0,3420 0,4096 
350 0,9848 0,1737 0,0875 0,1737 0,2148 
360 1,0000 0,0000 0,0000 0,0000 0,0000 

ue ut us 
: = 00 V 0,0000 0,0000 0,0000 
M 0,0000 0,0000 0,0000 
ees H 0,0000 0,0000 0,0000 

es, ie rt ies 
3,1416 1,5708 1,5708 


8. Zusammenfassung. Die Arbeit enthalt ein Verfahren zum Berechnen der reibungsfreien 
inkompressiblen Strémung durch gerade und radiale ebene Schaufelgitter im Rahmen der linearen 
Theorie. Behandelt wird das direkte Problem: Ermittlung der Druckverteilung, Umlenkung und 
Schaufelkraft aus Gittergeometrie und. Anstrémung. Als Hilfsmittel sind die bekannten Birn- 
baumschen Zirkulationsverteilungen des Plattenprofils auf gerade Plattengitter verallgemeinert 
und tabelliert. An den Beispielen eines geraden Gitters, dessen Geschwindigkeitsverteilung mit 
einem vorliegenden theoretischen Ergebnis verglichen wurde, und eines Kreisgitters, das berechnet 
und im Versuchsstand erprobt wurde, wird die Brauchbarkeit des Verfahrens gezeigt. In die 
Tabellen ist die konforme Abbildung des Plattengitters auf den Einheitskreis aufgenommen, wie 
sie zur Behandlung verschiedener Gitterprobleme niitzlich ist. 


(Eingegangen am 26. Oktober 1957.) 
Anschrift des Verfassers: Helmut Kriiger, Leverkusen-Wiesdorf, BreidenbachstraBe 22 
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Beitrag zur Berechnung von Einflufflachen schiefwinkliger Platten 
Von R. Heinen 


1. Einleitung. Die statischen Verhaltnisse schiefwinkliger Platten, die wegen ihrer vielfachen 
Verwendung besonders im Briickenbau groBe Bedeutung gewonnen haben, sind noch wenig geklart. 
Die meisten der in der Literatur vorhandenen Arbeiten beziehen sich auf schiefe Platten unter 
Gleichlast, wahrend nur wenige EinfluBflachen berechnet sind. Die zur Ermittlung der EinfluB- 
flachen rechteckiger Platten angewandten Methoden sind fiir schiefe Platten z. T. unbrauchbar oder 
erfordern fiir eine ausreichende Genauigkeit einen sehr groBen Rechenaufwand. 

In der vorliegenden Arbeit werden EinfluBflachen parallelogrammférmiger Platten aufgestellt 
nach dem Verfahren von F. Schultz-Grunow, mit welchem fiir andere Plattenformen bei verschie- 
denen Lagerungsarten und Aufpunktlagen mit einem ertraglichen Rechenaufwand Ergebnisse 
von hohem Genauigkeitsgrad erreicht wurden! * 3. Abb. 1 zeigt die im folgenden behandelte schief- 
winklige Platte. Sie hat die Form eines gleichseitigen Parallelogramms mit der Seitenlange 2a 
und den AuBenwinkeln a, = 27 und o, =ia. Sie ist ringsum voll eingespannt. Berechnet 
werden die BiegeeinfluBflache und die Momenteneinflubflachen m, und m, fiir die Biegemomente 
in Schnitten x =konst. und y=konst. fiir den Aufpunkt in Plattenmitte. 


2. Grundlagen. Die Differentialgleichung fiir die Durchbiegung w(x,y) einer Platte lautet 
unter Voraussetzung eines homogenen isotropen Baustoffes, eines konstanten Elastizitatsmoduls EF, 
einer konstanten Querkontraktionszahl y, der Giltigkeit des Hookeschen Gesetzes, einer im Vergleich 
zu den iibrigen Abmessungen kleinen und konstanten Plattendicke h und kleiner Formanderungen 


Adw = 2&9 | (1) 


Hierbei bedeutet p(x,y) eine beliebige 
stetig verteilte Belastung und 
N = Eh3/12(1—v?) 

die Plattensteifigkeit; x, y sind recht- 
winklige Koordinaten der Platten- 
ebene. Durch Differentiation der 
Durchbiegung werden nach _be- 
kannten Beziehungen die Schnitt- 
gréBen gefunden, z. B. die Biege- 
momente 


es Ow Cw \ 
emo | 
Cw Ow 
Abb, 1. — eee i ee 
m, N ay? y ont “ 


Je nach Lagerungsart der betrachteten Platte sind fiir jeden Punkt der Plattenrander zwei Bedin- 
gungen gestellt. Im Falle voller Einspannung des Randes mu8 die Durchbiegung am Rand und 
ihre Ableitung normal zum Rand null werden: 
w=—0 9 ai = 0 5 

wobei n die in die Normalenrichtung zum Rand weisende Koordinate bedeutet (Abb. 1). Die 
Lésung von (1) kann fiir beliebige Einzellasten P, und einer stetig verteilten Belastung p(x,y) mit | 
Hilfe Greenscher Funktionen G(x, yo; x, y) dargestellt werden, wobei xj, y, ebenfalls rechtwinklige | 
Koordinaten der Plattenebene sind: 


8 2 N w(x, ¥5 Xo %o) = 2 C(x, os x,y) P; + J JS G(%o. ¥o3 % ¥) p(x, y) dx dy . (3) 


1 F. Schultz-Grunow, Z. angew. Math. Mech. 33 (1953) S. 227. 

* A. Scherer, Allgemeine Herleitung singularer Lésungen der biharmonischen Gleichung, sowie Aufstellung 
Greenscher Funktionen am Beispiel einer Dreiecksplatte mit freiem Rand — Diss. Aachen 1955. 

° A. Bohning, Die Berechnung von EinfluBflachen der dreiseitig sowie der allseits eingespannten Quadrat- | 
platte unter Ausnutzung von Symmetrien, Diss. Aachen 1956. 
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Der Faktor 82N wird im folgenden der GréBe w zugeschlagen. Unter w ist daher weiterhin die 
8a N-fache Durchbiegung zu verstehen. Die BiegeeinfluBflache fiir den Punkt (%9, Yo) ist die durch 
eine Einzellast P = 1 im Punkt (x, yo) hervorgerufene Biegefliche. Fiir P = 1 geht (3) iiber in 


W(%> ¥5 Hor Yo) = G(X Yo % ¥) » (4) 
und auf Grund des Maxwellschen Vertauschungssatzes erhalt man aus (4) die EinfluBflache der 
Durchbiegung im Aufpunkt (x9, yo) fiir den Kraftangriff im Punkt (x, y) 


W(X» Yo XY) = G(x, ¥5 Xo» Yo) - (5) 
Die Funktion G setzt sich zusammen aus einem im Aufpunkt singularen Integral G)(x,y; %,79) 
und einem reguléren Integral G, = F,(x,y; %9,¥9). Die Funktionen F, sind so zu bestimmen, da} 
G den speziellen Randbedingungen des Problems geniigt. 

Durch Differentiation der BiegeeinfluBflache nach den Koordinaten des Aufpunktes (x, yo) 
erhalt man die EinfluBflachen der Schnittkrafte 1, Diese kénnen ebenfalls aus einem singularen 
und regularen Anteil zusammengesetzt werden, und fiir sie gelten die gleichenRandbedingungen wie 
fiir die BiegeeinfluBflachen. Die singularen Anteile der BiegemomenteneinfluBflachen kénnen durch 
Anwendung der Differentialoperation (2) auf den singuldren Teil G, der Biegeflache ermittelt 
werden, die regularen Anteile sind wieder so zu bestimmen, daf die speziellen Randbedingungen 
erfillt sind. Fur die KinfluBflache m, z. B. entsteht ein formal gleicher Ausdruck wie (5), wenn nun 
die Greensche Funktion mit M, bezeichnet wird: 


™.<(%o> Yos XY) = M(%, ¥3 Xo» Yo) + (6) 
Gleichung (6) liefert zu einer in (x, y) angreifenden Last P = 1 das Moment m, im Punkt (xo, yo). 


3. Die singularen Anteile der EinfluBflachen. Liegt der Aufpunkt (4, yo) im Inneren des Platten- 
bereiches, so ist? der singulare Anteil G, fiir die Biegeflache 


Gols 95 Xo» Yo) = P(r? In r) (7) 


r? = (x — %)? + (vy — 0)” 
Der Singularteil M,o fiir das Moment m,, leitet sich aus (2) ab: 


Cs aC, &G, 


mit 


(8) 


und der Singularteil M,, fiir das Moment m, ist 
2G, a 


aye |” dx8 


(9) 


Um von der Querkontraktionszahl y, die ja eine Materialkonstante ist, méglichst unabhangig zu 
sein, werden statt der Momente die Kriimmungen k,) = 0?G,/0x2 und ky. = 0G,/0y? ermittelt. 
Hiermit kénnen dann fiir jedes beliebige y mit (8) und (9) die Biegemomente ermittelt werden. 

Werden alle Langenkoordinaten mit der halben Seitenlange a dimensionslos gemacht und wieder 
mit den gleichen Buchstaben bezeichnet, wird also z. B. fiir x/a nun x und y/a nun y gesetzt, so 
erhalt man fiir den in Plattenmitte liegenden Aufpunkt x) = 0, yy = 0 


My = N( 


6 =@ +) 5+), (10) 
2 x2 
ep = In (9? +9) + 3 or i (11) 
9 v2 
hyo = In (x2 +9?) + a? saute (12) 
An den Randern interessieren die Normalableitungen dieser GréBen. Wegen 
LOSE CeO ey 
én. Ox On © Oy On 


ergibt sich bei den zugrunde gelegten Eckenwinkeln (Abb. 1) 
ae dy V3 


aa er aS 
1A. Pucher, Ing.-Arch. 12 (1941) S. 76. . 
2 Ph. Frank u. "R von Mises, Differential- und Integralgleichungen der Mechanik und Physik, Band I, 
S. 856ff. New York 1943. 
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Ho = F(x +39) [ine +¥*) +1], (13) 
te ogg + Baxt + /B—I OD “, 
Okyo 1 3 2 243 1 3 /3 2 15) 

in yy Se i) +3 of xy). ( 


4. Die reguliiren Anteile der Einflufflichen. Die Konstruktion des Regularteils erfolgt nach 
dem Verfahren von F’. Schultz-Grunow. Hiernach ist zunachst das Innere des Plattenbereichs kon- 
form auf das Innere des Einheitskreises abzubilden. Die in der Plattenebene (z-Ebene) definierten 
Fundamentalintegrale (10), (11), (12) und (13), (14), (15) weisen durch die konforme Abbildung 
auf die Kreisebene (¢-Ebene) in den Bildpunkten der Plattenecken Singularitaten auf, deren Ord- 
nung durch die AuBenwinkel der polygonalen Berandung festgelegt ist. Diese Singularitaten sind 
zuerst zu eliminieren. Dazu werden in der Ebene des Einheitskreises definierte Potentialfunktionen 
gebildet, die in den Ecken Singularitaten der gleichen Ordnung wie die Fundamentalintegrale be- 
sitzen. Bipotentialfunktionen erhalt man durch Multiplikation der so gewonnenen Potential- 
funktionen mit x, y oder r?. Die Funktionen und die Fundamentalintegrale werden dann in der 
Kreisebene in Potenzreihen nach @ entlang des Randes entwickelt, wobei 0 den Mittelpunkts- 
winkel des Einheitskreises bedeutet (Abb. 3). Da die Fundamentalintegrale urspriinglich fir die 
Plattenebene gegeben sind, miissen die Potenzreihenentwicklungen der Fundamentalintegrale mit 
Hilfe einer Beziehung s(Q), wobei s die Plattenrandkoordinate in der z-Ebene bedeutet (Abb. 1), die 
sich aus der konformen Abbildung ergibt, durchgefiihrt werden. Mit den Freiwerten der Funktionen 
werden so viele Glieder der Potenzreihen eliminiert, daB die Randwerte der Fundamentalintegrale 
nach Uberlagerung der Funktionen an den Ecken nicht mehr singular sind. Es miissen also die 
Koeffizienten der Reihenglieder mit den Exponenten < 1 beziehungsweise in der Normalableitung 
0()/én die Koeffizienten der Reihenglieder mit dem Exponent 0 ausgeglichen werden. Der Rand- 
verlauf der so ausgeglichenen Fundamentalintegrale hat in der z-Ebene in den Ecken den Wert Null 
und eine horizontale Tangente, da hier 0()/dn = 0 ist. Wegen der durch die konforme Abbildung 
hervorgerufenen Verzerrung ist aber auch die Randbedingung w = 0 in der Eckenumgebung nahezu 
exakt erfillt. 


Die im eben beschriebenen Ausgleich berechneten Anteile der Greenschen Funktionen, im fol- 
genden ,,ausgeglichene Fundamentalintegrale“ genannt, sind so zu erganzen, da die Randbedin- 
gungen fiir den tibrigen Randbereich erfiillt sind. Dies geschieht durch Hinzufiigen einer biharmoni- 
schen Funktionsreihe. Damit diese Funktionen das gleiche Verhalten aufweisen wie die ,,aus- 
geglichenen Fundamentalintegrale“, werden sie so konstruiert, daB ihre Randverlaufe in den 
Ecken den Wert Null und eine horizontale Tangente besitzen. Es kénnen hier Funktionen der ein- 
fachsten Art, wie spaiter gezeigt wird, herangezogen werden. Sie werden aus Potentialfunktionen, 
die in der ¢-Ebene definiert sind, hergestellt, indem aus ihnen durch Multiplikation mit x, y oder r? 
Bipotentialfunktionen gebildet werden. Die in der¢-Ebene gebildete Reihe von Potentialfunktionen 
und die aus ihnen gebildeten Bipotentialfunktionen weist gegeniiber den Funktionsreihen, die 

: ; unmittelbar in der Polygonebene (z-Ebene) angesetzt wer- 
sO) den kénnen, eine weit bessere Konvergenz auf, wie sich 

im Vergleich von auf beide Arten berechneten Beispielen 

zeigte'. Man kann bereits mit wenigen Funktionen, wenn 

man auBerdem noch die Symmetrien ausniitzt, Ergebnisse 

von hohem Genauigkeitsgrad erzielen. AuBerdem sind die 

g andwerte dieser Funktionen nur von einer Variablen @ 

abhangig im Gegensatz zu den entsprechenden von zwei 
Abb. 2. Koordinaten abhangigen Randwerten der Funktionen der 


z-Ebene. 

Erwahnt sei noch folgender Vorteil: Auf Grund der konformen Abbildung hat die Beziehung 
s(Q) etwa den in Abb. 2 gezeigten Verlauf. Der Anstieg von s(Q) in Eckennihe ist um so oroben 
je Kleiner der Innenwinkel der polygonalen Berandung ist. Da der Ausgleich der Eckensingulari- 
taten und die Uberlagerung der Funktionen des Regularteils in der ¢-Ebene erfolgt und da die 
Randverlaufe der ,ausgeglichenen F undamentalintegrale“‘ und der zu tiberlagernden Funktionen 


1 Siehe die Arbeit von F. Schultz-Grunow in FuGnote 1 von Seite 268. 
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in den Ecken den Wert Null und eine horizontale Tangente besitzen, ist die nahezu exakte Er- 
fillung der Randbedingung w = 0 in der z-Ebene in einem gréBeren Bereich der Eckenumgebung 
verifiziert als in der ¢-Ebene und zwar ist dieser Bereich um so gréBer, je spitzer die Ecke ist 
(Abb. 7 bis 12). Bei schiefwinkligen Platten bedeutet das einen nicht geringen Vorteil, da bei 
anderen Berechnungsverfahren gerade die Erfiillung der Randbedingungen in den spitzen Ecken 
Schwierigkeiten bereitet. 


Die oben beschriebenen Vorteile des Verfahrens von F’. Schultz-Grunow wiegen das Mehr an 
numerischer Rechnung fiir die konforme Abbildung und den Singularitiétenausgleich bei weitem auf. 


a) Konforme Abbildung. Fiir die konforme Abbildung des Innenbereiches eines Polygon- 
zuges mit den Eckpunkten 2, besteht die Schwarz-Christoffelsche Abbildungsformel 
dz C i el ae = 
aaa (C= Cr) ; 
Hier bedeutet « = x + ty oder z = re'*® den Ortsvektor in der Polygonebene, € = é + in oder 
€ =e e'® den Ortsvektor in der Ebene des Einheitskreises. Es sind Cx, die Randpunkte des Kreises, 
deren Bildpunkte in der z-Ebene Eckpunkte des Vielecks sind. Die ax sind die AuBenwinkel des 
Polygonzuges. Es muf das Integral 
g rn 
a eee 
z= G, | tx) 9 al + G (16) 


K 


berechnet werden. C,, C, und die untere Integrationsgrenze ¢, sind komplexe Konstanten, welche 
die Lage und GréBe des Vielecks in der Polygonebene festlegen. Es bereitet wesentliche Schwierig- 
keiten und einen umfangreichen Arbeitsaufwand, die konforme Abbildung fiir Polygone mit mehr 
als vier Ecken und mit beliebigen Winkeln x, anzugeben. Ist aber die Eckenzahl < 4 und a, ein 
rationales Vielfaches von z, so ist eine Lésung im ertraglichen Rahmen méglich. Abb. 3 zeigt die 


Zuordnung der Eckpunkte in der z- und in der ¢-Ebene. 


z-Ebene 


Cu 


Abb. 3. Zuordnung der Eckpunkte fiir die konforme Abbildung des gleichseitigen Parallelogramms. 


Wenn alle Langen mit der halben Seitenlange a dimensionslos gemacht werden und fir 2/a 
wieder z geschrieben wird, bestehen die Wertpaare (Abb. 3) 


a 0 oo == Of 
z,=/3>G=1, 


hes (oe ly TF, (17) 
4=—y3>t,——l1, 
i= 1 > (, = — 1. 


Diese ¢-Werte in (16) eingefiihrt, ergibt 
1p == (Co lhe ae eer 
+) ¢-yF e—n FP C+ E+ 9" 


On 
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ne ste 
1—x 


z=AF(,k)+ 8B 


Mit der Substitution 


erhalt man 


mit 
(1+ 3) 22—1 


ae 2B 
(=14 73) 2a | 


4 


und cos Y = 


Die Konstanten ergeben sich mit (17) zu 


sy : 
4 =}5 = 0,5418 9515 , B= 0; wobei K = F(%,, 4). 
Zur praktischen Durchfiihrung der Rechnung wird der Plattenbereich mit einem Gitternetz 
iiberzogen, dessen Punkte in Abb. 4 angegeben sind. Die Bilder dieser Punkte in der Ebene des 


Einheitskreises werden mit Hilfe folgender Beziehungen berechnet: 


Hieraus ergibt sich 


(18) 


= a a ane 5 
(73 4 1) (V3 ena] —li+en 4] 
Die Trennung von ¢ in Real- und Imaginarteil liefert die Koordinaten der Bildpunkte 


Ere (GC). 7) Lin) 


Zur Berechnung der Jakobischen Funktionen muf zunachst eine Hilfstafel nach einem von Le- 
gendre' angegebenen Verfahren aufgestellt werden. Hiernach kénnen die Argumente yp der 
Beziehung 


v F(p, k) = Fy, k) (v< 1) 
berechnet werden. Die Gitternetzeinteilung ist so gewahlt, daB die x- und y-Achse in 10 Teile 


ees ist. Den Gitterpunkten auf der x-Achse, die /3-mal langer ist als die y-Achse, entsprechen 
ie Werte 


y¥3 =v 4 F(x,k) mit i male (» = 0; 0,170.2; ., -0.9; 10) 


und den Gitterpunkten auf der y-Achse die Werte 


y-l=yAF(y,k) mit Kk? = 
Fir . 
v F(a, k) = Fy, k) 
und 
y Fg, k y= Fe & ) 

sind die Winkel y bzw.y’ in Tabelle 1 zusammengestellt. Mit ihnen werden die erforderlichen Werte 
der Jakobischen Funktionen ermittelt. Tabelle 2 und Abb. 4 enthalten die konforme Abbildung 
der Gitternetzpunkte. 

F ir die numerische Auswertung der Regularteile der Greenschen Funktion wird zweckmaBig 
noch eine spezielle Abbildungsrelation zwischen der Randkoordinate s in der z-Ebene und der 
Randkoordinate © des Hinheitskreises gebildet. Es wird wieder die Schwarz-Christoffelsche 


1 A. M. Legendre, Traite des Fonctions elliptiques et des Integrals Euleriennes, Paris 1826. 
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Punkt é n Punkt é 
LT 
00 0 0 33. i] 
01 0,18146272 0 34. 
02 0,35848496 0 35 
03 0,52527544 0 36 
04. 0,67445035 0 37 
05 0,79837599 0 40 0 
06 0,89156250 0 4l 
07 0,95279035 0 42 
08 0,98577690 0 43 
09 0,99816863 0 44 
010 1,00000000 0 45 
10 0 0,10529195 46 
1l 0,18217456 0,10411422 50 0 
12 0,36015601 0,10011060 51 
13 0,52820351 009223326 52 
14 0,67869903 0,07975142 53 
15 0,80358271 0,06314913 54 
16 0,89698526 | 0,04440833 55 
17 0,95757299 0,02642072 60 0 
18 0,98924241 0,01193116 61 
19 0,99999620 0,00273477 62 
20 0 0,21133037 63 
21 0,18421954 020914589 64 
22 0,36505833 0,20152405 710 0 
23 0,53695845 0,18604996 a, 
24 0,69156774 0,16093722 712 
25 0,81946127 0,12704764 73 
26 0,91355094 0,08855840 80 0 
27 0,97215179 004855630 81 
28 0,99975683 0,02205173 82 
30 0 0,31873534 90 0 
31 0,18730242 0,31593383 91 
32 0,37282283 0,30557756 100 0 


Tabelle 1 
y y y’ 
0,1 18,29289929° 10,51769762° 
0,2 36,47289112° 20,71302042° 
0,3 54,47162633° 30,31324833° 
0,4 72,29085279° 39,12810334° 
0,5 90,00000000° 47,05859617° 
0,6 107,70914721° 54,08619138° 
0,7 125,52837367° 60,25206774° 
0,8 143,52710888° 65,63530297° 
0,9 161,70710071° 70,33498191° 
1,0 180,00000000° 74,45773676° 


Tabelle 2. Konforme Abbildung der Gitternetzpunkte 


Gleichung zugrunde gelegt: 


Auf dem Rand ist € = e'? und ds = {dz| . 


Hiermit erhalt man nach Umformung aus (19) 


s = C, f (sin? O cos @)— 13 dO + C,. 


Mit der Substitution 


wird 


) in |) 
= are sin oe ae 


s=AF(,k) + B:; 


Fx GM — 1+ YI + Gy. 


0,55142670 
0,71343006 
0,84682220 
0,94218985 
0,99722201 


0,19084110 
0,38270469 
0,56888371 
0,74492738 
0,88714460 
0,98458900 


0,19382632 
0,39336035 
0,59602893 
0,78699694 
0,94280904 


0,19457772 
040247647 
0,62425774 
0,84090505 


0,19037977 
0,40619312 
0,65350208 


0,17692554 
0,39827448 


0,14749316 


0,28318553 
0,24518384 
0,19255962 
0,13217270 
0,07431732 
0,42784838 
0,42513281 
0,41367476 
0,38541682 
0,33454543 
0,26080151 
0,17488427 
0,53850128 
0,53706215 
0,52730980 
0,49624902 
0,43184046 
0,33333333 
0,64967116 
0,65145491 
0,64806861 
0,61835183 
0,54118268 
0,73982491 
0,76707991 
0,77757082 
0,75692472 
0,84890516 
0,88120481 
0,91726629 
0,94080430 
0,98906308 
1,00000000 


(20) 
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Tabelle 3. Konforme Abbildung der Randpunkte s(@) 


0 
0,29365860 
0,36998993 
0,50218551 
0,63285768 
0,79807712 
0,91496782 
1,00923342 
1,09023704. 
1,16263038 
1,22914911 
1,29159991 
1,35129380 


Abb. 4, Konforme Abbildung der Gitternetzpunkte. 


1,40926983 
1,46643030 
152364183 
1,58183538 
1,64214544 
1,70616862 
1,77660551 
1,85958427 
1,91159281 
1,95201266 
1,96977055 
2,00000000 


30 40 
Abb. 5. Funktionverlauf s(@). 


50 


70 80 


12 


XXVI. Band 1958 RR. Heinen: Beitrag zur Berechnung von EinfluBflachen schiefwinkliger Platten AES 


Der Modul k und das Argument gy ergeben sich zu 


24/3 
eo, 


een Ea yee 


~ = arc cos oe 


1+y3)v+1 


und die Konstanten A und B zu 
2 


A =——_ = 0,5418 9515 , 1 =D) 
F(z, k) 


Tabelle 3 und Abb. 5 zeigen die Zuordnung der Randpunkte in der z- und ¢-Ebene. 


Wie oben erwahnt, wird fiir den Singularitatenausgleich noch die Reihenentwicklung s(@) 
bendtigt. In der Bezichung (20) kann fiir die Umgebung der Ecke sin @ bzw. cos @ durch das 
erste Glied der Reihenentwicklung ersetzt werden, da lediglich die Koeffizienten der Reihen- 
glieder mit den Exponenten < 1 eliminiert werden. Fiir die Ecke I (Abb. 1) bei O = 0 ist fiir 
die Randkoordinate s, 

" dO 
UB 


=3C,08=K, 0%, 


Se 


Abb, 6. Zuordnung der Eckpunkte fiir die konforme Abbildung des ungleichseitigen Parallelogramms. 


Fiir die Ecke II bei 9 = x7/2 ist fiir die Randkoordinate s, 
d@ 3 | r 
= Cs | Sip 7 Ca OF = KO. 


Die Konstanten K, und K, ergeben sich zu 


ey Sy a ey EW Ip 
Der Vollstandigkeit halber sei noch die Zuriickfiihrung des Abbildungsintegrals fiir ein ungleich- 
seitiges Parallelogramm auf die Legendresche Normalform angegeben. Abb. 6 zeigt die Zuordnung 
einiger Randpunkte in der z- und ¢-Ebene. Es entsprechen sich 
ine > 65 = 05 
Pr faces oat a 
i, =) > Co = h, (21) 
Z,=—a—>l(,=—l, 
4 =—b>f,=-—), 


wobei die GréBen 
aX, b=%, +12, 


A = 4 +i, = cos O, + isin O, 
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vom Seitenverhiltnis des Parallelogramms abhangige komplexe Konstanten sind. Somit lautet.das 
Abbildungsintegral f 
ZG Es . 
| (2—1)'8 6—ayh E+ 917 6 +4)" 
und nach Umformung | 


2=G fIG—) C—PP- + G. 


eer oA 
pa fee) 


+ Cy 


Mit der Substitution 


wird 


Setzt man 


gh Nh . (23) 


so erhalt man 


(24) 


Das Integral (24) hat die Form 


du | ; du 
| YP? J V(a—u) (b—u) (c—u) (d—u) 
Es 1aBt sich! auf die Legendresche kanonische Form transformieren, und zwar auf ein elliptisches 
Integral erster Gattung 


Hier wird also (24) zu 
z=A [Fi ky) + F (Pz; k,)} + B. 


Die Wurzeln a, b, c, d der Polynome in den Integranden sind alle reell. Sie berechnen sich aus (24) zu 


a=+2, b= 2 cos, 
aes O, | % vA 0, Tt 
c= —2e0s(3 +5), d = —2e08(32—F). 


Die Konstanten A, B, das Argument ~ und der Modul k sind abhangig vom Eckenwinkel 0, 
(Abb. 6). Uber die zu (22) und (23) inverse Beziehung 
, FIC + 12» 
4 @—%) C1) 


und aus der Zuordnung der Randpunkte (21) kénnen der Winkel O, und damit der Modul k und 
die Konstanten A und B bestimmt werden. 


b) Der Ausgleich der Eckensingularitaten. Der Aufbau des Regularteils gliedert sich, 
wie erwahnt, in zwei Teile. Im ersten Teil, dem Singularitatenausgleich, werden die durch die 
konforme Abb. in den Eckpunkten auftretenden Singularitaéten ausgeglichen. Aufgrund der Ab- 
bildungsformel haben die zum Ausgleich zu benutzenden Funktionen zweckmaBig die Form 


TAC) Ce ere 


Fiir die numerische Auswertung werden die Real- baw. die Imaginarteile dieser Funktionen heran- 


gezogen: 
Re ’ cos n/3 0; yp 
C—C no G n/3 te Kn , 
Im 3 i ae | fe n/3 A tag 


1 P. F. Byrd u. D. Friedmann, Handbook of Elliptic Integrals for Engeneers and Physicists, Berlin- 
Gottingen-Heidelberg 1954. 
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Beim Durchgang durch die Ecke K wird (¢ —¢,) = 0 und das Argument @, erfahrt einen Sprung 
nm —z. Deshalb ist der Funktionsverlauf in verschiedene Definitionsbereiche zu unterteilen. 


Auf dem Kreisrand fiir 9 = 1 gilt dann der folgende Verlauf der Funktionen von K = 1 bis 
K = 4 gemaB Tabelle 4. 


Tabelle 4 
ee ee a ee 


Definitionsbereich 9<O0<5 5 <O<n n<O<—on an <O<2n 

Funktion | Betrag Argument 0; Argument 0; Argument 0; Argument @; 
_ n . O\n/3 © .& OQ ke (C) 2) 
€ —C. [3 Pin 2 ae \ TU 7 
G—oy) Yin ae aT 5 alae aa) a7 3 
¢—t,)"8 Pon a OD aN PN 0 1% (OQ! a OQ Sa O on 

Yon 2 4 | 2 4 2 4 2 4 2 4 

@ OL aN ini3 C) oO (C) C) 
a ea) TEN pea = ee —— a 
Cs) Ys n De aeanay 2 2 aes Buin 

+ Ky WS Oe aa Qs i QD. a C) 3 
ee [3 Pan 9 so re Es: | 
eet V4 n sl ae | 2 oA meg ar 2 he 


Fir die Erfillung der Randbedingungen interessieren auBerdem die Normalableitungen dieser 


Funktionen f(¢) am Rand 


FO) _ ae) a a 


on at G0 an* 
Hier ist 00/0n = |d¢/dz| der Ma®stabsfaktor der konformen Abbildung. Aus Gleichung (20) folgt 
patos (sin? O cos O)1/3 . 
2/3 A 


Die Normalableitung der Funktionen f(¢) lautet also 
Hee MoE ay aig 
on 


3 
Diese Funktionen haben noch nicht die Symmetrien der Fundamentalintegrale (10) bis (15). Sie . 
miissen Symmetrie zur x- und y-Achse besitzen. 


d¢ 
dz 


¢—2g(7) 


Die Symmetrieeigenschaften finden Ausdruck in den Bezichungen 
f(O) =f — 0) =f + 0) =fex—®). 
Bildet man folgende Kombinationen, die sich aus den Realteilen gx,, und den Imaginarteilen yx, der 
Funktionen zusammensetzen, so erhalt man Funktionen, die den Symmetriebedingungen geniigen: 
8&1 = Pu + 3 Yu + 2 Pay » 
82 = 3 Pa — Ya + V3 Pas + Par » 
83 = — Pet /3 We + 2 Pao » 
b4 — F227 3 P22 T Pag 1 3 Pax > 
os — aa /3 Yu — 2 Psa » 
f6 — Yea a /3 Vos “nA oe /3 Yas ? 
87 = 15 /3 Ys + 2 P35 » 
Cao" 


Es handelt sich hier um Potentialfunktionen. Man braucht aus ihnen keine Bipotentialfunktionen 
zu konstruieren, wenn man die Fundamentalintegrale ebenfalls in Potentialfunktionen U, V 
zerlegt, etwa in der Form 


G=U+rv. 


Dies gilt sowohl fiir die Durchbiegung wie fiir die Kriimmungen. Die aus (10) bis (15) gewonnenen 
Werte U und V der Fundamentalintegrale, die kurz ,,FundamentalgréBen“ genannt werden sollen, 
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2.0% 


U,,9 =mr+2, 


Ue nr a 


weeCo cy 
on 
OU ko mee V3 
Ge Mare? w 
CUiyy _ 
on re? 


Welche Funktionen g sind zum Ausgleich dieser FundamentalgréBen heranzuzieben ? Die Ord- 
nung der durch die konforme Abbildung entstandenen Singularitaten ist fiir die Ecke I (bei 0 = 0) 
eine andere als fiir die Ecke II (bei 0 = 7/2), da sie durch den AuBenwinkel der polygonalen 
Umrandung festgelegt ist. So haben die in der Eckenumgebung nach @ entwickelten Fundamental- | 
gréBen an Ecke I die Potenzen £n (n = 1, 2,3) und an Ecke II die Potenzen =n. 
gleich der FundamentalgréBen an Ecke I werden die aus (¢ —€,)"* und (€ —€,)"* gebildeten 
Funktionen, namlich g,, g5,g,... herangezogen, zum Ausgleich an Ecke IT die Funktionen gy, g¢, g,..., 
die aus (€ —€,)"? und (€ —€,)"° gebildet sind. Die Funktion g, entfallt, da an Ecke II keine 
Potenz 1/3 auftritt. Wie erwahnt, sollen die Koeffizienten der Reihenglieder mit den Potenzen } 
n <3 in der Funktion und mit der Potenz n = 0 in der Normalableitung eliminiert werden. Eine | 
Aufstellung der Reihenkoeffizienten aller FundamentalgréBen und der zu iiberlagernden Funktionen 
findet sich in der Tabelle 5. Hier ist zu bemerken, daB zwischen den K oeffizienten lineare Abhangig- 
keiten auftreten, die schon an anderer Stelle! begriindet wurden. Abhangig werden die Glei- 


chungen Nr. 2 und 4, 3 und 5, 8 und 9. 


Tabelle 5. Koeffizienten der Funktionen und der Fundamentalintegrale fiir den Singularitdtenausgleich 


Se V3 2y + 3xy2+ V3 ¥3 
: i soma - 7) 


r6 


Ving —x8—3 3 y+3xy2+ 34 
a — ee 


Zum Aus- 


cy Cy C3 C4 cS 6 Gy kyo» kyo 
Nr. 
&1 §3 84 &s 86 &s U i“ U Vv | 
3 ya= 3 — —~ = SS eee aaa 
1 @ 24 } 2 ya! ys |—y27 |_y28 |-_ you y/2_y3] 0 = In3+2 _ 
hal Ai oe y3 0 0 0 0 0 0 | =2/3 4 doors + 4 
| 
3 @7/3 0 1 0 0 0 0 0 an 24 | a 
/3 /3 y | 
1 
4 Ce ee 0 0 0 0 0 sete = y 
oF, 23 A 2y3 y3 313 
Sa ia | 
5 Qs ae z a | 
0 a 0 0 0 0 0) 4 | ~ A |——A] 
3 3 3 
| 6 Fanta 3,— _ — ja fee ; io. 
6 or |jary2+y3 | }2# yB ye |—j35|—yar] —jae | 0 2 7 
| os 
en Re 0 0 0 0 0 0 0 0 0 0 
a | 
8 2/3 3 
C) 0 0 i 0 0 0 0 —lA —\3 A |—y3 A | 
oF 13 : = | 
PIA @e 3 = | 
se || fo 0 0 sl i 0 0 0 0 1 y3 y3 | 


1 Siehe die Arbeit Bohning in FuBnote 3 von Seite 268. 
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Dadurch wird das Gleichungssystem auf 5 Gleichungen reduziert. Die Koeffizienten dieses 
Systems ergeben die Tabelle 6, und die daraus errechneten Randverlaufe der so ausgeglichenen 
Fundamentalintegrale sind in Abb. 7 bis 12 als Kurven mit der Bezeichnung II aufgetragen. 


Tabelle 6 
Nr. Funktion eh Up baw. Uns 9 View g bev —Vivg 
: gy | +0,4117 5127 + 0,8235 0255 —0,2745 0085 
- & 4) -+-0,1695 3911 + 0,3390 7822 —0,3390 7822 
£4 +0,3565 8706 + 0,7131 7413 + 0,7131 7413 
: gs | + 0,4530 6987 +1,1496 1207 + 0,8433 8273 
& | + 0,4403 7411 +1,1242 2054 —0,8103 2129 
4am 
| 2 on 
|e 
I Fundamentalintegrale | Ve ~ SAS 
II Singularititenausgleich ‘ WA : 


It Punktweise Ertiillung 


I (20-fach) 


S 


E L 
05 40 45 40 0 05 10 15 20 
Abb. 7. Abb. 8. 


Abb. 7—12. Randverlaufe der Fundamentalintegrale (I), nach dem Ausgleich der Eckensingularitaten (II), nach der punktweisen Erfiillung 
der Randbedingungen (II). 


c) Punktweise Erfiillung der Randbedingungen fiir den Rest des Randbereiches. 
Die so berechneten regularen Anteile der Greenschen Funktionen sind so zu erginzen, daB die 
Randbedingungen in den iibrigen Randteilen, und zwar fiir eine gewisse Anzahl von p einzelnen 
Randpunkten, erfiillt sind. Es werden in der Kreisebene definierte Potentialfunktionen der Art 


Vor 
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(¢) =" verwendet. Den Symmetrieverhaltnissen entsprechen die Funktionen 
f, = Re (22") = 9?™- cos2m@ (i == 082, Sea) s 


Die zugehérigen den Symmetrien geniigenden Bipotentialfunktionen uw,, erhalt man durch Mul- 
tiplikation mit r? 


Uo Toda 
ky 
mM 
\ 
S 
ae 4) 
= a 
SS 
os 
~S 
ye ee = = = 
BATS 
Sel re 
> 
1 es Bi le ey To a cae 
eg ante 
ea iim. 
ae | 
OS 
I (10-fach) 
[ @ 
0 10 20 50 40 50 60 70 80 40 
2 
¥ SSS 
as, 
Se 
os 
SS 
SS 
g =e 
“ SS 
SS 
™S 
SE 
g SS 
NS 
mae 
™~ 
7 | 
toy 
ete - 
; — ——— Pe 
It (10-fach) 
Nite eeceme es eh 
0 O85 40 i) 20 
Abb, 9. 


Die Funktionen Fe (¢?™~+1) erfiillen nicht die Symmetriebedingungen. Bildet man aber aus ihnen 


Bipotentialfunktionen der Art 
«© Re (cars 2) und a Hal Cee ne) 
so sind die gestellten Bedingungen erfullt. Folgende Funktionen kénnen also verwendet werden: 
jn 0. cosa, 
u,, = 1? 0°” cos 2mO> 


V_ = xO?™+) cos(2m+1)G, 


w, = y¥ 02" +) sin (2 m + le. 
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Die Normalableitungen dieser Funktionen ergeben sich, wenn lediglich die Randwerte (or=211) 
der Normalableitungen angegeben sind, zu 


2 P= 2m cos 2mO[E |, 


cus nO toon 


alg /3 |) 


om — eos (2m +1) 6 


OWm 
=sinQm+)6lyQm+1)[2) +3]. 


rQm+y/#\ 43], 


on 


Abb. 10. 


Wie erwahnt, miissen diese Funktionen, ehe sie den ,,ausgeglichenen F undamentalintegralen“ 

tberlagert werden, in den Ecken im Wert und in der ersten Ableitung zu Null werden, damit sie 

das gleiche Verhalten wie die ,,ausgeglichenen Fundamentalintegrale“ aufweisen. Sie werden daher 

in Potenzreihen nach. @ entwickelt und entsprechend den Fundamentalintegralen ausgeglichen. 
19 
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Werden nur die zum Aus 
bezeichnet der Index I bz 


hip ec alo % 


tim =3— 613 AQUS 4..., 


— 4 
fin =|3 ls AO ane 


Wim =O-+-+-, 


Stim =a ee 


Witm = + (1—13 A OF 1... | 


VITm =-+ (0 = 7 


Ne 


gleich notwendigen Glieder der Reihenentwicklungen angegeben, und 
w. II die Entwicklung an Ecke I bzw. Ecke II, so ergibt sich 


Sint 

pet at ie ; 
OU, 1 a 

a = 3 eT, Boxe 
OV, 1 a ih 

an gabe 
Ow, 1 a 

on i | “) 
Pn 
Fe =O, 
Un 1 aS 


Fir die Reihenentwicklun 


gen an Kcke II ist fiir 
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Exponenten 0 an beiden Ecken in der Funktion und in der N ormalableitung eliminiert, so sind gleich- 
zeitig die Koeffizienten der Reihenglieder mit dem Exponent 1/3 an Ecke I bzw. 2/3 an Ecke II 
in den Funktionen eliminiert, d. h. die gestellte Bedingung, da® der Randverlauf der Funktionen 
in den Ecken den Wert Null und eine horizontale Tangente haben sollen, ist erfiillt. 


“Gy 05 10 
Abb. 12. 


Dieser Ausgleich wird gleichzeitig mit der Erfiillung der Randbedingungen in einzelnen Rand- 
punkten, und zwar fiir die Punkte 


GO =415 7, A5e2 15" finaG; bund. —k 
Oa, 45", To” fiir 


durchgefiihrt. Es ergeben sich die Koeffizienten von Tabelle 7 fiir das Gleichungssystem der 
BiegeeinfluBflache. 


Tabelle 7 
INES Funktion Nr. Funktion 
1 f, —0,4396 3120 6 uy —0,3177 3964 
2 Ug + 0,1816 2121 7 v, —0,2072 3731 
3 v% + 0,8385 0626 8 w, + 0,8776 9985 
4 Wo —0,0287 9557 9 Vp + 0,3929 0754 
5 fi —0,4999 7146 10 Uy —0,0461 1314 


De 
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Die auf diese Weise erreichte Naherung ist bereits ausreichend. Die grote Randabweichung 
fiir die Biegeflache z. B. betragt noch 0,0341. Der Durchbiegungswert im Aufpunkt ergibt sich zu 
0.43963. Durch Uberlagerung der Funktion v, konnten die Randabweichungen allgemein ver- 
kleinert werden. Man erhalt nunmehr fiir die Biegeflache als gréBte Randabweichung 0,00871 und 
fiir die Durchbiegung im Aufpunkt den Wert 0,43943. Damit ist die gute Konvergenz des Ver- 
fahrens bewiesen. Weiterhin zeigt sich, daB die Randabweichungen sich auf die groBen mab- 
gebenden Ordinaten der HinfluBflache im Platteninneren nur in geringem Mae auswirken. 

Nach Uberlagerung der Funktion v, ergeben sich die Konstanten von Tabelle 8 und die nach 
vorhandenen Randabweichungen sind in den Abb. 7 und 8 als Kurven mit der Bezeichnung III 
dargestellt. 


Abb. 13. BiegeeinfluBflache. 


Tabelle 8 

Nr. Funktion Nr. Funktion 
| | 

1 bar oe —0,4394 3128 7 v —0,3897 8550 
2 Up | 00380 4503 8 Ww, + 0,2982 0654 
3 v +0,4799 6729 9 fr + 0,2620 3212 
4 Wo +0,7472 6331 10 Us | —0,0959 3079 
5 fi | 40,0471 2918 11 vy 40,0602 4096 
6 uy +0,0905 5258 

Tabelle 9 
Nr. | Funktion Nr. Funktion 
1 camer —0,9009 9424 7 v | 2,0211 

| E [Reese 1906 

2 u |  —0,0120 2558 8 i, — 0.5359 5679 
3 vp + 1,7308 4736 9 ve +0,5928 3854 
4 b matt +1,1287 2467 10 Uy —0,3681 4183 
5 eee +0,5241 9066 11 vy +0,0408 1633 
6 | uw |  +0,4521 7910 | 


Nachdem die einzelnen Funktionen fiir verschiedene Punkte de i i i 
die r Kreisebene ermittelt d, 
lassen sich die EinfluBflachen fiir die Punkte berechnen, und sie kénnen dann in die Bildpunkte 
der Plattenebene iibertragen werden. Die so ermittelte BiegeeinfluBflache ist in Abb. 13 dargestellt. 


Die Durchbiegung in Plattenmitte unter der Einwirkung einer Einzellast P ergibt sich hier- 
aus zu 


~ Pa@ 
w = 0,0175 Ww? 
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die Durchbiegung unter Kinwirkung einer stetig verteilten konstanten Last p ergibt sich zu 
at 
w = 0,0124 7, 


Zur Erfiillung der Randbedingungen fiir die EinfluBflachen der Kriimmungen k, und k, kénnen 
die gleichen Funktionen des Gleichungssystems der BiegeeinfluBflache benutzt werden. Es andern 
sich lediglich die auf der rechten Seite stehenden ausgeglichenen Fundamentalintegrale. Hier 
konnten ebenfalls durch Uberlagerung der Funktion v, die Randabweichungen verkleinert werden. 
Die gréBten Abweichungen ergeben sich fiir die k,,-EinfluBflache zu 0,00972 und fiir die k,-HinfluB- 
flache zu 0,01834. Der Verlauf der Randabweichungen findet sich in den Abb. 9 bis 12 als Kurven 


Abb. 14. m,-EinfluBflache. 


mit der Bezeichnung III. Die zurk,- EinfluBflache gehérenden Konstanten erhalten die Werte von 
Tabelle 9 und die zur k,-EinfluBflache gehérenden die Werte von Tabelle 10. In Abb. 14 und 15 
sind die EinfluBflachen der Momente m, und m, fiir eine Querkontraktionszahl von y = 0,3 dar- 
gestellt. Tabelle 11 enthalt die Ordinaten der drei berechneten EinfluBflachen. 


Das Moment m, unter Einwirkung einer stetig verteilten konstanten Last p ergibt sich (fir 
r = 0, Seitenlange a) zu 


4 
m, = 0,0084 7 
und das Moment m, Zu 


= pa 
m, = 0,0183 7. 


Tabelle 10 

Nr. Funktion | Nr. | Funktion 

1 h —0,6148 4136 7 i e ie. 20,2815" 7877 
2 Up +0,1456 4245 8 | Ww, +2,2118 5244 
3 vp +.0,7473 5758 9 Ne. of |  +0,4479 3770 
4 Wo + 0,6907 2122 10 | Ug + 0,0522 4480 
5 hi —— —0,9274 6943 ll vy + 0,2347 4187 
6 Uy —0,5626 7831 
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Tabelle 11. Ordinaten der EinfluBfléchen (8 N-fach) 


Punkt w | x ly 
00 + 0,43943 oo) oo) 
01 + 0,37661 EPG 3.070115 
02 + 0,27089 | OS AED koisvarl 
03 + 0,16761 | + 0,04416 + 1,284.98 
04 -- 0,08670 — 0,23451 | + 0,75741 
05 + 0,03524 — 0,284.29 | + 0,44068 
06 + 0,00779 082085 BIa()-977230 
07 + 0,00363 — 0,13098 + 0,12888 
08 | + 0,00265 — 0,04089 + 0,06285 
09 | + 0,00112 — 0,00866 | + 0,00568 
010 0 0 0 
20 + 0,35438 + 3,30604 + 2,15582 
21 + 0,31410 ak wl + 2.03705 
22 + 0,22931 | + 0,73804 = 1.52357 
23 + 0,13938 + 0,07207 + 0,98810 
24 + 0,06836 — 0,19390 + 0,59771 
25 + 0,02514 == (0223958 | ++ 0,34278 
26 + 0,00410 | — 0,17466 + 0,19565 
DT + 0,00239 | — 0,05691 | + 0,08719 
28 + 0,00308 | — 0,00208 — 0,00069 
40 | + 0,21393 | + 1,58587 + 0,59155 
Al + 0,19050 + 1,20762 + 0,61731 
42 + 0,13420 + 0,53781 2 0.57714 
43 -+ 0,08333 -+ 0,02827 + 0,39322 
4A + 0,02672 —0,10790 + 0,25314 
45 + 0,00413 — 0,09180 + 0,11417 
46 | + 0,004.92 — 0,00405 | — 0,00539 
60 + 0,08822 + 0,70400 | — 0,03458 
61 + 0,07354 + 0,55014 | + 0,00129 
62 + 0,03124. + 0,29267 + 0,09223 
63 + 0,01108 + 0,01206 + 0,05263 
= _ (of — 0,00545 ++ 0,00537 — 0,00143 
80 + 0,02533 + 0,19673 — 0,18849 
81 + 0,00935 ++ 0,14000 — 0,10472 
82 — 0,00486 + 0,00777 + 0,01623 
100 | 0 0 0 


Abb. 15. my-EinfluBfliche. 
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Graudenz! erhilt hier die Werte 
m, = 0,0076 P* 
bzw. 
ae pat 
my = 0,0191 NO" 


Die Anregung zum vorliegenden Aufsatz und wertvolle Hinweise bei der Ausarbeitung gab 
Herr Professor Dr. F’. Schultz-Grunow, dem ich an dieser Stelle meinen ganz besonderen Dank aus- 
sprechen méchte. 


Aus dem Institut fiir Mechanik der Technischen Hochschule Aachen. 
(Eingegangen am 11, November 1957.) 


Anschrift des Verfassers: Dr.-Ing. Richard Heinen, Disseldorf, Kurfiirstenstr. 10. 


1 H. Graudenz, Beitrag zur Berechnung der MomenteneinfluBfelder schiefwinkliger Platten, Diss. Hanno- 
ver 1948. 
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Beitrag zur Berechnung von Translationsschalen* 


Von G. Mittelmann 


1. Allgemeines. Eine Schale laBt sich bekanntlich als zweidimensionales Kontinuum auffassen, 
und dementsprechend werden auch der Spannungs- und Verformungszustand zweidimensional dar- 
gestellt. Mit Hilfe der Tensorrechnung in der Schreibweise des Ricci-Kalkiils, ohne den man heute 
bei der Untersuchung solcher Probleme nicht mehr auskommt, ist es méglich, die Gleichgewichts- 
bedingungen und die Spannungs- und Verzerrungsgleichungen in allgemeiner Form exakt auf- 
zustellen. Die praktische Lésung wird aber auf diesem Wege auBerst schwierig und ist nicht durch- 
fiihrbar. Es miissen daher Naherungen eingefiihrt 
werden. Dies gilt auch fiir die Berechnung der 
Translationsschale. 

Es werde ein im allgemeinen schiefwinkliges Ko- 
ordinatennetz in die Mittelflache M der Schale S 
gelegt. Mit den beiden krummlinigen Koordinaten 
Q, und @,, sowie der Koordinate @, in Richtung 
der Normalen zur Schalenmittelflache M 1aBt 
sich der Ortsvektor Si eines beliebigen Schalen- 
punktes A wie folgt ausdriicken (Abb. 1): 

yt = ter a (1) 
Hierbei sind = t6,, o,) der Ortsvektor eines Punk- 
tes der Schalenmittelfliche M, ferner z = tO, 
der Abstand eines Punktes A der Schale S von 
der Schalenmittelflache M in Richtung der Nor- 
malen (9, lauft demnach von — 0,5 bis + 0,5). 
Weiter ist t die Dicke der Schale S in jedem Punkt, ebenfalls in Richtung der Normalen und 
(3 = A3(0,,0,) der Normalenvektor der Schalenmittelflache M. 

Die metrischen FundamentalgréBen (kovariant) g;,, der Schale S lautent: 


Sap yp —2tO, by 6 ar (t)? (O3)? Cx B 2 
8x3 = 9, (2) 
&33 = (i). 


Darin bedeuten a, Koeffizienten der ersten Fundamentalform der Mittelflache M oder ko- 
variante MaBzahlen des MaB tensors, b, , Koeffizienten der zweiten Fundamentalform der Mittel- 
fliche M oder kovariante MaBzahlen des Kriimmungstensors, c,, Koeffizienten der dritten 
Fundamentalform der Mittelflache M. 

Ks interessiert zundchst die Geometrie der Schalenmittelflache M (zweidimensional). Die oben 
angegebenen Koeffizienten der entsprechenden Fundamentalformen lassen sich mit Hilfe des Orts- 
vektors t fiir die Mittelflache wie folgt berechnen: Mit Hilfe der kovarianten Basisvektoren der 
Mittelflache 


Abb. 1. Geometrische Beziehungen. 


or 
a, = 00x (3) 
erhalt man 
ay gp = Ay > Ag - (4) 
AuBerdem gilt 
oa OA. 
bus = — On 26, = + Gq, (5) 


* Dissertation der Technischen Hochschule Hannover. Referent: Professor Dr.-Ing. Zerna; Korreferent 
Professor Dr.-Ing. Pfliiger. 

1 Es werde vereinbart, daB lateinische Indizes die Zahlen 1, 2, 3 und griechische Indizes die Zahlen 1, 2 
verkorpern. Weiter soll tiber doppelt auftretende gleiche Indizes summiert werden. Potenzexponenten werden 
durch Klammern abgetrennt, um eine Verwechslung mit obenstehenden Indizes zu vermeiden. 
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wobei 
1 
Use ie (a, X Mg) (6) 
mit 
@ = Ay gg — (419)? > 0 (7) 


ist. Von den kontravarianten MaBzahlen werden bendbtigt 
ace Ay a 
qi — 2 | p= Ca (8) 
a a a 
Jedem Punkt des Raumes lassen sich Vektoren bzw. Tensoren zuordnen, die dann Feldvektoren 
oder Feldtensoren heiBen. Erstere lassen sich ausdriicken durch die Form 


a a I at (9) 
Daraus ist ersichtlich, daB die V' bzw. V; als Komponenten in Richtung der kovarianten Basis a, 
bzw. der kontravarianten Basis a‘ gedeutet werden kénnen. 

Die SchnittgréBen (also Normalkrafte) von Schalen werden als die Komponenten in Richtung 
der kovarianten Basis des in Frage kommenden SchnittgréBen-Tensors, bezogen auf die ent- 
sprechende Langeneinheit der Schalenmittelflache, definiert. 

Da die kovarianten Basisvektoren a; im allgemeinen 
keine Einheitsvektoren sind, muf Gleichung (9) durch 
den Betrag von a; dividiert werden. Fiir den auf die 


tee Oy 


yy See 


/ ae o b oN 


a N, 
: “ey i 
Abb. 2. SchnittgréBen in Richtung der kovarianten Abb. 3, Abmessungen der Translationsschale und Schnittkrafte am 
Basisvektoren. Schalenelement (9 = x/L; p =6&/l). 


Langeneinheit bezogenen Wert werden kleine Buchstaben eingefiihrt. Phanomenologisch 1aBt 
sich schreiben 


DP)" (10) 


eas 

Yair 

Hierbei entsprechen die v,, auf Grund der Definition den SchnittgréBen (Abb. 2). 
Im folgenden soll iiber drei wesentliche Punkte berichtet werden: Ableitung und Lésung der 

Differentialgleichung fiir den Membranzustand, Verhdltnisse in den Ecken und das Randstérungs- 

problem. 


2. Membranzustand. Die Schnittkrafte des Membranzustandes werden am zeckmaBigsten tiber 
die von Zerna! angegebene Differentialgleichung ermittelt: 


one « On xo = 
Oy 00 + 30,N “+ pya=90, (11) 


1 W.Zerna, Ing.-Arch. 17 (1949), S. 149. 
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oder ausgeschrieben: 


oN eN?1 oN?” GINPE O01 nr , Or arto ; [2 nro 1 Oe N22 Pes 
+0, 56, + % 36, + % Go, + % Go, + 30, %” + 20, NP 5p Neng cae \ 


Aus den N°” lassen sich die Schnittkrafte wie folgt ermitteln: 


Map) = N°? ee , (12) 


oder ausgeschrieben: 


le / Ayz9 
ay = NOY/ 2, = Nn |/t Na nas) = N® |) ey == INP ih 
Se aa a 
ney = ay), Nas 92) = N® ee —_ sf ne |/ a =e 
Analog gilt fiir den Belastungsvektor 
p = Pi ts + Pe ip -- P3 te D) wobei P; = Va P; . (13) 
AuBerdem ist i 
NU= Ny, Neato Nea. (14) 


Nach Einsetzen der a,-Werte in Gleichung (11) und Multiplikation mit den einzelnen Vektoren 
ergeben sich unter Beriicksichtigung der Bezeichnungsweise von Abb. 3 folgende Gleichgewichts- 
bedingungen: 


No + T +7=9, | 
Yee iy, | (15) 


y’ Naty N, ia Dern | 


Ableitungen nach ? werden durch Striche, solche nach gm durch Punkte angezeigt. 

Es werde vorausgesetzt, da} die beiden Querschnittskurven quadratische Parabeln sind. AuBer- 
dem werde dem Belastungsgesetz eine konstante Schalendicke ¢ (in a;-Richtung) zu Grunde gelegt; 
die sich daraus ergebende Belastung wirke nur in 1,-Richtung. 

Es gilt demnach 


t = konst., Pi=Pp.=0, De +P ; (16) 
Hierbei ist p die Last bezogen auf die Flacheneinheit der Schalenflache. pe) 
Die Gleichgewichtsbedingungen (15) lassen sich auf eine einzige Differentialgleichung fiir N, 

zuriickfiihren; diese lautet unter Beriicksichtigung der eingefiihrten Vereinbarungen (16) 
O +20" 425 p=0, (17) 


wobei 
N, = 8f ® (18) 


ist. Diese Differentialgleichung 1aBt sich geschlossen lésen (vgl. z.B. Pucher!); bei Einfithrung der 
Randbedingungen 


F Nolo= +05 = 0 
un 
N,, go= +05 = 0 
ergibt sich sda 
/D tie 
Sr CVoj x |/ 
© Ife aoa © 
me 910" 51 Lg pale ae BE erin —_——_— cosa? 
¥ can 0,5 Co} a us 0,5 
D Mh, 
na |/ — (|0| — 0,5) nx |/ — (|o| —0,! 
+ & Cre V5 cosnagy + X& CE e a orga (19) 


1 A, Pucher, Beton u. Eisen 33 (1934) S. 298; Bauing. 18 (1937) S. 118. 
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Fiir die Konstanten gilt 


ali SIE) +h 
a= [A(t s) UD eee 


—— 


pee 


o> {sl D 4 Lf 64fD 
er ta : (20) 
vo a Bi +5}; 

2 hee if L cD ct 
Gane 5) 4 Lf ei)" 


Fir die Vorzeichen der mit n behafteten Konstanten gilt 


das vorhandene Vorzeichen ci) oeek hece-cit a ee 


das entgegengesetzte Vorzeichen fiir n= 3,7,... 


Daraus ergeben sich die SchnittgréBen zu 
N, = +8f)/% 6, (21) 


N,=—(8D 0+ sf). 2) 


T=—y [@'dp. (23) 


Damit ist der Membranspannungszustand beschrieben. 


Zur iiberschlaglichen Berechnung interessieren die Schnittkrafte entlang der Rander. Hierfiir gilt 
am Rand ? =~ 0,5: 


fee 2 LWA IDE a 
i tis aru? V a - Nolp = 0 =—gypV Ll? + 16D P> IN 0% (24) 
am Rand g =+ 0,5: 
*. : =— ~~ P+ 6p = 25 
Zs mri \ Nolo = 0 =—gpjil? + lf P> No =. (25) 


Fiir die Schubkrafte lassen sich auf Grund der Lésungen fiir ® keine einfachen Formeln ableiten. 
Zur Berechnung ist es erforderlich, samtliche Integrationskonstanten zu kennen. 

Unter Beriicksichtigung der Ergebnisse des Zahlenbeispiels laBt sich folgende grobe Naherung 
angeben: 


Te Oa wober O22 0,375 (26) 
pts ose 82 
Der Eckwert ergibt sich nach Csonka! mit einem Korrekturfaktor (Erfahrungswert) von 2,0 zu 
0,5 10,51 
tt ie F+D- 
oder 
Ll 


Diese sehr grobe Naherung muB an Hand anderer Zahlenbeispiele noch iiberpriift werden. 
3. Formanderungszustand. Der Verschiebungsvektor )) eines Punktes der Schalenmittelflache sei 
p= ut + vt, + Wz, (28) 
worin Uu, V, Wy seine Komponenten in Richtung der Koordinatenachsen 1,, i, und der Flachen- 
normalen sind. Der Vektor a, ist der Normalenvektor der Mittelflache [siehe (6)]. Fiir die ko- 
varianten MaBzahlen des Verzerrungstensors gilt nach Zerna 
Oa == O° Ds Se Ue Ti ae aa Oy (29) 
Zur Lésung der im folgenden behandelten Fragen (Eckproblem, Randstérungen) ist es erforderlich, 
eine Aussage iiber wy, entlang der Rander zu machen. 


1 P. Csonka, Acta Technica Hung. 10 (1955) S.59; Acta Technica Hung. 11 (1955) S. 231. 
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Nach Zerna’ gilt . 
%,—-Lw 
P= as 025 aru = be ? Wolo = 0,5 a Oey ebay ee * (30) 


Wo 


Fiir uw’ und v’ werde vereinbart, dafB das Randglied eine Verschiebung der Schale am Rand in 
u- baw. v-Richtung verhindert; es gilt demnach 


pe hOs SM p= ws aes vig— £05 =U pewoa— le (31) 


u 


Daraus ergibt sich 


(04 
Woe = 08 ea? Wolo = +405 = —_ - (32) 


Das Elastizitatsgesetz bleibt in seiner allgemeinen Form erhalten?: 
Eton, = Hager? s (33) 


wobei 
1 
Apo, => [4.2 2p + Fx Mpa — M (Ex Bo + Eno Epa) 
mit 
‘a =0 
&9 = — a5 =| 2 &11 — £92 — V > 
oder umgeschrieben 


Etya Oy = (4%)? Ns + 2 dy a. T + [(a42)? —p a] Ne, es 
Et ya Og = Ay Ao Ny + [(Gqg)? + Gy Gog + wa] T + ay age No ’ (34) 
Etja Xe = [(A2)? — wa] Ny + 2 aye Age T + (4p)? ive . 


Damit ist der Formanderungszustand, wie er fiir die Behandlung der in dieser Arbeit behandelten 
Probleme benétigt wird, beschrieben. 


4. Eckproblem. An den vier Ecken der Schale zeigt sich eine eigentiimliche Singularitat, auf 
die Fliigge® bereits hingewiesen hat. Die vorgeschriebenen Randbedingungen kénnen nicht erfillt 
werden. Am Rand ) =-- 0,5 soll N, =0 sein. Aus der Gleichgewichtsbedingung folgt, daB N+ 0 
ist. Betrachten wir die Ecke ? =-+- 0,5, p=+ 0,5, die zu beiden Randern gehért, so ist ersichtlich, 
daB hier ein Widerspruch besteht. 

Die Lésung der Differentialgleichung ergibt in den Ecken bestimmte GréBen fiir die Normal- 
krafte, also nicht Null. Das liegt daran, daB der Teil der Lésung, der aus Fourier-Reihen besteht, 
eben in den Ecken zu Null wird und somit die partikulare Lésung nicht mehr beeinflussen kann. 

Fiir die Schubkrafte ergeben sich in den Ecken unendlich groBe Werte, da die Fourier-Ent- 
wicklung, wie aus Gleichung (23) ersichtlich, wegen der sin-Glieder nicht konvergiert. 

Dies ist auch richtig, da die Schale den in den Eckpunkten auftretenden Konflikt zwischen der 
Differentialgleichung und den eingefiihrten Randbedingungen so lést, daB sie hier die Last mit 
unendlich groBen Schubkraften an unendlich kleinem Kontingenzwinkel tragt; die beiden Druck- 
gewolbe sind infolge N,—N, =0 nicht in der Lage, die értliche Last zu tragen, Ne sie es 
unmittelbar en dem Rande noch tun. 

Bei der Untersuchung der quadratischen, doppeltsymmetrischen Translationsschale ergibt die 
Rechnung unterschiedliche Werte fiir die Normalkrafte in den Ecken, was im Gegensatz zur physi- 
kalischen Tatsache steht. 

Fir p = 1 gilt 


G=—0,25, N,=+8DG=—2D, N, ==(8 PDs 


mit 
a = [4+ .32:L2 D? 
ergibt sich 


L? L2 

a Hip es 4p\-—2D—5, 
also 
Ny + No- 


1 W. Zerna, Osterr. Ing.-Arch. 7 (1953) S. 181. 
2 Green-Zerna, Theoretical Elasticity, Oxford 1954. 


4 W. Fliigge, Das Relaxationsverfahren in der Schalenstatik, Federhofer-Girkmann-Festschrift, Wien 1950, 
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Es wird daher zur Klarung der Verhaltnisse in den Ecken vorgeschlagen, den Bereich, in dem die 


Fourver-Lisung keine richtigen Ergebnisse bringen kann, getrennt zu behandeln und den Form- 
anderungszustand mit heranzuziehen. 


Ate s Durchbiegung wy entlang der Rander gilt unter Beriicksichtigung von (32) und (34) 
und w= 


am Rand ? = + 0,5: 
1 — eo 
E Wolo = + 0,5 = ~ B8LIft [(412)? Ng + 2 ay, Ag, T +- (yo)? Ne ’ 


am Rand g = + 0,5: C) 


E Wop =o ep TDi [(@11)? Ng + 2 Ay, Ayn T + (a49)? Nol : | 


Auch hieraus ist ersichtlich, daB die Schubkraft einen endlichen Wert aufweisen muB, um den 
physikalischen Verhaltnissen Rechnung zu tragen, da fiir T = co auch w) = 0, was nicht még- 
lich ist. 

Zur Lésung des Problems wird wy entlang des Randes (z. B. ®=0,5) graphisch aufgetragen 
und unter Voraussetzung gleicher GesetzmaBigkeit sinnvoll bis zur Ecke y=+0,5 verlangert; 
die sich aus der graphischen Auftragung ergebende Durchbiegung in den Ecken werde mit w 
bezeichnet. 

Als weitere Bedingung, die fiir die Berechnung der SchnittgréBen benétigt wird, gilt die Tat- 
sache, da w in der Ecke gleich gro8 sein muB, unabhangig davon, ob diese als zum Rand # = + 0,5 
oder p=+0,5 gehérig angesehen wird. 

Unter Beriicksichtigung dieser Bedingungen (also T ein endlicher Wert, Ermittlung von w 
aus ‘der graphischen Auftragung an beiden Randern, und wo|, = w|, fiir die gemeinsame Ecke) 
lassen sich die Schnittkrafte in den Ecken eindeutig bestimmen. Es gilt demnach 


Walp eos = (Ny, T,N,)—w, ~~ wolp— 405 = (Nos Ty N,) = 0 (36) 
aus graphischer Auftragung. 
AuBerdem ergibt sich aus der Gleichgewichtsbedingung (15.3) 
Ny =f(N») - (37) 
Man erhalt somit zwei Gleichungen mit den beiden Unbekannten T und Ny, womit das Problem 
gelést ist. 
Fir die Schnittkrafte ergibt sich 
Ng (71 fe — ma ty) = Wt, — Wt, , (38) 
4h (n, ee Ty t,) = W, Ny — We Ny , (39) 
wobei 
wv, = (a2)* gy jpP 8 EWLItf, W, = (143)" gy ph 8 EWLItD, 


(40) 


-, D D 
1 = (ay)? — ve (dg9)? , Ng = (ay)? mee 2 (442)*, 
ty = 2 Ayp Aye , ty = 2 dy Ap - 


5. Randstérungen. Langs der Rander einer Schale lassen sich mit Hilfe der Membrantheorie 
im allgemeinen nicht alle gewiinschten Randbedingungen erfiillen, so daB die lings ihrer Rander 
beanspruchte, im tibrigen aber unbelastete Schale untersucht werden muf. 

Es werde ein Rand ? =, betrachtet. Die Randstérungen, die von diesem Rand ausgehen, 
sind dadurch gekennzeichnet, daB die GréBenordnung aller Funktionen, die diese Stérung beschrei- 
ben, durch Differentiation nach ? zunimmt, durch Differentiation nach py jedoch unverandert 
bleibt. Es kénnen daher naherungsweise alle Funktionen auSer denen mit den héchsten Ableitungen 
nach # vernachlassigt werden. Hinsichtlich eines Randes y =qy gilt Entsprechendes. 

Hierzu ist festzustellen, daB sich alle interessierenden SchnittgréBen (Abb. 4) in Abhangigkeit 
von der Verschiebung w darstellen lassen’. 


1 Siehe FuRnote 2 von Seite 292. 
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Fiir den Rand # = @, gilt 
(1), ; 0. A99 ay Ay 
muy = + Bw ahs / ee 


Or Pll 
N12) = aay 8} w (CE) 


a bo» Ya : 
ae 
Na Baa ay ce 
muy = 9, 


(1),, 
ma = +Bw p, 
Gra 
1112) = — Bw’ cae 


en eye 
Ht a ay 


r11 A99 Va,, 


a a 


(),, 
122) = -- B Ww (1 
(1) 
qa) = —Bw 


° 


(1) Any @ 
=> B we” ee . 
zie e a Va ay 


a; 


Abb, 4, SchnittgréBen der Translationsschale. 


Fir den Rand ~ = q gilt 


2)b 
na) = Eros /_*_, 
Oy F yy Apo 

@).., (4)? 1 
Ni2) = —Bw AES Te : 

Qeng Cr iG Gap 
N22) = -+ Bw” ae 22 - 

= ().. 72 Le 

muy) =—Bw' (l—p) ra 


mei = +B oe mee , 


2), 
m2) = —Bwy-, 


M22) = 0. 


2 
qa) = + Bp 


q(2) =-—B )... ay Vay : 


ee oe 


——S PQs 
Ge) Fy ™m, Ty) 
(1) Hi (21) Ty) 


(41) 
(42) 
(43) 
(44) 
(45) 
(46) 
(47) 


(48) 


(49) 


(50) 
(51) 
(52) 
(53) 
(54) 


(55) 
(56) 
(57) 


(58) 
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Nach Zerna' und Green-Zerna? laBt sich das Problem der Randstérungen auf eine Differential- 
gleichung fiir die Normalverschiebung w zuriickfithren. Fiir den Rand @ = 0, gilt 


w'’ t+Akkw=0. (59) 
Darin ist 
3 (1 — pw?) 
ki = # (a)! ha Das): (60) 


Die Lésung der Differentialgleichung laBt sich wesentlich vereinfachen, wenn fiir die Schalen- 
dicke t folgende Naherung eingefiihrt wird: 


Va s 
t=,—, Ll = t. 
(4y2)" Rites Cy 
Hierdurch erhalt man 
3(1—p?) po 
oh = SO) (2 — konst . (62) 
0 - 


Unter Beriicksichtigung dieser Naiherung gilt unter der Voraussetzung des schnellen Abklingens 


(1) k = 


w= +e"? (C* cosk, 0 + Ct sink, 0). (63) 
Hier sind C* und C* beliebige Funktionen von Q: 
Fiir den Rand g = q gilt analog 


- Va y” 
= (ay)? i vie is (64) 
| Os her 
by a Ve 2 
erie ee 
Co ae (7)? = Kkonst.; (65) 
eo o (Cf cos k, p + Cf sin k, ) . (66) 


Hierbei zahlen die mit einem Querstrich versehenen Koordinaten ? und g vom Rand her zum Null- 
punkt hin. 

Zur Lésung des Randstérungsproblems am Rand #=W, wird die Stelle p=0 betrachtet und 
angenommen, das sich die Verformungen cos- bzw. sin-formig verteilen; dies kann auf Grund 
des Verlaufes von w in erster Naherung an- 
genommen werden (Abb. 5). 


Mit 
Cr=Gfo, G=Gfo (67) 


A = Wy COS P (68) 


und 


gilt 


P 


l | | | ! | ! ! 
-05 -O% -03 -Q2 -Of +0|0 +01 +02 +03 +04 +05 
10, aus Schnittkraften 
———My"=Ny 00s 0p 


Wye =+G, 9” cos kg P+ C, eo” sin ky? . - 
(69) +, [mmly 


Der Aufbau der Formeln fiir die verschiedenen 


; 5 Abb. 5. Verlauf von w entlang des Randes ® = 0) = + 0,5. 
Ableitungen von wy, ist so, da man allgemein 


fiir die SchnittgréBen Y, und die Verformung schreiben kann 
Vx =f vo {+ Ce” [+hysinkyd + kycos ky] + Cye ‘9? [—k, coskyd +hysin ky O]}. (70) 


Die Werte f(,), sowie x, und k, , kénnen fir die einzelnen SchnittgréBen angegeben und in einer 
Ubersicht zusammengestellt werden (siehe Tabelle 1). 
Zur Ermittlung der Konstanten C; miissen die Klaffungen (d. i. die Differenz aus den Verfor- 
mungen der Schale am Rande und denjenigen des entsprechenden Randtragers) zu Null werden. 
Die sich aus den Einheitsbelastungsfallen C,=1 ergebenden Schnittkrafte der Schale werden 
als Belastung des Randtragers aufgefaBt und die sich hieraus ergebenden Verformungen ermittelt. 


1 Siehe FuBnote | von Seite 292. 
2 Siehe FuBnote 2 von Seite 292. 
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Wy Up 4 yp 


ena re 


et 8S 


io 20 


ed el 


v| 
resi; 


BS il dp 


T 


T U 


T 


Sin 


IT 


417A Dente : 
Eas Geng 
17 aes ce 
LET eae a 
07 ane | 
41 TI — 
ial tal(r— D9 


[ei =D ef] 2 q 


(<r'— 1) 9 ' 
(i—Dehag ' 


°@ = G puny 1ystsseg”~ “T T]9qR J, 


d 22 $00 + 


du soo + 


du s09 + 


uw soo + 


db 2 s09 + 


dh 1 809 . 2 — 


b 2 UIs . 2 — 


(Op 


(Z22)ay 


(ZDay 
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Als rechte Seiten des sich aus diesen Uberlegungen ergebenden Gleichungssystems miissen die 
Verformungen der Schale unter Eigengewicht (Membranzustand) und die Verformungen des 
Randtragers infolge der Membrankrifte der Schale und seines Kigengewichtes beriicksichtigt werden. 

Um schnell zu einem anschaulichen Ergebnis zu gelangen, erscheint es angebracht, zwei Sonder- 
falle zu untersuchen. Hierfiir werden folgende Annahmen gemacht: 


a) Der Randtrager sei ideal starr, d. h. 


w=0, —=0 


b) Fir den Membranzustand sei @w/@ der Schale an der Stelle # =-+0,5, y =0 vernachlassig- 
bar klein. 


Fir den Sonderfall der Einspannung ergibt sich folgende Matrix: 


o— —_ Membran- 
Swe Gs ee 
Schale + 1,0 0 Obie 
oF Trager 0 0 Oo 
A + 1,0 0 | —wiy, 
Ow Schale Leese kg —ky» 0 
— Trager 0 | 0 | 0 
oo ee ee ct | eo eS 8 eS > 
A + keg —kg 0 


* Uber die Bedeutung von wag vgl. Gleichung (80) und Bild 5. 


Ks ist sofort zu sehen, daB 


GC=—wy, und C,=+G (71) 
sind. Damit vereinfacht sich Gleichung (70) zu 
YE = f(g) x Cpe? [+ (ky + hy) sin ky F  (Itg — hy) -c08 hy 8] - (72) 
Fir das Quermoment gilt 
ME= — wy, 2) k,e'8° (sin ky & — cos ky 9) - (73) 
An der Stelle ? = + 0,5, mp = 0 ergibt sich 
- E®y3 Lin” 1 9 EetLf 
,=— ee 31 ——___ *_— : TA 
. 3 6 Va ayy ay: 2 YL? + 16 D? me 


Damit kann abgeschatzt werden, bei welcher GréBenordnung von wy das Quermoment so grof 
wird, daB es bei der Bemessung beriicksichtigt werden mub. 


Der andere GréBtwert des Quermomentes liegt bei J” = m/2 ky und betragt 


7 


1 
a — <M. (75) 


ME = — wy %y ky e— 7 sin 


Fiir den Sonderfall Gelenk ergibt sich nach Auflésung der Matrix 
Yo = fi xy Cre *” (hy sin ky D + hy cos ky 8) , (76) 


M§ = —wy%,hkye “0” sin ky @. (77) 
‘ : 1 : aS I 
Das gréBte Quermoment betragt hierbei etwa —-- M}, und liegt bei 0” = | ese 


Es darf noch erwahnt werden, daB die Dehnungen ¢ am Schalenanschnittspunkt bei Schale 
und Randtrager iibereinstimmen miissen. Bei Untersuchung der angefithrten Sonderfalle ergibt 


sich folgendes: 


| | 


} : | Membran- 
| C | C, | zustand 
w (Schale) + 1,0 0 | —iwyr 
N22) = No=é -t (Schale) %9(-+ ky) —k,=0 — Now 


20 
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Es zeigt sich, daB C, beide Gleichungen erfiillt [vgl. (32) und (34)]. 
Aus der ersten Gleichung folgt 


boo (Ge2)” a7 Ja ee 1 a 21/8 
ae ben Et fa MLE Et Lr (722) Opp PME 


oe A59 Vay, aa ae 1 Nee 
pena eH, n” PM #9 (+ ko) %M 
Damit ist die zweite Gleichung erfiillt. 
Die endgiiltigen Schnittkrafte sind dann 


Map) = Map) + MG p) > (78) 
Fhe, 1B 
nay = Ny + niin. 


6. Zahlenbeispiel. Die Berechnung der Schnittkrafte des Membranzustandes bereitet keine 
Schwierigkeit; als Beispiel wird der Verlauf von N, am Rand ) =) angegeben (Abb. Oh Ab- 
messungen: L = 27,50m; 1 = 21,26m; D=3,34m; f=—1,95m; t= 6cm; p = 0,35 t/m?. 


-05 -O+ -03 -O2 -01 +00 +01 +02 +03 +04 +*05——? 
T T T T 


ee 
7) 
i i ! i ¥ 

“05 “03 -Q2 -O01 +010 +01 +02 +03 +04 Fp) Pie 

Nepit/m] Poel E=35-10 thn 
Abb. 6. Verlauf von NV entlang des Randes Abb. 7. Verlauf von w entlang beider Rander (sinnvolle Verlangerung 

pele ears bis in die Ecken). 
9 = ee 0 


Der Verlauf der nach Gleichung (32) unter Benutzung von Gleichung (34) berechneten Ver- 
schiebung wy ist aus der graphischen Darstellung Abb. 7 ersichtlich. Hieraus ergab sich der Wert 
fiir w in den Ecken, welcher der Berechnung der SchnittgréBen in den Ecken zugrunde gelegt wurde. 


Es ist 
5: (79) 


Fir die Berechnung der Randstérungen gilt (vgl. auch Abb. 5) 


* 
Wy = Wo |o= + 0,5» (80) 
P > 0 
wobei 
hg = 
wy = Wy—w 
1st. 


Den Berechnungen fiir die Randstérungen wurde die Ubersicht Tabelle 1 zugrunde gelegt und 
die beiden angefiihrten Sonderfalle behandelt; der besonders interessierende Verlauf des Quer- 
momentes ist in Abb. 8 dargestellt. In der Tabelle 2 sind die Ergebnisse zusammengefaBt. 

Aus Abb. 9 ist ersichtlich, daB die eingefiihrte Veranderlichkeit die Schalendicke t in vertret- 
baren Grenzen gegeniiber der konstanten Dicke der Membrantheorie variieren laBt. AuGerdem 
zeigen die Ergebnisse, da® durch die Lésung diejenige Randbedingung nicht beeinfluBt wird, 
welche durch die Membrantheorie erfiillt worden ist, also mui) = Ng = 90 entlang des Randes. 
Uber ni2) = N, macht die Membrantheorie fiir diesen Rand keine einschrankende Aussage, so 
daB die Randstérung diese Gré&e beeinflussen kann, was auch tatsachlich der Fall ist. 
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Wie man leicht zeigen kann, betrigt das GréBenordnungsverhiiltnis (Fa) w fir das Zahlen- 
beispiel, fiir Kinspannung etwa = = 17,5 und fiir Gelenk etwa Boe es 27,7, wahrend das Ver- 


haltnis (@w/0y~)/w in beiden Fallen 2/1 = 3,14 ausmacht. 


5,7 


~ 60,0 


~ 300 


a 


+0,025 +6,050 +0075 +0700 +0925 ia 0150 +0175 +0200 
M2) [mkg f m 


S 


+ 600 +0025 +0050 O05 +000 +0105 +000 +01 +0200~ 
+7172) [mkg/ m| 


Abb. 8. Verlauf des Quermomentes bei Einspannung in den Randtrager (oben) und gelenkigem AnschluB (unten). 


Abb. 9. Veranderliche Schalendicke t, aufgetragen iiber der Grundriffliche der Schale. 


20* 
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Tabelle 2 
a EEE ee 
chnittgréB Punk Randstérungen Summe 
es i= as beat Einsp. Gelenk Einsp. Gelenk 
WeaiNe Mitte » = 0 0 Eeaiwe: 0 Le=e013 0 
ae Ecke = 0,5 — 6,30 0 0. | = 6300 =noso 
UP Mitte 0 0 0 0 0 
ie Ecke — 27,3 4 1,65:60 0,83 | 2 25005 eee 
ee Mitte —11,27 == 10,56 | = 10:56) |--"21,83 1) sess 
t[m Ecke 74 0 0 | Sena | ae 
wy) | Ewk baw. Ew | Mitte +1,73 | —5670 | +5670 | + 5670 0 0 
mm. | t/m Ecke +>3,39 0 0 0 0 0 
MeO al: Mitte 0 40,183 | 0  |emonigay or 0 
mt/m 
S ‘ y a 7 | 
Hoos Mitte 0 +0,35 | +0,17 | +035 | +0,17 
m7 
Ws 
Po - x es 2 
| Se 
116 =a so IE 
| | TT Witelnann P=05 
= ; | 
110 a Z 
405 


Abb, 10. Verschiedene Belastungsfunktionen entlang des Randes 0 = 0. = + 0,5. 


Dieses die Voraussetzung befriedigend erfiillende Ergebnis ist abhangig von dem Wert k,, 
der im betrachtenden Beispiel 27,7 betragt, also sehr groB ist und demnach auch die Randstérungen 
sehr schnell abklingen laBt, was ebenfalls angenommen worden war. 


Der Vergleich mit anderen Berechnungsverfahren (z. B. Fliigge! und Csonka? in bezug auf die 


geometrische Form, auf die Belastung und die sich daraus ergebenden Schnittkrafte ist sehr auf- 
schluBreich. Fiir die Belastungsverteilung gilt z. B. nach Fliigge 


Da 1 


Po Y1—sin?Osin2p 
(# und @ sind hier Offnungswinkel), nach Csonka 


= i 
Po 
nach Mittelmann 
Pa! Va 
Po LI° 


1 W. Fliigge, Statik und Dynamik der Schalen, S. 99, Berlin 1934. 
2 Siehe FuSnote 2 von Seite 291. 
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Die graphische Auftragung dieser Funktionen entlang des Randes ? = 0,5 zeigt Abb. 10. Der 
Verlauf der Schnittkraft Ny ist fiir die Schnitte ? = 0 und ? = 0,5 in Abb. 11 dargestellt. 

Eine Untersuchung des Formanderungszustandes und des Einflusses von Randstérungen ist 
nicht méglich, da deren Beschreibung hier erstmals praktisch durchgefiihrt wurde. 

Eine ausfiihrliche Gegeniiberstellung der bisher bekannten Berechnungsverfahren (u. a. 
Pucher', Tungl?, Salvadori*) diirfte weitere aufschluBreiche Ergebnisse zeitigen. 


| a 47100 a 


P 


hae ee ESL 
“05 -04 -03 -02 -01 +40)0 +07 +02 +03 +04% +05 


+70 
N. Verlauf nach 
f Fligge 
ST 
—— — Mifelmann 


Abb. 11. Verlauf von N,, entlang des Randes } = #, = + 0,5, — Vergleich zwischen verschiedenen Berechnungsverfahren. 


7. SchluBbetrachtung. Die bisher bekannten Berechnungsmethoden fiir Translationsschalen 
befriedigen nicht in allen Teilen; einmal wegen ihres erheblichen Rechenaufwandes, zum anderen 
wegen ihrer Unvollstandigkeit in bezug auf den Formanderungszustand und das Randstérungs- 
problem. Mit der vorliegenden Arbeit ist versucht worden, einen neuen Weg zur Berechnung 
dieser Schalenform aufzuzeigen, dessen Gleichungen und Ergebnisse sowohl alle Gleichgewichts- 
als auch die Formanderungsbedingungen erfillen. Dies konnte durch die Einfiihrung verschiedener 
Naherungen erreicht werden, deren berechtigte Verwendung nachgewiesen wurde. Ein Vergleich 
mit den anderen Verfahren ergibt wissenswerte Aufschliisse und zeigt die ZweckmaBigkeit der an- 
gewandten Methode. AuBerdem ist es nunmehr méglich, die GréBe der durch das Randstérungs- 
problem auftretenden SchnittgréBen zu berechnen und ihren Auswirkungsbereich zu ermitteln, 
so daB auch Ausfiihrungen mit schwierigen Randgliedausbildungen damit der Berechnung zugang- 
lich werden. 


(Eingegangen am 14. November 1957.) 
Anschrift des Verfassers: Dip].-Ing. Goswin Mittelmann, Frankfurt (Main), ZiegenhainerstraBe 89. 


1 Siehe FuBnote 1 von Seite 290. 7 
2 E. Tungl, Osterr. Ing. Arch. 10 (1956) S. 308; Osterr. Bauz. 11 (1956) S. 274, 
3 M.G. Salvadori, ACI Proceedings 52 (1956) S. 1099. 


302 W. Kellenberger: Biegeschwingungen einer unrunden, rotierenden Welle Ingenieur-Archiv 


Biegeschwingungen einer unrunden, rotierenden Welle in horizontaler Lage 
Von W. Kellenberger 


1. Einleitung. Der Querschnitt einer unrunden Welle besitzt zwei senkrecht aufeinander 
stehende Haupttragheitsachsen durch den Schwerpunkt, denen zwei verschiedene Biegesteifig- 
keiten zugeordnet sind. Die durch eine sogenannte Unwucht erzeugten Schwingungen werden 
dadurch abhangig von der Lage der Unwucht beziiglich des Querschnittes. Ist die Welle hori- 
zontal gelagert, so werden zusatzlich doppelfrequente Biegeschwingungen angeregt, die bei gut 
ausgewuchteten Wellen durchaus die Gré®enordnung der durch Unwucht erzeugten Schwin- 
gungen erreichen kénnen. Sie sind im Gegensatz zu diesen durch mechanisches Auswuchten nicht 
beeinfluBbar. 

Das Problem gewinnt praktisch immer mehr an Bedeutung durch die enorme Leistungssteige- 
rung an Turbogeneratoren, deren zwei verschiedene Biegesteifigkeiten bedingt sind durch die 
iiber den Umfang der ,,Welle“‘ 
ungleichmaBig verteilten Nuten 
fir die Stromleiter. Die mecha- 
nische Festigkeit solcher schnell- 
laufender Rotoren setzen einer 
nennenswerten VergréBerung ihres 
Durchmessers enge Grenzen, wo- 
durch sie sehr schlank werden. 
Man ,,fahrt** heute schon zwischen 
der zweiten und dritten normalen 
kritischen Drehzahl bei Nenn- 
betrieb. 

Uber das entsprechende Pro- 
blem mit nur einem Freiheits- 
grad, der masselosen Welle ver- 
schiedener Biegesteifigkeit mit 
in der Mitte aufgesetzter Scheibe, 
wurde verschiedentlich und mehr 
oder weniger ausfiihrlich geschrie- 
ben. Im Folgenden soll das Pro- 
blem im Sinne der Kontinuums- 
mechanik an der gleichmabig mit 

Masse belegten Welle verschiede- 
Abb. 1. Orientierung der Welle und Krafte im ruhenden System. (2) = Auslenkung ner tiber die Lange jedoch kon- 
des Wellenmittelpunktes W; “e(2) = Exzentrizitiit des Schwerpunktes S. stanter Biegesteifigkeit behandelt 
werden. Damit ergibt sich dann 

auch der Kinflu® der ungleichmafigen Biegesteifigkeit auf die Oberschwingungen. 


2. Differentialgleichungen. In Abb. 1 ist ein Querschnitt durch die Welle dargestellt. Dabei 
ist x, y, 2 ein ortsfestes, u, v, z das mit der Welle rotierende System, wobei u und v die Richtungen 
der Hauptschwerachsen des Querschnittes angeben. Damit kann sich der Querschnitt nur trans- 
latorisch beziiglich u, v bewegen. Die y-Achse habe die Richtung des Eigengewichtes der Welle 
und die z-Achse die horizontale Lage der unverformten Wellenachse. Der Punkt 0 ist der Durch- 
stoBpunkt der unverbogenen Achse, der Punkt W ist der geometrische Mittelpunkt der Welle 
und S der Schwerpunkt, somit e = WS die > Exzentrizitat‘. 

Bezeichnet man mit E J, = B, baw. E J, = B, die zwei verschiedenen Biegesteifigkeiten 
(Federzahlen) der Welle in der (u, z)- bzw. (v, z)-Ebene, und bleibt man im Rahmen der Balken- 


1 A. Stodola, Dampf- und Gasturbinen, S. 931, 5. Aufl., Berlin 1922. D. M. Smith, Proc. Roy. Society A 142 
(1933) s. 92. H.D. Taylor, Journ. Appl. Mech. (1940) A-71. W. R. Foote, H. Poritsky u. J. J. Slade, Jr., Journ. 
Appl. Mech. (1943) A 77. C. B. Biezeno u. R. Grammel, Technische Dynamik (Bd. II, S. 64, Berlin, Gottingen, 
Heidelberg 1953. O. Kraemer, Neuere Feststellungen iiber Biegeschwingungen an Kurbelwellen, Miinchen 1956. 
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biegetheorie von Bernoulli, so schreibt sich der elastische Widerstand je Langeneinheit Fund 
F, in Richtung u und v 
Fi=—B,ulv 


; BeBe. (1) 
F = — Bp ov ‘ : 


Rémische Exponenten bezeichnen im folgenden Ableitungen nach der Stabachse z Punkte nach 
der Zeit t. Die allgemeinen Transformationsformeln hei®en 


x =ucoswt—vsinwt, v y 


y =usinwt+vecoswt, 


(2) 


u=xcoswt+ysinwt, 
v =—xsinwt+ycoswt. 


Es ist zweckmabig, die Bewegungs- 
differentialgleichungen in einem mit der 
konstanten Winkelgeschwindigkeit w dre- 
henden K oordinatensystem zu formulieren. 
Dabei sind neben den absoluten Kraften 
die fiktiven Zusatzkrafte der Relativbewe- 


gung, die Fliehkraft Z und die Corioliskraft 
C zu beriicksichtigen. Die Aufspaltung in 


Komponenten der an der Welle angreifen- (2) 
den Krafte (Abb. 2) fiihrt dann zu folgen- Abb. 2. Krafte und Scheinkrafte im rotierenden System. —r(z) = Auslenkung 
der Tabelle!: des Wellenmittelpunktes W; e(z) = Exzentrizitat des Schwerpunktes S. 


Krifte je Langeneinheit parca in patra: eter 
Eigengewicht G = mg gmsinwt | gmcoswt 
Fuhrungskraft F(F., Fy) — B,ulV — B, viv By, < By 
(Federkraft) 
a ss 
Tragheitskraft T = — mr, — miis —mis 
Flichkraft Z = m o?7, 2 2 el ara i 
iehkraft Z = mw? fF, m W* Us m @ vs vs =v +esinf£ 
oe 
Corioliskraft C = 2 m [vs ret ©] 2movz —2mou, 


Dabei ist m die Masse der Welle j a Langeneinheit, g die Erdbeschleunigung und f(z, t) der Lage- 
winkel der Exzentrizitat WS = ¢ gegen die u-Achse. Das dynamische Gleichgewicht Pe Krafte 


C + F + T + 7 — C = 0 liefert in Komponenten des rotierenden Systems ausgedrickt die 
Gleichungen 


IV 


: ; B ; 
ti —2wv +—uV —o2u=ew* cosh + gsinwt, 
m 


(3) 


is 5 B : 
v+2wu+—vV — o?v =ew*sinf + gcoswt. 
™m 


Die zweite dynamische Gleichgewichtsbedingung fordert das Verschwinden der Drehmo- 
mente aller beteiligten Krafte. Bedeutet 9 = mk? das Tragheitsmoment je Langeneinheit, k 
also den Tragheitsradius und M das aufere Drehmoment der Welle, so gilt 


OB = M—F esinf + F,e cos B 
oder mit Riicksicht auf obige Tabelle 
ee M | e 1 (B, bi sin B — B, v'¥ cos f) . (4) 


mk? hk? em 


Die Gleichungen (3) und (4) stellen ein System von drei linearen, gekoppelten, partiellen 
Differentialgleichungen sechster Ordnung mit konstanten Koeffizienten in u, v und f dar. Das 
entsprechende System in ortsfesten Koordinaten wiirde zeitlich harmonische Koeffizienten auf- 


1 Schubkrafte, rotatorische Tragheit, Kreiselwirkung und Dampfungskrafte sind vernachlassigt. 
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weisen, aihnlich der Mathieuschen Gleichung im Gebiete der gewéhnlichen Differentialglei- 
chungen. Der Vorteil einer Transformation liegt in der Vereinfachung der Differentialgleichungen 
durch die Vermeidung harmonischer Koeffizienten. 

Wir beschranken uns noch auf den stationaren Betricb M = 0. Bleiben wir auBerdem im 
Bereich der GroBmaschinen, so ist das Verhaltnis Exzentrizitat zu Tragheitsradius von der 
GréBenordnung < 10-4, also der Ausdruck e?/k? < 10-8 eine sehr kleine Zahl. Man darf deshalb 
den zweiten Summanden auf der rechten Seite der Gleichung (4) in erster Naherung vernach- 


lassigen’. Gleichung (4) fiihrt damit auf p= 0 oder 6 = at + b, wobei a und 6 zeitlich kon- 
stante GréBen bedeuten. Im stationiren Betrieb besteht somit in erster Naherung auch ein 
stationirer Drehzustand. Setzt man noch voraus, daB w die mittlere, konstante Winkelge- 
schwindigkeit der Welle bedeutet, so verschwindet offenbar a, und § (= b) wird eine zeitlich 
konstante GréBe. Die Bewegung der rotierenden Welle wird, und diese ist das wichtigst Ergebnis 
der Naherung, allein durch die Differentialgleichung (3) beschrieben, in der jetzt B eine zeitlich 
kontante, im allgemeinen noch von z abhangige GréBe darstellt. 

Die beiden Gleichungen (3) mit denen wir unsere Untersuchung allein weiterfiihren wollen, 
sind immer noch gekoppelt. Zusammen mit den noch zu formulierenden Lager- und Anfangshe- 
dingungen wird ein Randwertproblem achter Ordnung in rotierenden Koordinaten definiert. 


3. Homogene Gleichungen (Freie Schwingungen). Der vereinfachte Produktansatz 


u=Z(z) Tt), v= 4) T(t) (5) 
fihrt auf folgende Separation der Variabeln: 
eds fis TAS m 7 Ts ype 
@ (t 20 at) 4 Z 0 a -2w T, ie eo == 0); 
was jederzeit und an jedem Ort nur méglich ist unter den drei Bedingungen 
giv 
ae hae (6a) 
Te z 
u v 2 Se ya 
T, 20 a k ae (6b) 
T, Ts 9 By, 
7 oes ed ae (60) 


wo k eine noch zu bestimmende Konstante sein soll. 
Die Gleichung (6a) fiihrt auf die bekannte Lésung 
Z=Asinkz+ Beoskz + Csinhkz + Dcoshkz (2) 
mit den vier Integrationskonstanten A, B, C, D, welche aus den Rand- oder Lagerbedingungen 
zu bestimmen sind. 


Im Falle der an den Enden starr, aber frei-drehbar gelagerten Welle lauten die Lagerbe- 
dingungen? 


Z0)=0,-Z20)=0, ~“Z2'0)=0, ZO) =0: (8) 


Sie ergeben zusammen mit (7) vier homogene Gleichungen mit der Vertraglichkeitsbedingung 


na 
k=— (n= 15.273 eee) (9) 


“a 


und der einzigen, von Null verschiedenen Konstanten A, wodurch 
NMI 
—s% 


Z = A, sin 1 (1 == 1525308 8) (10) 


1 Greift man der Untersuchung etwas vor und beniitzt die Gleich a i 
in der etwas tibersichtlicheren Form memibmeneins aac r ier or ee 
[ib ces Bom Be 


on e2 @2 e ; 7B 
Dy == = inp — = 
Bais |g )uame— (+ g)ooe Se KELE 

darstellen, in der u, v und e Langen von eleicher GroR d 2 /}p2 Assi ; 

sehen werden kénnen, so dab f e O20. zt, 0 wird. enordnung, ¢*/h* dagegen als vernachlassigbar klein ange- 


* Andere Lagerbedingungen lassen sich ganz analog behandeln, bringen aber keine neuen Effekte. 
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wird. Mit Riicksicht auf die Superponierbarkeit der Lésungen kommt 


Z Seu 8 sin“* 2. (11) 


i 


Die Gleichungen (6b) und (6c) lassen sich mit den Abkiirzungen! 


und mit (9) auf die Form bringen 


Tae orto (a: a) T= 0., T, +20 T, + (w?—w) T,=0- (13) 


Der Schwingungsansatz 
T,, = Aele'= A (cosot+jsingt), T, = Beie* = B(cosgt + sing t), (14) 


wobei A, B und ¢ im allgemeinen komplexe, von t unabhangige Werte annehmen kénnen, fiihrt 
auf das homogene Gleichungssystem 


4 | : | 
= (15) 
(@?2 — w?) — 0? —j2we | 0 
+j2o@ (w2—ow?)— oe? | 0 
mit der Vertraglichkeitsbedingung 
Ae 2 
Det = 04 — 2 0? (* Loe x sd (a? — w?) (w* — w?) = 0, (16) 


die zugleich die gesuchte Frequenzgleichung im rotierenden System darstellt. Sie 
liefert vier und nur vier verschiedene Eigenwerte 0, welche entsprechend dem Ansatz (14) als 
Kreisfrequenzen in Betracht SLES Aufgelést nach o? wird (16) 


wz - we 


gh at + EEF W202 (wt + 08) + 4 (808). (17) 

a) Diskussion der Heat tor tb Die Diskussion der Frequenzgleichungen 
im rotierenden System ist relativ einfach, da sie lediglich biquadratisch in @ ist. 

Das Kriterium fiir die Stabilitat der Lésungen folgt unmittelbar aus (14) mit der charakteri- 
stischen GréBe o. 

Zunachst erkennt man aus (16), daB sich fiir jedes ~ keine komplexen, sondern nur reelle, 
positive oder negative Werte fiir 07, somit nur reelle oder imaginare Werte fiir 9 ergeben. Ist 0 
reell, so fiihrt der Ansatz (14) auf periodische Schwingungen mit der Kreisfrequenz @ im rotie- 
renden System, die Lésung ist stabil. 

Ist o hingegen imaginar, so ergeben sich Kreisfunktionen mit komplexem Argument j @ ¢, 
die durch Hyperbelfunktionen mit reellem Argument 0 t zu ersetzen sind. Da letztere mit der 
Zeit unbeschrankt zunehmen, ist die Lésung instabil. 

Fir 0 = 0 geniigt der Ansatz (14) nicht mehr. Es mu8 ein weiteres partikulares Integral ge- 
funden werden, das dann zu instabilen Lésungen fihrt. 

Die in 0? quadratische Gleichung (16) betiece zwei reelle positive Wurzeln o? fiir 


[w +- Ou 2 ae aes > (w? — w?) (ow? — ww?) > 0. 


Es ergeben sich dann zweimal zwei s t’abile Lésungen fiir 
(i), 0) = bYe 
Weil in (1) B, > B,, vorausgesetzt wurde, ist nach (12) w, > w,,. 
Andererseits hesitzt (16) eine reelle positive und eine reelle negative Wurzel 0? fiir 


(wo? — w?) (w? — wz) < 0 
1@,, und my muBten mit einem weiteren Index n versehen sein. Es soll einfachheitshalber davon abgesehen 
werden, wenn keine Verwechsiung moglich ist. 
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und somit zwei stabile und zwei instabile Lésungen fir 
O, << W<, 


womit dieser ganze Bereich instabil wird. 

Fir w = +, undw = +,, d.h. fiir die Grenzen zwischen stabilem und instabilem Ge- 
biet wird 9 = 0, der Ansatz (14) liefert eine (zeitunabhaingige) Konstante im rotierenden System. 
Er geniigt jedoch nicht und ist durch ein zweites partikulares Integral zu erganzen, das mit der 
Zeit unbeschriankt anwachst. Die Grenzen sind deshalb ebenfalls zum instabilen Gebiet zu 
zihlen. Damit sind auch samtliche Méglichkeiten ausgeschépft. 

Der einzige instabile Bereich liegt vor fiir Winkelgeschwindigkeiten w, =< w= w,. Da je- 
doch den Gréen w, und w, nach (12) wegen n = 1, 2, 3,... verschiedene Werte zukommen, 
wiederholen sich die instabilen Bereiche entsprechend n. Die physikalische Bedeutung von w, 
und w, geht aus (12) hervor. Es sind die Kreisfrequenzen der Welle fiir ebene Biegeschwingungen 
in Richtung der Hauptachsen des Querschnittes, wobei n die Ordnung der Oberschwingung 
angibt. 

Die Frequenzen der méglichen Eigenschwingungen im ortsfesten System folgen aus 
der Beziehung 


yv=ot+e, (18) 


Da es sich ja um eine Drehung des bewegten Systems mit der konstanten Winkelgeschwindig- 
keit w handelt, und zwar ist vy die Kreisfrequenz im ortsfesten, 0 diejenige im rotierenden System 
nach Gleichung (17). 

Fiir die graphische Darstellung des Zusammenhanges vy = »(~) empfiehlt es sich, von folgen- 
den drei definierten GréBen Gebrauch zu machen. Es sei 


By + Bu 


B=—5 (19a) 
der arithmetische Mittelwert der Biegesteifigkeiten, 
SRB B, See Bea Be 
a Ee ean a ele 
die bezogene Unrundheit der Welle, und 
ee ee ial Mgr ae 19 
ae? | M=L23,.., (196) 


die kritische Winkelgeschwindigkeit einer runden Welle der Biegesteifigkeit B auf starr-frei- 


drehbaren Lagern. n bedeutet die Ordnung der Oberschwingung. 


g 
Daraus leiten sich die Beziehungen 


= é = & 
B, = B(1—+), B, = B(1 +4) (20) 
ab und mit (12) auch 
le € w? — w2 
of = B(1—4), ale ae pa 
of =B(1 +5] 50 wit wt ee 


Die Frequenzgleichung (18) im ortsfesten System lautet mit diesen neuen Werten in 
dimensionsloser Form alsdann 


ye Ons ae ea, sae 
oO | Op Oo 1 +(3) +2 (By +e - (22) 
Sie ist in Abb. 3 fiir die Parameter € = 0,0.2, 0.5, 1.0 und 1.5 graphisch dargestellt. Es lassen sich 
unmittelbar folgende Aussagen machen: 
1. Es ist zu beachten, daB @ nach (19c) mit einem Index n zu versehen ware, der die Ordnung 
der Oberschwingung angibt. Das Diagramm wiederholt sich somit fiir alle Oberschwingungen. 
2. Die Kreisfrequenzen » im ortsfesten System sind drehzahlabhangig und damit auch die 
EKigenschwingungen. 
; 3. Hiner festen Winkelgeschwindigkeit w der Welle und einer gegebenen Unrundheit ¢ sind 
im Ganzen vier Kigenkreisfrequenzen y bezogen auf @ zugeordnet; somit lassen sich vier Typen 
von Eigenschwingungen, jede mit Oberschwingungen unterscheiden. 
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A. Die Bewegung ist im allgemeinen nicht mehr periodisch, da die Kreisfrequenzen y in 
keinem rationalen Verhaltnis zueinander stehen. 


3. Die Strahlen vom Nullpunkt aus deuten ein konstantes Verhaltnis der Kigenkreisfrequenz y 
zur Winkelgeschwindigkeit w an. Man erkennt, daB z. B. Eigenschwingungen mit Drehzahlfre- 


quenz nur méglich sind fiir @ = ay/1 —+= =o, und w=@y/1 + > =w, wahrend Eigen- 


She 
x 


bezogene Eigenkreistrequenz 


§ 
iS 
Sile 


SS 
-) 
Sy 

a) 


Wr Ae 
: bezogene Winkelgeschwindigkeit der Welle —— 


Abb. 3. Eigenkreisfrequenzen » im ruhenden x—y-System in Funktion der Wellen-Drehzabl. 


05 


Ss 
S 


Unrundheit der Welle 


48 


10 


Abb. 4, Resonanzlagen und instabile Bereiche. — Kurve la und lb = es 
w 
Resonanzlagen fiir Schwingung mit doppelter Drehzahlfrequenz (Schwere). 


Resonanzlagen fiir Schwingung mit Drehzahlfrequenz 
v=o 


(Unwucht); Kurve 2 = ee 
wo ]v=0;20 
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schwingungen mit doppelter Frequenz der Drehzahl » = 2 w, sog. doppelfrequente Kigenschwin- 
gungen, nur fiir eine einzige, bestimmte Winkelgeschwindigkeit w < @/2 méglich sind. 


b) Instabiler Bereich. Zwischen m, und w, liegt der in Abb. 4 schraffierte, instabile Be- 
reich, dessen Breite mit zunehmender Unrundheit ¢ der Welle zunimmt. Mach (21) ist 


o,—o, =a(Yi +4 —fi—=), (23) 


was in Abb. 4 graphisch zum Ausdruck kommt. 
Fiir geringen Unterschied der Biegesteifigkeiten, wie etwa bei Turbogeneratoren, ist nach 
(19b) auch die Unrundheit der Welle klein. Man schreibt dann fiir € zweckmaBiger ¢€ 
e>e<l, 


so daB in erster Naherung ¢? gegen 1 vernachlassigt werden darf. Die Breite des instabilen Be- 
reiches nach (23) wird dann in erster Naherung 


oa 
On iO pig 
und ist somit unbedeutend. Fiir die Oberschwingungen wachst der Bereich allerdings nach 
(19¢) mit n?. 

Die den Grenzkurven selbst entsprechenden Werte fiihren ebenfalls zu instabilen Lésungen, 
was noch gezeigt wird. 


c) Stabile Lésungen. Die zwei reellen, positiven Wurzeln der Frequenzgleichung (16) fir 
0? seien mit 9? und 03 bezeichnet, so daf sich vier Kigenwerte +e, und + @, fir @ ergeben. Der 
Ansatz (14) fihrt dann auf vier Kigenfunktionen mit je vier Konstanten fiir T,, und T,. Die zu 
T, und T,, gehérigen Konstanten sind aber durch die Beziehung (15) verknipft, so daB nur 


noch vier Integrationskonstanten — entsprechend dem vierten Grade der Differentialgleichung 
(13) — C, bis C, bleiben. Die Rechnung ergibt 
T,, = C, cos 0, t + C, sing, t + C, cos o,t + C, sinost, (24) 
T, = 0, C, sin 0, t — 6, C, cos 0, t + 6, C, sin 0, t — 0, Cy cos Ont, 
mit 
_ @=(@3 =o) _ 200; ey 
emery a ceety Cs 


Nun interessieren vor allem die Koordinaten x, y im ortsfesten System. Die Transformation 
nach (2) liefert vorerst 


A/a (oe oe 


x 


y = Zz) Tit), T, = T,sinwt + T,coswt 


(25) 


und mit (24) nach einigen Umformungen 

T,, = D, cos (w — 9,) t — Dy sin (w — @,) t + Dg cos (w — Q) t — Dy sin (w — 0.) t 

+ 0, D, cos (w + 0) t + 0, Dy sin (w + 0,) t + 8, D3 cos (w +0) t + J, D, sin (w + 05) t, 

T,, = D, sin (w — 9;) t + D, cos (w — e,) t + Dj sin (w — 9.) t + D, cos (w — ,) t 

+ 8, D, sin (m + Q;) t — 8, Dy cos (w + @,) t + B, Dy sin (w + 02) t — 0, D, cos (w + 05), 
. — % —@,— o’)—200; _ 20; —o? + @2—o%) 


iB 0? —— (2 — «w?) + PO) 0; =>  Q; + 0? a= (2 — w?) (u == Le 2) > (26) 


wobei D, bis D, neue, beliebige Konstanten bedeuten. 


Der Gleichung (26) kann mit Riicksicht auf die Kreisfunktionen mit gleichem Argument auch 
die Form gegeben werden 


T, = H, cos (vt + 9) + Hy cos (vt + op) | (vv, =w — 0) 
+ @, H, cos (v3 t —q,) + D, Hy cos(v, t — qo) , (vy = w@ — 9Q,) 
. z 2 
T, = Hy, sin (vy, t + ¢,) + Hy sin (vy, t + os) (v3 =@ + Q,) “ 
+ 0, Hy sin (v3 t — g,) + 0, Hy sin (vy, t — q,) (v4 = + 2) 


mit neuen Integrationskonstanten H,, H,, y, und Po 
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Es handelt sich offenbar um je vier harmonische Schwingungen mit den vier verschiedenen 
Kreisfrequenzen »y = w + 0; (i = 1,2). In einem Zeigerdiagramm mit T, und T\., als Ko- 
ordinaten, kann die Eigenbewegung (27) der Welle auch gedeutet werden als vektorielle Summe 
von vier Zeigern, die sich mit verschiedenen Winkelgeschwindigkeiten »y um den Ursprung 
drehen. Die GroSe und Lage der Zeiger hangt von den Anfangsbedingungen ab (vgl. Abb. 5). 

Die Eigenbewegung ist im allgemeinen, wie schon aus Abb. 3 ersichtlich, nicht periodisch, 
da die Kreisfrequenzen y nur in Spezialfallen rationale Verhaltnisse aufweisen. 


a) Spezialfalle. 1. Kleine Unrundheit. Die Gleichung (22) liefert die vier Eigenkreis- 
frequenzen 9 im rotierenden, baw. y im ortsfesten System und Gleichung (26) J, und #,: 


0, =o+o, 4=—oO, 9 6 We 
0. =o—o, == s0ON, be 8 w(w@ +o)’ (28) 
C3019 gy =20W+0, oe é a 
= a) =F —\ ¢ 
04 = — 02> Y%=20—oO, 8 w(w— @) | 
{ y clipes 
| 25=W- 
mse 
Ip 
Tr 
es 
A 
V;=O*0, 
Abb. 5. Zeigerdiagramm der stabilen Eigenbewegung Abb. 6. Ortskurve der Zeiger bei stabiler Eigenbewegung fiir m = @/2 
der Welle. und kleine Unrundheit ¢ <1. 


Abb. 3 zeigt die Abhangigkeit »y von w. Es sind dieselben Kurven wie fiir ¢ = 0. Periodische 
Lésungen treten auf, wenn w/@ eine rationale Zahl ist, z. B. ...1/3, 1/2, 1, 2,3,... 


Fur w = @/2 wird (27) zu 


T,, = H, cos (~@t + ¢,) + H, cos (Wt + Gp) —~= H, cos (2wt—q,) + =H, cos (— @p) , 
(29) 


T, = H, sin (—@t + 9) + Hy sin @t + y,) —— H, sin (2 Ot —,) + +H, sin (— gy). 


Es treten Schwingungen mit doppelter und vierfacher Frequenz der Drehzahl auf, wobei die 
,,doppelfrequenten® Schwingungen amplitudenmafig iberwiegen. 

Ein Ortsdiagramm der Gleichung (29) ist in Abb. 6 dargestellt. Die Form der durchlaufenen 
Kurve ist geschlossen und leicht nierenférmig; sie degenieriert fiir ¢ > 0 zu einer Ellipse. 


2. Kreisrunde Welle, ¢ =0. Die allgemeine Lésung (27) laBt sich auf den Fall der 


kreisrunden Welle spezialisieren. Mit ¢ = 0 wird B, = B, = B und w2 =? =@ Die Kreis- 
frequenzen sind (28) zu entnehmen und die Lésung folgt aus (29) mit ¢ = 0: 
T,, = H, cos (— ®t + 9,) + H, cos (wt + @) ; ! 


= (30) 
T, = H, sin (— wt (1) + Hy, sin (@t + 2) . 


Das Zeigerdiagramm besteht aus zwei Zeigern der Lange H bzw. H,, die mit den Winkelge- 
schwindigkeiten —@ bzw. w, also entgegengesetzt drehen. Es handelt sich somit um eine 
ellipsenformige Bewegung mit der Periode 2 7/@ fir z = konstant, was bekannt und zu zeigen 
war. 

6) Anpassung an die Anfangsbedingungen: Es soll in diesem Abschnitt gezeigt werden, 
daB sich die stabile Lésung (26), (27) beliebigen Anfangsbedingungen anpassen JaBt und deshalb 
eine vollstandige Lésung darstellt. In einem beliebigen Zeitpunkt (t = 0 bedeutet keine Ein- 
- schrankung der Allgemeinheit) soll Lage und Geschwindigkeit der Welle als Funktion der Lange z 
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vegeben sein, also 

= v2.0) =Sle)> He 0) = hl), 
yz, 0) =f) » yz, 0) = hy) . 

Mit Bezug auf (25) und (11) ergeben sich die 4. n Gleichungen’ 


fo] 


>, sine T,(0) =f, D> ae? T (0) ao h,, » 


ru nw 
co 


- - na 
>) sin “m2 T,(0) =f, 2 nt T,(0) =h,. 
n= n= 


f Sp oe. WE ; , oe ene & 5 
Multipliziert man beide Seiten mit sin > e dz und integriert von z = 0 bis z = 1, so wird 


I 
oo 


1 
; 2 Aa UE 
>) 7.00) = | fasin de, 70) = Ff hesin ede, 
0 0 


n=1 é 7 
l l 
= 2 _ nse eis 2 . na 
> T,(0) =a Te sin —- dz , > T,(0) = rae sin —-% dz. 
n=1 6 n= 6 


Setzt man an Stelle von JT, und T,, die Integrationskonstanten D; nach den Beziehungen 
(26) oder (27) ein, so wird man fiir jedes n auf vier inhomogene Gleichungen zur Bestimmung 
der vier Kontanten D, gefiihrt. Ist die Determinante dieses Systems von Null verschieden, 
so existieren tatsdchlich nicht triviale Lésungen. Es handelt sich dann um eine vollstandige 
Lésung im stabilen Bereich die den Rand- und Anfangsbedingungen geniigt. 

Im speziellen Falle kleiner Unrundheit ¢ > ¢< 1, e?-> 0 lautet das auf die zweckmaBigere 
Form (26) bezogene Gleichungssystem 


D, D, D; D, | n= 192s Sync | 
! 
2 . Wie 
1+ 3%, 0 I te Wi 0 +f fesin ede 
6 
I 
0 1— #8, 0 1— 3#, +f fysin 2 de (31) 
6 
l 
23 s. = a2 5 SS 
0 @ + 0, (20 + @) 0 —o-+ 0, (20—o) Ff hesin pe 
0 
l 
= fa = 2 Te 
—@+0,(20+ 0) 0 ® + @, (20 —@) 0 4 f hy sin ™2 de 
i) 
mit 
; € w? ae € @? 
a, 8 w(w+ oO)’ Os 8 w(w—@)° 


Die Ausrechnung der Determinante ergibt 

Det = 46, (420). 
Damit ist gezcigt, daB im stabilen Bereich fir alle «@ die Determinante von Null verschieden 
ist, somit mindestens fiir kleine Unrundheit die Konstanten D, nicht trivial werden. 

Als Beispiel sollen fiir den Fall w = @/2 die Anfangshedingungen so gewahlt werden, daB 
fir t = 0 die Welle gerade ist und infolge eines StoBes eine Anfangsgeschwindigkeit f, = v, 
auf ein kurzes Stiick 6) im Abstand z = z, in Richtung x ausgeiibt wird. Weil far t = 0 ruhendes 
und bewegtes Koordinatensystem zusammenfallen, ist die x-Richtung Hauptachsenrichtung fiir 
das gréBere Tragheitsmoment des Querschnittes, — der StoB erfolgt also auf die ,,flache Seite“ 
der Welle. Es gilt 

f.A2) = 0 2 h,(%o) = % (h,(z) a 0, © a 2) ? 
f(z) = 0, hy(z%) = 0. 


* Die hier belanglose Konstante A in (11) wurde eins gesetzt. 
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Die absoluten Glieder der Gleichung (31) lauten dann in der Reihenfolge der Kolonne: 0, 0, 


2 5 LEE : 
T Yo % Sin —— 2% = K,0, und die Konstanten geben 


fem : a K E€ Bs K é 
Di e= 02 Dy =7a(1—): Da 0), Diez! x) 3 
aus Gleichung (26) wird mit (28) endlich 
i (Ce — : — 
(oe =F ltin@s + e(Gsinde—pysin2o)], 
KI = 
Be = Ge (— zoos Be + pyeos 2 we -- z)|: 


In (25) eingesetzt und iiberlagert ergibt mit dem obigen Wert von K fiir w = @/2: 


< B \ 
QQ = wae | Fa ° 


va Deen 5 WES || oe = off te ka. 1 9 

lo 7 Sin oe D)). (32) 
I 1 

5 ie > ee icin te po sin "T= Je 7 cos Mt + 75 beer Ese z) 


In einer Ebene z = konst. durchlauft der Wellenmittelpunkt, wie schon bemerkt, eine ge- 
schlossene, nierenférmige Kurve, die fiir ¢—+ 0 zu einer Geraden in x-Richtung degeneriert. 
Die Komponenten schwingen mit doppelter und vierfacher Frequenz der Drehzahl. 

Beziglich der Oberschwingungen ist zu sagen, daB durch den Sto8 im allgemeinen alle 
Eigenformen (n = 1, 2, 3,...) angeregt werden, wahrend im speziellen fiir einen StoB in Wellen- 
mitte z —1/2 nur symmetrische Formen (n = 1,3,5,...) auftreten. Die Amplituden der 
Oberschwingungen nehmen in jedem Fall mit 1/n? ab. Die Verhaltnisse sind im wesentlichen 


ahnlich denen der runden Welle (€ = 0). 


d) Instabile Lésungen. Die allgemeine Lésung von (13) mit vier Integrationskonstanten 
fiir den instabilen Bereich erfaBt man ebenfalls mit dem Ansatz (14). Die sich ergebenden 
imaginéren Argumente j gt der Kreisfunktionen sind dann, wie anfangs bemerkt, durch Hy- 
perbelfunktionen mit reellem Argument 9 t zu ersetzen. Im iibrigen bleibt die Rechnung die- 
selbe wie fiir den stabilen Fall. Auf die Wiedergabe der Resultate darf deshalb verzichtet werden. 

Hingegen bediirfen die Grenzfalle w =m, und w =w, noch einiger Erlauterungen. Geht 
man mit w =w,, baw. w =w, in die Frequenzgleichung (17) ein, so wird 


ot = 3.08 + ot, f= ott 308, 
Eames (Uae oO=w 
a 
; eo, =0, ae 
und im ortsfesten System 
Fira Uy Oi Cea One ae 
Yo =, — 02 = W Pea 00) 10-0. 
w® =w,:2? a Pe CieaG es ae < if (33) 
¥3 =O, T Q1> Vg = We T Or» 


Vee Oo Op On ¥4 =O, + 02 = %,- 


Infolge des Zusammenfallens zweier Wurzeln liefert der Ansatz (14) nur drei partikulare Inte- 
grale der Gleichung (13) mit den drei Integrationskonstanten H,,y, H. Transformiert ins ro- 
tierende System wird 


T, = H, [cos (x,t. + gy) + 0, cos (v3 t — y)] + H cosa, t,} 
T, = H, [sin (4, t + ¢) + 0, sin (¥3t —@)] + Hsing, t, 


= 0,3 
b __ 20ue1— ej 
eras 2? 
2 @u 1 + Qj (34) 
T,, = H, [cos (y, t + y) + 0, cos (»gt — @)] + H sin (— 
Ee BN oe area coxa 9, 
v* ae 
6,=2 21 @v Qi | 


0? + 21 My 01 
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Es fehlt somit ein weiteres partikulares Integral zur vollstandigen Lésung der Gleichung (13); 
da sie vom 4. Grade in ¢ ist. 
Man erhalt nun mit 
fier he igs sore GG" 
Qs 2 Wu WO =W,: Ms OO, 


LS  @28—o2 ~’ Laat 
tatsichlich das fehlende Integral mit der Konstanten C. Transformiert ins ortsfeste System 


ergibt sich! 


20On 
= C|t cos w, t ——,— sin (—), t 
fie. e u Ope ( u ) 9 
oO =,: 9 : 
, a) 
Dy ==, f sin W, t — ¢ 5— 4 cos (— w, 2) > 
v u 
2 Wy : \ Cy 
== G — a a C08 Mot + t sin (—w, t)], 
v “Ub 
aye nO io 
T,=C E Fer sin w, t + tcos (—a, o| ; 
v u J 


Da hier die Zeit t als Faktor vor einem Summanden auftritt, nimmt diese Teillésung unbe- 
schrankt zu mit der Zeit. Die Bewegung ist instabil und damit gehért auch die vollstandige 
Lésung von (13) fir mw =w, und m =w, zum instabilen Typ. Die Grenzen zwischen stabilem 
und instabilem Bereich sind deshalb zum instabilen Bereich zu zahlen. 


4. Inhomogene Gleichungen (Erzwungene Schwingungen). a) Einflu®B der Unwucht 
(Exzentrizitit des Schwerpunktes). Gesucht wird ein partikulares Integral der inhomogenen 
Gleichung (3) mit dem ,,Stérungsglied“ (absoluten Glied) em? cos 6 baw. ew? sin B, das durch 
eine Verlagerung des Schwerpunktes gegeniiber dem Wellenmittelpunkt, die sogenannte Exzen- 
trizitat e entsteht. 

Ein zeitunabhangiger Ansatz im rotierenden System fihrt die entsprechende Gleichung (3) 
iiber in zwei ungekoppelte, gewéhnliche Differentialgleichungen 
= — wu = e(z) w* cos B(z) , * vlY — wv = e(z) w* sin f(z) . (36) 
Man kann wu und v in Form einer Fourier-Reihe ansetzen, wobei mit Riicksicht auf die gewahlten 
Randbedingungen (8) nur sinus-Glieder zu beriicksichtigen sind: 


= - NUZ = - NUZ 
a > a, sin ——, U = > b,, sin a ae (37) 


n= 


Diese Ausdriicke in (36) eingefiihrt, ergeben mit den Bezeichnungen (12) 


co fe °) 
. naz = nz 
: . ae , , 
a, (w2,, — w?) sin ee a w* cos B , > 5, (w?2,, — w*) sin = ¢ (Ose 
n=1 nw 
Multipliziert man beide Seiten mit sin —— und integriert von 0 bis ] so findet man die zu be- 
stimmenden Koeffizienten: 
l ! 
2 2 7 
a) 2 5 ee jw? 2 NZ 
Oh se eS | ONS /G Shim dz = === = || ein Oo aim dz. 
Oe Ona | P Liveaa, Gey p ine 
0 iy 
Das gesuchte partikulare Integral im rotierenden System lautet dann 
3 , 5 ! 
w . naz . nz 
u =u, = sin Ss e(z) cos f(z) - sin —— dz 
0 
Dy af, —antin FT |e) 008 Ble) sin” dz, 
6 
a 1 (38) 
w* a Tae 4 ; LORS 
v=vU, = > Se a) . e(z) sin B(z) - sin dz. 
0 sow @2,— w? I l | ( ) B( ) I 
6 


' Die Lésungen entsprechen in ihrer Bedeutung denen der Mathieuschen. Die Mathieuschen Funktionen 


sind ebenfalls Losungen der Mathieuschen Gleichung fiir Grenzpunkte zwischen den verschiedenen Gebieten der 
Ince-Siruttschen Karte. 
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oder tansformiert ins ortsfeste System 


x = Uy Coswt + vo sin (—w t), ¥ = Up Sinwt + v9 cos (—wt). (39) 
Es handelt sich in jedem Fall um eine kreisférmige Bewegung der Welle mit der Winkelgeschwin- 
digkeit w im ortsfesten System — wie bei der isotropen Welle (vgl. Abb. 7a). Der Unterschied 
besteht darin, daB erstens der Ausschlag der Welle von der Lage der Unwucht beziiglich dem 
unrunden Querschnitt abhangt und zweitens derjenige Punkt der Welle, der bei der Drehung 
den gréBten Radius beschreibt, der sogenannte ,,hohe Punkt“ im allgemeinen nicht auf dem 


Strahl OS liegt, wie bei der isotropen Welle, was eine gewisse Bedeutung fiir das Auswuchten hat. 


Abb. 7. Erzwungene Unwuchtschwingung. EinfluB der Unwuchtlage. — a) Unwucht beliebig; b) Unwucht auf 1. Hauptschwerachse 


(6 = 0); ec) Unwucht auf 2, Hauptschwerachse (6 = 5) . 


Die erste Tatsache ist unmittelbar aus den folgenden zwei Extremfallen zu ersehen. 


1) Mit e(z) = ey = konst., B(z) = 0 wird aus (38) (vgl. Abb. 7b) 


ies) 


uy =e, wala sin“, x =U, coswt, 
n=1,3, 5 ue 
v0 Vi, Si (it. 
2) Mit e(z) = ey = konst., f(z) = 2/2 wird aus (38) (vgl. Abb. 7c) 
uy, —0; 
~ 3 4A .- WZ x% =v, sin (—wt), 
po eace > na (w2,, — 0%) = Ie a ae 
n=1, 3,5 y =v, cos (— wt). 
Die zweite Tatsache zeigt der Spezialfall 
e = e, = konst, 6 =p, —-konst: 
Aus (38) folgt namlich 
Sie a 2 GOS 
Uo 1 a? 4 TEA 
nes” w2, — w il 


Da der zweite Faktor nur im Falle der isotropen Welle (w?,, = 3,) eins wird, ist die zweite 
Erscheinung offensichtlich (vgl. Abb. 7a). 


Aus Gleichung (38) geht deutlich hervorvor, daB fir w =w,, und m =«,,, die Ausschlage 
gegen oo streben. Der Grund liegt einfach darin, daB bei diesen Drehzahlen die Stérfunktionen 
in ,,Resonanz“ geraten mit den EHigenfunktionen. (Die homogenen Gleichungen enthalten die 
Stérfunktionen als Summanden.) 

Beziiglich der Oberschwingungen ist zu sagen, daB sie im wesentlichen durch Gréfe und 
Lage der Exzentrizitat ¢ langs z beeinfluBt werden. Fiir langs z konstantes é@ nehmen sie mit 
dem reziproken Wert der Ordnungszahlen n ab (gleiche bezogene Verstimmung w/w,, baw. 
w/o,, Vorausgesetzt). 

21 
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b) EinfluB der Schwere. Es ist ein partikulares Integral der inhomogenen Gleichung (3) 
mit dem Stérungsglied g sinjwt bzw. gcoswt zu finden, indem der EinfluB der Erdbeschleu- 
nigung g zufolge der horizontalen Lage der Welle zum Ausdruck kommt. Der Produktansatz (5) 


liefert 


P. 7 5 Bae ee eae 

T, 2w T. 40) oe eT sinwt, (40) 
Ty oo Tu we By ZY aoe Fetoogt | 

Te cae Oa a Pa ee DAU Bs : 


Setzt man wie in (6a) Z!¥/Z = k4 so wird entsprechend den Beziehungen (9), (10) und (12) 


nN Itz 


i [eo 2 o) Tt + T, (a2 —o*)] =gsinwt, 


sin 


aren Paar ai ove 


sin aoe “es +2w ae + T, (@?— w*)] =geoswt, 


nN I % 


Beide Seiten mit sin i 


multipliziert und integriert von 0 bis / ergibt wegen 


fom" Feas (aves 13, nee) 


ie —2w ah Se ACO Os) ae = sin wt, 
us (n= 1, 3, 5,.--) (41) 


. 4. 
T, +20 T, + T, (#2 —o?*) = —S coswt, 


eine gewohnliche inhomogene Differentialgleichung vierter Ordnung in t. 
Der Ansatz fiir ein partikulares Integral von (41) 
i =A sin 6 lb eosin (42) 
fiihrt auf ein inhomogenes Gleichungssystem zur Bestimmung der Konstanten A und B. Die 


Rechnung ergibt 
o ( + Al —4@? 


Aas NO NG 
~ nx (8 —2 0%) (02 20%) —4ot na a|or(1— 7) — 404 | 
" (n = 1,3,5,...). (43) 
@* {1 ——-|— 40? 
ES ee Beas : 
na (wy — 2 w*) (03 —2 0°) —404 nx at|oe(1—T)— 4a | 


Fiir die Transformation ins ortsfeste System ist es zweckmaBig A und B in folgender Weise zu 
zerlegen und von der zweiten Form in (43) Gebrauch zu machen: 


ay 2 
AA, = CE Doe a ae eee 
n n n n NO? _ e2 e 
@? ( hee 4) —402 
= : (44) 
BB = Ce ee pee ss : 
2 a LS GE OR e 
ow? ( — 4) —4Aa? 
Es wird dann nach ausgefithrter Transformation 
To = Di sinwus Tn = G, =D, coa2 ots (45) 


Der anfangliche Produktansatz (5) liefert alsdann mit Z = sin” = die gesuchte Bewegungs- 
gleichung in Form einfacher Reihen: 
to) co) i 


ie.2) 


x == Diz) sin? 0.0% DG Nie > D,, sin ae 
isa SS : 
2 (46) 
Bo ae pls wilt Tee 
y (z) (z) cos2mt, C(z) > C,, sin 7 
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mit 
_ & 4gmit oo 4g mlIt* @* — 4 w? 
n 9 5 =: 1 eR aera: 
(na)? B ae ( eat (nx) Bae ( et ee 


Der sich ergebende Bewegungszustand ist einfach zu tiberschen. Fiir zs — konst. durchlauft der 
Wellenmittelpunkt einen Kreis mit Radius D um einen Mittelpunkt, der yom Ursprung um den 
Betrag C in Richtung y der Schwere verschoben ist. Die Winkelgeschwindigkeit ist die doppelte 
derjenigen der Welle. Man spricht deshalb von Schwingungen mit doppelter Frequenz der Dreh- 
zahl oder kurz von doppelfrequenten Schwingungen. Der Kreis wird bei einmaliger Umdrehung 


der Welle zweimal durchlaufen. Vel. Abb. 8. 


Es werden nur symmetrische Schwingungsformen angeregt, 
da die Summation in (46) nur tiber n = 1,3,5,... zu erstrecken 
ist. Fir die bestimmten Drehzahlen 

oS w2 w? _ @ e 

a Vi; (w2 +02) 2 yj 4 0) 
wird der Nenner der zeitunabhangigen Konstanten A, B, C und 
D Null. Die Stérfunktionen geraten bei diesen Drehzahlen in 
Resonanz mit einer Eigenfunktion. Man erkennt dies rasch im 
Falle einer kleinen Unrundheit ¢ > ¢<1,¢? + 0. Es ist fiir diesen 
Fall die kritische Winkelgeschwindigkeit der Welle fiir die doppel- 
frequente Schwingung nach (47) @ = @/2, also (mit Riicksicht auf 
die Naherung <? ~ 0) ungefahr die Halfte derjenigen einer un- 


runden Welle mittlerer Biegefestigkeit B—=(B,-+ B,)/2. Das 
partikulare Integral (45) der inhomogenen Gleichung ist daher fiir “Pb. 8. Erzwungene Schwingung infolge 

=e i 2 a 3 der Schwere bei horizontaler Wellenlage. 
® = w/2 in der allgemeinen Lésung (29) der homogenen Gleichung — — y hat die Richtung der Erdbeschleu- 


enthalten, und somit ist die Resonanzbedingung erfillt. ee 


Die Bedeutung der Mittelpunktsverschiebung C(z) erkennt man wiederum rasch im Falle der 
kleinen Unrundheit ¢ + ¢< 1, c2 + 0. Wie aus (46) hervorgeht ist fiir diesen Fall C(z) unab- 
hangig von der Drehzahl und betragt 


_mgl 4 —y 1 . naz 
C(z) = 3 vee sin —— 


Dies ist eine der méglichen Formen fiir die Gleichung der elastischen Linie der Welle infolge des 
Eigengewichtes!. Fiir geringe Unrundheit beschreibt also der Wellenmittelpunkt angenahert 
einen Kreis mit dem Radius D(z) um die statische Gleichgewichtslage einer Welle mittlerer Biege- 
steifigkeit B. 

Die Lésung in Form der Summe von Oberschwingungen gestattet sofort das Verhalten 
der Oberschwingungen zu beurteilen. Aus Gleichung (46) erkennt man unschwer, daf die 
Oberschwingungen (gleiche bezogene Verstimmung w/@® vorausgesetzt) im Verhaltnis 1/n® 
(n = 1,3,5,...) abnehmen, wenn n die Ordnung angibt. 


Im Gegensatz zu den durch eine Unwucht oder Schwerpunktsexzentrizitat bedingten Ober- 
schwingungen, welche fiir konstante Exzentrizitat @ im Verhaltnis 1/n(n = 1, 2, 3,...) ab- 
nehmen, diirften somit die doppelfrequenten Oberschwingungen von ungleich geringerer Be- 
deutung auf das Schwingungsverhalten eines Rotors sein. 


c) Lagerkrafte. Durch die erzwungenen Schwingungen werden radial gerichtete Krafte 
auf die Lager wirksam. Ihre Berechnung kann getrennt werden fiir diejenigen, welche durch 
Unwucht oder Schwerpunktsexzentrizitat, und fiir diejenigen, welche durch die Schwere her- 
vorgerufen werden. Die beiden Anteile sind in ihrer Wirkung zu tiberlagern. 


1 Vel. S. Timoshenko, Strength of Materials Bd. II, S. 47. Fiir z = 1/2 wird unter Beriicksichtigung der 
ersten drei Glieder 


l —2 ie, 
C() = 1302-10 Sy eee 


Palle 
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1) EinfluB der Unwucht. Die Summe beider Lagerkrafte im rotierenden System kann mit 
Gleichung (1) berechnet werden. Es ist allgemein 


l i ! 1 
Kp =| Fodz — Balad, K,=—[-F dv— B, [a dz. (48) 
0 0 0 0 

Mit wu und v aus (38) wird 
1 


l 

0° > p 4 

ow 2 d 5. aE 5 YE 

Ae ie == > ee - |sin —— | e(z) cos B(z) sin —— dza| dz 

K, Kuo B, { w2,—o? | | dz l (2) Ble) l 4 
ee 


0 
l l 


sd s = w? 2 Final gees 2 : . WIZ 
Ky e=-530-= Bae Sai] fa sin | e(z) sin A(z) sin ie dz| dz 
Fal on 


0 0 


oder einfacher mit Riicksicht darauf, daB das zweite Integral von z unabhiangig ist, 
bs l 
Aa? Bu aw 3 GBH 
ae 3 
Kio = — a > Ceres e(z) cos f(z) sin j dz, 
6 
! 


47° B wo? ; 5 UE 
K, 0 =a zt > [es [ ee sin f(z) sin a. 


Se 
n (00) } 
0 


Es ist hier zweckmaBig, die Komponenten der Exzentrizitat é in sinus-Reihen zu entwickeln. Mit 


e(z) cos 6(z) = a, + DEC sin “i ; e(z) sin f(z) = by + DA sin “_ 


iu n= 


und wegen der Orthogonalitatsrelation 


wird schlieBlich im rotierenden System 


2igte Big wo? 4 
Rios ree > are emer Fe a, + “,). 
eines | (49) 
K, 0 ser > i. 2 bo + by] 
ate o2,— wo \n a 


co 


dy By 8 x? ni Oe 
Oe cana ea Lee saan 
SIR, un @ (e cos hb = a), 
3 b B82 << w? (esinf = by). 
Ky,» =—2 2 eee re 0 
ae! et eT ol oe 
n=1,3,5 NG 


Im ortsfesten System sind die Lagerkrafte zeitlich veranderlich, entsprechend der Transfor- 
mation (2) 

K, = Kyo cosmt + Kyo sin (— wt), Ky, = Kyo sinwt + Kyo cos (— wt). (50) 
Resonanzen treten auf bei den kritischen Drehzahlen 


O=Wyn> O=O,,° 


Im Gegensatz zu den Auslenkungen macht sich der EinfluB der Oberschwingungen auf die 
Lagerkrafte starker bemerkbar. Er ist verstandlicherweise abhangig von GréBe und Lage der 
Exzentrizitat e(z) langs z. Fiir ein langs z konstantes @ wachst der Einflu® mit n2 (gleiche Ver- 
stimmung (/,,, baw. w/w,,, vorausgesetzt) 
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2) EinfluB der Schwere. Die Summe beider Lagerkrafte im rotierenden System wird durch 


Gleichung (48) und Gleichung (5) mit Z = sin a wie folgt bestimmt: 
| 
a aie u= S'sin"**. 7 


l ) 
a : . nwzZ 
K, = B, | o™ dz, ‘= y er peak ee 
0 
Nach durchgefiihrter Differentiation und Integration ist 


Le =2 B, > eal (ee i, =2 ae oe Cr) Loe ; 
n=153, 5 TSS. 


Mit Hilfe der Transformationsformeln 


K, = K, coswt— K,sinwt, Ne = K, sinwt-+ K,coswt 
sowie den Beziehungen (42) und (44), namlich 
T, =(¢, 4+ D,) sinot, T, =(C€, — D,) coswt 
lassen sich die Lagerkrafte im ortsfesten System berechnen. Nach einigen Umformungen entsteht 
Dein 2 ot, K, = C— Deos2ot, (51) 


wobei 


D— > Gel PAC (Bean ae (bee By. 


a RS 


c= > Cr) (Gn (B, + B,) —D, (B, — Bl 


n=—1,3, 5 


oder mit den GréBen B, @ und @ entsprechend (19) ausgedriickt 


n2 
H—a35.5 


= 8gm ely tl 4 w? — = B 
ore wees » wie 2 ; O = @, = (nm)? |/ a» 
sme Sa wal | hat 
(52) 
c bam “1 oe 2 oo 
gh 2 


oo aT es 62 
a1 —7)— 408 *a(1—F)—4o 


Fir kleine Unrundheit ¢ > ¢ <1, e*?-> 0 wird aus (52) ,wenn mit G = g mI das Eigengewicht 
der Welle bezeichnet wird!, 
4 Sy 1 4 w? See) pee 
D=<c, > ( ); C=¢5 >) 4=6. (53) 


n2\ @? — 4 w? TU J 1 
n=1, 3,5 n=1, 3,5 


Man erkennt aus der Form (51), dafB es sich um einer Schwingung mit doppelter Frequenz der 
Drehzahl handelt. Im Falle kleiner Unrundheit beschreibt der Schwingungszeiger einen Kreis 


mit Radius D um den Vektor G. Die Resonanzen sind durch Gleichung (51) bestimmt. Insbe- 
sondere ist die kritische Drehzahl fiir kleine Unrundheit w % @/2. 

Gleiche Verstimmung w/@ vorausgesetzt, wachsen die Oberschwingungen im Verhaltnis 
1/n?. Ihr Einflu8 auf niedere Drehzahlen ist deshalb bedeutend geringer als etwa im Falle einer 
(langs z konstanten) Exzentrizitat des Schwerpunktes. Ein Vergleich 1a8t den Schlu8 zu, dab, 
gleiche Verstimmung vorausgesetzt, die Oberschwingungen der doppelfrequenten Schwingungen 
in ihrer Wirkung n* mal kleiner sind als etwa die Schwingungen mit Drehzahlfrequenz infolge 
einer tiber die Lange z gleichmaBig verteilten Unwucht. Ihre technische Bedeutung ist deshalb 
ungleich geringer. 


1 Fiir eine Lésung in geschlossener Form siehe W. Kellenberger, Brown Boveri Mitteilungen 42 (1955) S. 79. 
Der Vergleich ergibt 


4 wl 4 w?/@? 1 /coshL—1 cosL—1 e 
a ot (A a : se ES, Laon Vos 
Ge > (aa orat | 2L( sinh L sin L a a? 
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5. Zusammenfassung. Es wird das Biegeschwingungsverhalten einer glatten, rotierenden 
Welle untersucht, die einen unrunden Querschnitt aufweist. Das Problem wiirde mathematisch 
die Auflésung eines Systems von drei linearen, partiellen Differentialgleichungen sechster Ord- 
nung mit harmonischen Koeffizienten erfordern. Es gelingt nun aber mit den tblichen Voraus- 
setzungen und durch Formulierung der Gleichgewichtsbedingungen in rotierenden Koordinaten 
ein System von zwei Differentialgleichungen vierter Ordnung mit konstanten Koeffizienten zu 
erhalten. Zusammen mit den vorgeschriebenen Rand- und Anfangsbedingungen ist damit ein 
Randwertproblem achter Ordnung definiert. Die vollstandige Lésung fiir vereinfachte Randbe- 
dingungen wird in Form von Reihen angegeben. Auf die gegeniiber einer runden Welle sich er- 
gebenden Besonderheiten wird hingewiesen. Insbesondere werden die Eigenschwingungen dreh- 
zahlabhangig. Es ergeben sich instabile Bereiche und neue, vermehrte Resonanzlagen. Die Un- 
wuchtschwingung wird abhangig von der Lage der Unwucht beziiglich des Querschnittes, und 
nicht zuletzt, bei horizontaler Lage der Welle werden Schwingungen mit doppelter Frequenz der 
Drehzahl angeregt, sog. ,,doppelfrequente Schwingungen“ mit eigenen, neuen kritischen Dreh- 
zahlen. Der Sonderfall geringer Unrundheit, wie er in der Technik haufig ist, vereinfacht die 
Resultate erheblich. 


(Eingegangen am 19. Dezember 1957.) 
Anschrift des Verfassers: Dr. Walter Kellenberger, Wettingen (Aargau) (Schweiz), Flurweg 14. 
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Inhaltstbersicht: Elastostatik: Ein Elastizitatsgesetz fiir kleine Verzerrungen. Der ebene Spannungs- 
zustand und das Biegeproblem. Drehsymmetrische Spannungszustande. Torsion zylindrischer Stabe mit 
beliebigem Querschnitt. Biegung diinner Platten. — Schwingungslehre: A. Schwingungen mit einem 
Freiheitsgrad: I. Autonome Bewegungen: Allgemeine Untersuchungen. Konservative Schwingungen. 
Gedampfte Schwingungen. Selbsterregte Schwingungen. — II. Heteronome Bewegungen: Zwangsschwin- 
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Die meisten Probleme, die man in der klassischen Mechanik mathematisch durch eine lineare Differential- 
gleichung zu formulieren pflegt, sind streng genommen nichtlinear. Beriicksichtigt man dies und geht von 
den nichtlinearen Differentialgleichungen aus, so lassen sich zahlreiche, von der linearen Theorie nicht 
erfaBbare Erscheinungen beschreiben, die heute in der Technik eine wichtige Rolle spielen. Dies gilt be- 
sonders von der Elastostatik und der Schwingungslehre. 

Im ersten Teil des Buches wird eine Elastizitatstheorie der kleinen Verzerrungen entwickelt, bei der an 
die Stelle des Hookeschen Gesetzes ein nichtlineares Spannungs-Dehnungsgesetz tritt, das den tatsach- 
lichen Eigenschaften vieler Werkstoffe besser gerecht wird. Die Theorie wird auf viele, von der klassischen 
Elastomechanik her bekannte Probleme angewandt. Der zweite Teil befaBt sich mit den wichtigsten Ver- 
fahren und Ergebnissen der in den letzten Jahren schon sehr weitgehend entwickelten Theorie der nicht- 
linearen mechanischen Schwingungen. 

Besonderer Wert wurde gelegt auf die Darstellung der mathematischen Methoden, um den Leser instand 
zu setzen, auch neue Anwendungsbeispiele, die im Buch noch nicht behandelt worden sind, selbstandig zu 
bearbeiten. 


Das Buch setzt die Kenntnis der Grundlagen der Mechanik und der héheren Mechanik voraus und ist vor 


allem fiir Mathematiker, Physiker und Ingenieure bestimmt, die mit diesen neuen Anwendungsgebieten der 


Mechanik vertraut werden wollen. 
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der magnetischen Spannungen und der Konzentrationsspannungen. Die spannungsfreien Strukturkriim- 
mungen. Die Grenzflachenbedingungen fiir die Distorsionen. Die Grenzflachenbedingungen fiir die Defor- 
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methode zur Berechnung der Selbstenergie singularer Versetzungen. Fremdatome als elastische Dipole 
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Die Versetzung hat als Trager der plastischen Verformung und als elementare Eigenspannungsquelle in den 
letzten Jahren gréStes Interesse gewonnen. Im Gegensatz zu den bisherigen zusammenfassenden Dar- 
stellungen, die einen atomistischen Standpunkt beziehen, wird hier die Stellung der Versetzung in der 
Kontinuumsmechanik besprochen, in der sie eine Rolle spielt, die in vieler Hinsicht derjenigen des Stromes 
in der Elektrodynamik analog ist. Der Leser wird in die neuen Begriffe eingefiihrt. Der plastische und 
elastische Anteil der gesamten Verformung werden beide durch asymmetrische Tensorfelder beschrieben, 
deren Rotation die Versetzungsdichte ist. Letzteres bedingt im allgemeinen die Nichterfiillung der Kom- 
patibilitatshedingungen von de St. Venant bei Anwesenheit von Versetzungen. Die spezifischen Hilfa- 
mittel von Eigenspannungszustinden sind der Inkompatibilitatstensor und der Spannungsfunktions- 
tensor. Die hiermit arbeitende neue Theorie der Eigenspannungsbestimmung wird konsequent dargestellt, 
Die Kontiniuumstheorie wird auf die realen Kristalle angewandt. Die punktférmigen Gitterfehlstellen 
werden als elastische Singularitaten behandelt, fiir deren Wechselwirkung mit den elastischen Feldern im 
Kristall sehr weitreichende, einfache Formeln gelten. Ferner wird die nicht-Riemannsche Geometrie der 
Versetzungen von Kondo und Bilby und Mitarbeitern kurz eae die eine schéne Anwendung des 
Cartanschen Begriffes der Torsion auf den realen Korper ist. 
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